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A  COMPUTER  PROGRAM  TO  RAPIDLY  SCREEN 
EXPLOSION-PROOF  ENCLOSURES 


By  Frank  T.  Duda1 


ABSTRACT 

A  U.S.  Bureau  of  Mines  interactive  FORTRAN  computer  program  can  screen  explosion-proof 
enclosures  intended  for  underground  service.  The  program  incorporates  important  information  and 
criteria  to  assist  enclosure  manufacturers  and  Mine  Safety  and  Health  Administration  (MSHA)  personnel 
by  providing  a  first-order  approximation  of  rectangular  explosion-proof  enclosure  characteristics  prior 
to  actual  certification  testing. 

The  program  uses  information  and  data  that  can  be  retrieved  automatically  during  program  execution. 
The  program  also  performs  the  calculations  for  a  series  solution  algorithm  and  various  equations  to 
evaluate  the  enclosure  for  both  a  static  internal  pressure  and  an  internal  explosion.  The  program  writes 
the  intermediate  and  final  calculations  and  information  concerning  the  status  of  the  screening  process 
to  a  summary  file  on  a  floppy  disk. 

One  data  base  contains  MSHA  certification  requirements  given  in  the  Code  of  Federal  Regulations 
(30  CFR  18).  A  second  data  base  contains  material  properties  of  acceptable  construction  materials  and 
weld  joint  efficiencies.  A  series  algorithm  determines  bolt  and  cover  stresses  for  a  simply  supported 
rectangular  cover  with  arbitrary  bolt  locations.  Limit  and  yield  line  analysis  is  used  to  calculate  the 
strength  of  body  panels. 


'Electrical  engineer,  Pittsburgh  Research  Center,  U.S.  Bureau  of  Mines,  Pittsburgh,  PA. 


INTRODUCTION 


To  receive  certification,  an  explosion-proof  (XP)  en- 
closure must  meet  the  requirements  of  part  18,  title  30  of 
the  Code  of  Federal  Regulations  (CFR)  and  other  policy 
guidelines  of  the  Mine  Safety  and  Health  Administration 
(MSHA).3  The  MSHA  investigator  first  visually  checks  the 
enclosure  to  make  sure  that  it  is  properly  constructed  from 
acceptable  materials.  If  the  enclosure  meets  the  visual 
requirements,  it  is  then  put  through  a  lengthy  and  costly 
series  of  inspections  and  explosion  tests.  Finally,  the 
applicant  must  perform  static  pressure  tests  on  enclosures 
of  a  specific  design.  The  test  procedures  must  be 
submitted  to  MSHA  for  approval. 

A  methodology  to  prescreen  enclosures  before  sub- 
mitting them  to  the  actual  certification  process  was  de- 
veloped by  the  U.S.  Bureau  of  Mines  under  contract  with 
the  Southwest  Research  Institute.4  The  method  addresses 
primarily  the  strength  requirements.  Southwest  Research 
Institute  had  previously  evaluated  welding  procedures, 
construction  techniques,  and  the  most  common  materials 
of  construction  used  in  the  manufacture  of  XP  enclosures.5 
It  had  also  performed  numerous  tests  and  analyses  on 
enclosures  including 

1.  Finite-element  analysis  to  determine  stresses,  strains, 
and  deflections  in  the  enclosure  produced  by  internal  static 
pressures. 

2.  Analysis  to  predict  the  response  of  enclosures  to 
dynamic  loading. 

3.  Hydrostatic  testing  of  enclosures  with  measurements 
of  strains,  deflections,  and  strain  patterns. 

4.  An  analysis  of  stresses  during  impact  and  thermal 
shock  tests  for  lenses  (30  CFR  18.66). 


2U.S.  Code  of  Federal  Regulations.  Title  30-Mineral  Resources; 
Chapter  I-Mine  Safety  and  Health  Administration,  Department  of 
Labor;  Subchapter  D-Electrical  Equipment,  Lamps,  Methane  Detectors; 
Tests  for  Permissibility;  Fees;  Part  18-Electric  Motor-Driven  Mine 
Equipment  and  Accessories;  July  1,  1985. 

3Mitchell  D.  W.  (Chief,  Approval  and  Certification  Center,  MSHA). 
Letter  to  All  Interested  Parties,  Apr.  22,  1977;  available  upon  request 
from  F.  T.  Duda,  BuMines,  Pittsburgh,  PA. 

4Cox,  P.  A.  Estimating  Rectangular  XP  Enclosure  Performance, 
Final  Engineering  Report.   SW  Res.  Inst.  Project  06-8188,  1984,  35  pp. 

5Cox,  P.  A.,  O.  H.  Burnside,  E.  D.  Esparza,  F.  D.  Lin,  and  R  E. 
White.  A  Study  of  Explosion-Proof  Enclosures  (contract  H0377052, 
SW  Res.  Inst.).  BuMines  OFR  96-83,  Dec.  1982,  426  pp.;  NTIS  PB  83- 
205450. 


5.  An  experimental  study  of  the  response  of  XP  en- 
closures to  vibrations. 

These  surveys,  tests,  and  analyses  contributed  directly  to 
the  development  and  verification  of  a  methodology  for  the 
rapid  approximate  evaluation  of  rectangular  XP  enclosure 
performance. 

The  Bureau  adapted  that  methodology  into  an  inter- 
active computer  package  written  in  Microsoft6  FORTRAN. 
The  package  is  designed  to  be  run  on  an  IBM  personal 
computer  (PC)  or  compatible  PC  having  two  floppy  disk 
drives  and  a  256-Kb  memory.  It  can  also  be  run  on  a  PC 
equipped  with  a  hard  disk  drive,  one  floppy  disk  drive,  and 
a  256-Kb  memory.  An  8087  coprocessor  will  shorten 
execution  time  but  is  not  mandatory.  The  package  is 
intended  to  be  used  by  both  the  manufacturers  of  XP 
enclosures  and  MSHA  to  provide  a  first- order 
approximation  of  rectangular  XP  enclosure  characteristics 
prior  to  actual  certification  testing.  Enclosures  that  are 
clearly  inadequate  can  be  redesigned  by  the  manufacturer 
before  expensive  approval  testing  is  begun.  The  computer 
program  also  formalizes  a  procedure  currently  being 
followed  by  MSHA  and  provides  the  inspector  with  a 
checklist  to  ensure  that  no  evaluation  criteria  are 
overlooked. 

The  first  part  of  this  report  describes  the  computer 
program,  data  bases,  and  system  requirements.  Also  in- 
cluded is  the  procedure  for  using  the  program.  The  sec- 
ond part  of  the  report  contains  an  example  of  the  strength- 
screening  aspect  of  the  program  to  determine  deflections 
and  stresses  on  a  10.75-  by  13.5-in  bolted  cover  resulting 
from  an  internal  static  pressure  of  100  psi.  The  input  data 
and  output  results  are  presented.  The  results  of  calcu- 
lations performed  on  a  32-bit  mainframe  computer  using 
a  classical  finite-element  solution  and  results  of  the 
Bureau-developed  program  are  compared  to  ensure  the  va- 
lidity of  the  results.  Finally,  the  results  of  using  the 
program  to  determine  deflections  on  other  panels  and 
covers  are  presented. 


6Reference  to  specific  products  does  not  imply  endorsement  by  the 
U.S.  Bureau  of  Mines. 


DISCLAIMER  OF  LIABILITY 


The  U.S.  Bureau  of  Mines  expressly  declares  that  there 
are  no  warranties  express  or  implied  which  apply  to  the 
software.  By  acceptance  and  use  of  said  software,  which 
is  conveyed  to  the  user  without  consideration  by  the  U.S. 
Bureau  of  Mines,  the  user  thereof  expressly  waives  any 


and  all  claims  for  damage  and/or  suits  for  or  by  reason  of 
personal  injury,  or  property  damage,  including  special, 
consequential  or  other  similar  damages  arising  out  of  or 
in  any  way  connected  with  the  use  of  the  software. 
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DESCRIPTION  OF  EXPLOSION-PROOF  ENCLOSURE  SCREENING  PROGRAM 


PACKAGE  OVERVIEW  AND  SYSTEM 
REQUIREMENTS 

The  XP  enclosure  screening  package  consists  of  two, 
double-sided,  double- density  IBM  PC  or  compatible  5-1/4- 
in  floppy  diskettes  and  this  report.  Copies  of  the  two 
diskettes  can  obtained  from 

Pittsburgh  Research  Center 
U.S.  Bureau  of  Mines 

P.O.  Box  18070 
Pittsburgh,  PA   15236 

Requests  for  the  package  should  reference  this  report  and 
include  two,  double-sided,  double-density,  5-1/4-in  floppy 
diskettes.  The  user  is  expected  to  have  available  a  blank, 
formatted,  5-1/4-in  floppy  diskette. 

This  report  is  designed  to  serve  as  a  user's  manual.  It 
contains  illustrations  and  tables  that  the  user  must  refer  to 
when  using  sections  of  the  screening  program.  The 
screening  program  is  designed  for  use  on  an  IBM  PC  or 
compatible  PC  equipped  with  dual  double-sided,  5-1/4-in 
floppy  disk  drives  and  a  minimum  of  256-Kb  of  random 
access  memory  (RAM).  The  program  can  also  be  run  on 
a  PC  equipped  with  one  floppy  disk  drive,  a  hard  disk 
drive,  and  a  256-Kb  memory,  if  appropriate  changes  are 
made  in  the  loading  procedure.  An  8087  coprocessing  chip 
is  desirable  but  not  mandatory.  The  time  required  for 
calculations  will  be  decreased  significantly  if  an  8087  is 
used.  The  XP  program  is  designed  to  run  as  an  executable 
task  after  the  Microsoft  disk  operating  system  (MS-DOS), 
any  version,  is  loaded.  The  ANSI  device  driver  must  be 
installed  at  the  time  that  the  MS-DOS  operating  system  is 
loaded. 

The  first  diskette,  marked  "XP  Enclosure  Screening 
Package,"  contains  the  executable  image  of  the  interactive 
XP  program  called  XP.EXE  and  all  the  support  files  that 
make   up  the  data  bases.      This  is  the   only  diskette 


necessary  to  run  the  program.  The  second  diskette, 
marked  "XP  Screening  Package  Source  Code,"  contains  the 
FORTRAN  source  code  for  the  main  module  program 
called  XPMAIN.FOR,  the  FORTRAN  source  code  for  all 
the  necessary  subroutines  called  by  XPMAIN.FOR,  and  a 
copy  of  this  report  as  a  DOS  text  file.  Both  diskettes 
should  be  write  protected. 

The  information  for  the  data  bases,  the  equations  and 
sequence  of  calculations,  and  the  graphs  and  figures  were 
taken  directly  from  the  work  cited  in  footnote  4.  In  most 
cases  the  data  from  the  graphs  and  figures  have  been  in- 
corporated in  the  data  base  or  the  computer  code.  Where 
such  incorporation  was  not  practical  or  possible,  the  user 
is  directed  to  the  tables  or  figures  in  this  report. 

Appendix  A  contains  a  catalog  of  all  the  files  on  the 
two  diskettes.  Appendix  B  contains  listings  of  the 
FORTRAN  source  code  for  XPMAIN.FOR  and  all  the 
called  subroutines. 

OPERATING  PROCEDURE 

In  order  to  run  the  screening  program  on  a  PC 
equipped  with  two  floppy  disk  drives,  MS-DOS  must  first 
be  loaded  in  the  PC.  After  the  correct  day  and  time  is 
entered,  the  diskette  marked  "XP  Enclosure  Screening 
Program"  should  be  installed  in  disk  drive  A  and  a  blank 
formatted  diskette  should  be  installed  in  disk  drive  B.  In 
response  to  the  A:  prompt,  type  the  letters  XP  and  then 
press  the  return  key.  The  screening  program  will  auto- 
matically load,  and  the  following  title  display  will  appear  in 
reverse  video: 

XP  ENCLOSURE  SCREENING  PROGRAM 

V5.0  08/11/88 

WRITTEN  BY:   F.  T.  DUDA 

If  the  screen  display  is  not  working  properly,  the  ANSI 
device  driver  probably  was  not  installed  when  MS-DOS 


was  loaded.  If  for  some  reason  the  ANSI  device  driver  is 
not  installed,  the  program  will  print  an  escape  sequence 
to  the  screen  instead  of  affecting  the  cursor  movement. 
The  file  ANSI.SYS  and  the  configuration  file 
CONFIG.SYS  must  be  on  the  diskette  that  is  used  to  load 
MS-DOS.  Also,  the  file  CONFIG.SYS  must  contain  the 
statement:   DEVICE = ANSI.SYS. 

If  the  screen  display  is  not  working  properly,  the  XP 
screen  program  should  be  terminated  by  striking  the 
<CNTRL>  and  C  keys  at  the  same  time.  MS-DOS 
should  be  reloaded,  making  sure  that  the  ANSI  device 
driver  is  correctly  installed. 

In  order  to  run  the  screening  program  on  a  PC 
equipped  with  a  hard  disk  drive  and  one  floppy  disk  drive, 
MS-DOS  must  first  be  loaded  in  the  PC.  The  file 
ANSI.SYS  and  the  configuration  file  CONFIG.SYS  must 
be  in  the  same  directory  that  contains  the  MS-DOS  op- 
erating system.  All  the  files  from  the  diskette  marked  "XP 
Enclosure  Screening  Program"  should  be  copied  to  the 
hard  disk.  A  blank  formatted  diskette  should  be  installed 
in  the  floppy  disk  drive.  In  order  to  run  the  program,  at 
the  C:  prompt  type  the  letters  XP  and  then  press  the 
return  key. 

INTERACTIVE  MAIN  MODULE-XPMAIN 

The  computer  program  described  here  is  an  interactive, 
menu  driven  program  offering  the  user  choices  for  per- 
forming a  number  of  different  tasks.  It  prompts  the  user 
for  input  data  as  well  as  appropriate  choices  from  the  main 
menu  at  every  level  of  the  program.  Users  have  the  option 
of  only  performing  one  task  or  being  directed  through  a 
logical  sequence  of  the  screening  process  to  ensure  that 
they  have  not  overlooked  an  important  task. 

The  program  was  written  following  a  methodology  de- 
veloped through  a  U.S.  Bureau  of  Mines  contract.  The 
method  addresses,  primarily,  enclosure  strength.  The  tasks 
available  to  the  user  are  the  following: 

1.  Screen  for  flame  path  requirements. 

2.  Screen  for  weld  classification. 

3.  Screen  for  acceptable  construction  materials. 

4.  Screen  for  enclosure  and  window  strength. 

5.  Screen  for  ruggedness. 

Intermediate  and  final  results  of  calculations  as  well  as 
information  about  the  status  of  the  enclosure  are  written 
to  a  file  called  RESULT.DOC  on  the  B  drive  diskette. 

After  the  program  is  loaded,  an  introduction  describing 
the  scope  of  the  package  and  the  operation  instructions 
regarding  the  use  of  the  package  are  displayed.  Before  the 
first  menu  is  displayed,  the  user  is  asked  for  the  iden- 
tification of  the  enclosure.  This  can  be  any  label  up  to  64 
characters  in  length.  This  identification  is  used  to  ref- 
erence the  enclosure  throughout  the  screening  process. 
If  users  want  to  screen  another  enclosure,  they  must  exit 
the  XP  program  by  choosing  the  appropriate  menu  item 
and  then  execute  XP  again.   Each  time  that  the  program 


XP  is  executed,  the  summary  file  RESULT.DOC  is  written 
over.  To  save  the  summary  file  from  a  screening  process, 
use  the  MS-DOS  command  RENAME  and  enter 
RESULT.DOC  NEWFILE.DOC  at  the  prompt  before  XP 
is  executed  again. 

At  this  point,  users  are  asked  if  they  want  to  perform  a 
single  screening  task  or  if  they  want  to  be  directed  through 
a  logical  sequence  of  the  total  screening  process. 

SCREEN  FOR  FLAME  PATH 
REQUIREMENTS  MODULE 

This  feature  of  the  package  provides  the  users  with  the 
capability  to  view  the  applicable  information  for  their  en- 
closures concerning  flame  path  requirements  as  specified 
in  30  CFR  18.31.  After  viewing  the  applicable  flame  path 
requirements,  the  user  is  asked  to  check  the  enclosure 
being  screened  against  the  clearance  requirements.  The 
user  should  measure  the  appropriate  widths  and  clearances 
and  also  determine  them  from  the  enclosure  drawings.  If 
the  enclosure  does  not  meet  the  requirements,  the 
information  is  written  to  the  file  RESULT.DOC. 

Flame  path  requirements  are  specified  for  the  following 
three  volume  ranges  for  empty  enclosures:  less  than  45  cu 
in,  45  through  123  cu  in,  more  than  124  cu  in.  The 
program  will  ask  the  user  to  either  enter  the  internal 
volume  or  to  provide  the  dimensions  of  the  enclosure 
panels  so  that  the  internal  volume  can  be  calculated.  The 
program  calculates  the  internal  volume  by  first  determining 
the  inside  dimensions  of  the  enclosure  by  subtracting  the 
thicknesses  of  the  top  face  and  bottom  face  from  the  front 
face  and  one  of  the  side  faces,  and  by  subtracting  two 
times  the  thickness  of  one  side  panel  from  the  front  face. 
The  internal  volume  is  then  determined  by  multiplying  the 
remaining  area  of  the  front  face  by  the  effective 
nonadjacent  length  of  the  side  face.  This  procedure  only 
works  for  a  rectangular  or  square  enclosure.  The  user  is 
given  the  option  to  either  accept  the  result  of  this 
calculation  or  to  enter  the  volume. 

Flame  path  requirements  specify  both  the  joint  width 
and  joint  clearance.  The  appropriate  requirements  keyed 
to  one  of  the  three  volume  ranges  have  been  entered  into 
a  data  base  and  can  be  retrieved  by  the  user  and  viewed 
on  the  display.  There  are  footnotes  for  some  of  the  re- 
quirements which  can  be  viewed  also.  Users  are  also 
provided  with  a  procedure  that  they  should  follow  for 
other  then  metal-to-metal  joints  such  as  gasketed  surfaces. 
Because  it  takes  several  screens  to  view  the  data  base  for 
any  one  of  the  volume  ranges,  the  user  can  freeze  the 
screen  at  an  appropriate  section  or  review  any  previous 
screen. 

If  the  enclosure  joints  satisfy  30  CFR  18  requirements, 
it  is  unlikely  that  the  requirements  will  be  exceeded  after 
an  internal  explosion.  The  bolted  covers  are  checked  in 
the  task  4,  screen  for  enclosure  strength,  to  determine 
whether  or  not  permanent  bolt  elongation  will  occur, 
which  could  alter  the  flame  path. 


SCREEN  FOR  ACCEPTABLE  MATERIALS 
OF  CONSTRUCTION  MODULE 

Material  Properties  Data  Base 

A  data  base  of  acceptable  construction  materials  and 
their  corresponding  mechanical  properties  are  utilized  by 
this  task.  There  are  four  categories  of  materials  that  make 
up  the  data  base.  One  lists  acceptable  steels,  another  lists 
acceptable  aluminums,  the  third  lists  glasses  and  plastics, 
and  the  fourth  lists  adhesives  and  sealants.  The  material 
properties  of  yield  strength,  elastic  modulus,  Poisson's 
ratio,  and  density  are  stored  in  the  data  base  keyed  to  the 
material.  The  user  is  asked  to  choose  one  of  the  four  ma- 
terials. If  the  data  base  for  either  the  acceptable  steels  or 
aluminums  is  chosen,  the  list  of  materials  and  their  prop- 
erties is  shown.  The  user  is  then  asked  to  indicate  which 
material  is  used  in  the  construction  of  the  enclosure.  If 
the  material  is  listed,  the  enclosure  is  rated  acceptable  for 
this  part  of  the  screening.  Appropriate  material  properties 
for  the  steel,  aluminum,  and  glass  and  plastic  listings,  such 
as  the  yield  strength,  can  be  retrieved  and  used  if  required 
by  another  task.  The  data  bases  for  adhesives  and  sealants 
can  be  viewed  also,  although  the  material  properties 
cannot  be  retrieved. 


Seven  of  the  most  commonly  used  steels,  either  wrought 
or  cast,  and  their  physical  properties  of  yield  strength, 
elastic  modulus,  Poisson's  ratio,  and  density  are  stored  in 
the  STEELDES.DAT  file.  The  acceptable  steels  and  then- 
material  properties  are  shown  in  table  1. 

Seven  of  the  most  commonly  used  aluminums,  either 
wrought  or  cast,  and  their  physical  properties  of  yield 
strength,  elastic  modulus,  Poisson's  ratio,  and  density  are 
stored  in  the  ALDES.DAT  files.  The  seven  acceptable 
aluminums  and  their  material  properties  are  shown  in 
table  2. 

Four  of  the  most  commonly  used  glasses  and  plastics 
and  their  physical  properties  of  yield  strength,  elastic 
modulus,  Poisson's  ratio,  and  density  are  included  in  the 
GLASSDES.DAT  file.  The  acceptable  glasses  and  plastics 
and  their  material  properties  are  shown  in  table  3. 

Two  materials  that  are  commonly  used  as  both  ad- 
hesives and  sealants  in  the  construction  of  explosion-proof 
enclosures  are  Hysol  934  bonding  epoxy  and  GE108  RTV 
elastomeric  adhesive  and  sealant.  Their  physical  prop- 
erties of  yield  strength,  elastic  modulus,  Poisson's  ratio, 
and  density  are  stored  in  a  file  called  SEALDES.DAT  and 
are  shown  in  table  4. 


Table  1  .-Characteristics  of  steels  commonly  used  in  explosion-proof  enclosure  design1 


Yield 
strength, 
kips/in2 


Elastic 

modulus, 

kips/in2 


Poisson's 
ratio 


Density, 
lb/iny 


A-36 

AISI  1025  low-carbon  steel   

SA-516,  grade  60 

SA-573,  grade  65 

AISI  Type  304  annealed  stainless  steel 
AISI  410  chromium  stainless  steel  .  .  .  . 
ASTM  A48  cast  iron,  grade  30 

STEELDES.DAT  disk  file. 


36 

29,000 

0.33 

0.286 

55 

29,000 

.33 

.284 

32 

29,000 

.33 

.286 

35 

29,000 

.33 

.286 

35 

29,000 

.33 

.286 

35 

29,000 

.33 

.286 

30 

29,000 

.33 

.260 

Table  2.-Characteristics  of  aluminums  commonly  used  in  explosion-proof  enclosure  design1 


Yield 
strength, 
kips/in2 


Elastic 

modulus, 

kips/in2 


Poisson's 
ratio 


Density, 
lb/in7 


Aluminum  alloy  2024  (T351) 

5083-H32  (H116)    

5456-H111    

6061-T6  (T651) 

A201.0  T4  cast  aluminum2   . 

A319.0T63 

A413.0  cast  aluminum4  .  .  .  . 

1  ALDES.DAT  disk  file, 
formerly  ASTM  .CQ51A 
formerly  ASTM  SC64D. 
formerly  ASTM  413.0. 


42 

10,700 

0.34 

0.100 

30 

10,200 

.33 

.096 

33 

10,200 

.33 

.096 

40 

10,700 

.34 

.100 

31 

10,300 

.33 

.101 

24 

10,700 

.33 

.100 

19 

10,300 

.33 

.101 

Table  3. -Characteristics  of  glasses  and  resins  commonly  used  In  explosion-proof  enclosure  design1 


Yield  strength,  kips/in2, 
at  various  temperatures 

Elastic  modulus,  kips/in2, 
at  various  temperatures 

Poisson's 
ratio 

Density, 
lb/inr 

68°  F 

150°  F    240°  F    300°  F 

68°  F     150°  F    240°  F    300°  F 

Glass: 

Soda-lime   NA 

Borosilicate   ....                NA 

Resin: 

Lexan  101   9.0 

Merlon  3113    ...                9.2 

NA         NA          6.6 
NA         NA          6.3 

7.9         6.4           NA 
8.0         6.6           NA 

NA        NA          NA        10,000 
NA       NA          NA         9,100 

345      304         240        NA 
335       302         220        NA 

0.24 
.25 

.25 

.25 

0.0914 
.0806 

.0433 
.0433 

NA        Not  available. 
'GLASSDES.DAT  disk  file. 

Table  4. -Characteristics  of  adhesives  and  sealants  commonly  used  in  explosion 

-proof  enclosure  design1 

Yield  strength,  kips/in2, 
at  various  temperatures 

Elastic  modulus,  kips/in2, 
at  various  temperatures 

Poisson's 
ratio 

Density, 
lb/inr 

68°  F 

150°  F    240°  F    300°  F 

68°  F    150°  F    240°  F    300°  F 

Hysol  934  epoxy    ...           10.0 
GE108  RTV    .05 

8.0         6.9           NA 
NA        6.3           NA 

300            NA         NA           100 
.190     NA         NA          NA 

0.49 
.49 

0.040 
.039 

NA       Not  available. 
^LASSDES.DAT  disk  file. 

Changing  the  Material  Properties  Data  Base 

If  a  material  and  its  properties  are  not  found,  but  are 
judged  to  be  more  relevant  than  other  materials  listed  in 
the  data  base,  the  appropriate  file  can  be  changed  to  in- 
clude the  material  and  its  properties  using  an  editor  such 
as  the  MS-DOS  EDLIN  editor  or  a  word  processor  with 
ASCII  text  file  capability.  The  user  must  exit  the  XP 
screening  program  in  order  to  change  the  data  base. 

The  number  of  entries  in  a  particular  data  base  must 
remain  fixed.  For  example,  the  data  base  for  acceptable 
steels,  STEELDES.DAT,  must  always  have  seven  types  of 
materials.  In  addition,  the  name  of  the  material  and  its 
mechanical  properties  must  always  occupy  the  same  col- 
umns in  the  data  base.  The  best  method  to  add  a  new 
material  is  to  edit  or  write  over  one  that  is  not  relevant. 
It  is  also  possible  to  add  a  new  line  or  lines  to  the  file 
providing  that  an  equal  number  are  deleted.  Every  time 
a  new  entry  is  added  to  the  data  base,  one  must  be  deleted 
to  keep  the  total  number  equal  to  the  original  number. 
The  user  should  note  that  some  of  the  entries  require  two 
lines.  To  change  one  of  those  entries,  both  lines  must  be 
replaced  by  an  entry  requiring  two  lines.  It  is  always 
possible  to  add  a  blank  line  if  required. 

WELD  JOINT  CLASSIFICATION  MODULE 

It  is  necessary  to  correctly  classify  all  weld  joints 
according  to  the  American  Welding  Society  welding  code 
DW14.4  in  order  to  assure  that  the  joints  will  exhibit 
sufficient  strength  to  withstand  static  loads  and  dynamic 
loading  produced  by  an  internal  explosion. 

The  main  purposes  of  this  section  are  to 

1.  Provide  the  user  with  a  formalized  procedure  for  the 
correct  classification  of  the  welded  joints. 


2.  Determine  whether  any  class  VI  welds  are  used. 

3.  Determine  the  static  and  dynamic  joint  efficiencies 
based  on  the  weld  classifications. 

Welds  are  classified  in  DW14.4  according  to  the  joint 
configurations  and  inspection  requirements.  The  user 
should  classify  the  welds  either  from  taking  the  infor- 
mation from  the  drawings  or  from  visual  examination  of 
the  welded  joints  of  the  actual  enclosure.  The  user  must 
first  determine  the  class  by  the  configuration  of  the  welded 
joint  as  shown  in  figure  1.  (The  weld  configurations  could 
not  be  accurately  displayed  by  the  computer  program.) 
The  user  then  must  make  sure  that  the  enclosure  weld 
joints  were  properly  inspected.  Inspection  requirements 
for  the  different  classes  are  displayed  for  the  user.  It 
should  be  noted  that  classes  I,  II,  and  III  require 
ultrasonic,  radiographic,  or  dye  penetrant  inspections, 
which  are  seldom  performed  by  enclosure  manufacturers. 
Therefore,  the  vast  majority  of  enclosures  will  fall  within 
classes  IV  and  V.  Finally,  the  inspector  must  check  the 
visual  and  dimensional  requirements  of  the  welds  as  given 
in  table  5.  It  was  not  possible  to  display  efficiently  all  the 
allowable  structural  discontinuities  by  class;  therefore,  the 
user  is  directed  to  use  table  5. 

The  welds  must  meet  all  the  requirements  of  the  class 
in  order  to  be  so  classified.  For  example,  if  a  weld  has  a 
joint  configuration  of  class  I  as  determined  from  figure  1, 
but  only  meets  visual  and  dimensional  requirements  of 
class  III,  determined  from  table  5,  the  weld  must  be  con- 
sidered a  class  III  weld.  The  program  determines  the 
highest  classification  of  the  three  categories  and  applies 
that  classification  to  the  weld  joint.  In  addition,  the  user 
is  asked  if  any  class  VI  joints  are  used  in  the  enclosure.  A 
class  VI  joint  is  unacceptable  for  dynamic  loading.  The 
results  of  the  classification  procedure  are  written  to  the 
RESULT.DOC  file.   If  any  class  VI  joints  were  used,  the 


enclosure  would  fail  this  section  of  the  screening  process. 
The  static  and  dynamic  efficiencies  of  the  weld  joint  that 
are  used  in  the  body  and  panel  strength  section  are 
determined   by  weld   class.      The   static   and   dynamic 


efficiencies  of  the  weld  joints  are  stored  by  class  in  the 
WELDEF.DOC  file.  They  are  displayed  to  the  user  by 
class  and  can  be  extracted  from  the  data  base  and  used  if 
required  by  another  task. 


Table  5.-Allowable  structural  discontinuities 


AWS 
weld 
class1 


Extent  of  visual  weld  defects  allowed 


Undercut 


Inadequate 
joint  penetration 


Surface  holes 


Weld  bead 
irregularities 


IV 


Visually  free  from  undercut  in 
direction  transverse  to  primary 
stress,  0.01-in  max  depth, 
smoothly  contoured  undercut 
parallel  to  direction  of 
primary  stress. 


Undercut  0.01-in  max  depth 
when  direction  transverse 
to  primary  stress.    1/32-in 
max  depth  when  direction 
parallel  to  primary  stress. 

.  .   do    


Undercut  1/32-in  max  depth 
when  direction  transverse  to 
primary  stress.   1/4  thickness 
of  member  or  1/16  in  max 
depth  when  direction  parallel 
to  primary  stress. 


Visually  free  from  all 
evidence  of  inadequate 
joint  penetration  or 
types  of  fusion  defects 
at  root  of  weld. 


1/3  thickness  of  member 
or  1/2-in  max  length 
inadequate  penetration. 
Sum  of  all  defects  not 
to  exceed  1  in  per  12  in. 

2/3  thickness  of  member 
or  3/4-in  max  length 
inadequate  penetration. 
Sum  of  all  defects  not 
to  exceed  1  in  per  6  in. 

Thickness  of  member  or 
1-in  max  length  inadequate 
penetration.   Sum  of  all 
defects  not  to  exceed  1  in 
per  6  in. 


Sum  of  1/32-in  max  diam 
scattered  or  randomly 
distributed  surface 
holes  not  to  exceed 
1/8  in/in2  or  1/4  in 
per  12-in  length. 
Linear  orientation  of 
holes. 

Sum  of  1/16-in  max  diam 
surface  holes  not  to 
exceed  1/4  in/in2  or 
1/2  in  per  12  in. 


Sum  of  3/32-in  max  diam 
surface  holes  not  to 
exceed  3/8  in/in2  or 
3/4  in  per  12  in. 


No  visual  inspection 
specified. 


Surface  of  weld  to  be 
ground  or  machined 
to  produce  smooth, 
uniform  profile, 
blending  with  base 
metal,  without 
notches,  depressions 
or  sharp  edges. 

Rough  irregular  welds 
and  excess  rein- 
forcement to  be 
ground  smooth. 


No  visual  inspection 
specified. 


Do. 


VI 


do 


No  visual  inspection 
specified. 


do 


Do. 


AWS     American  Welding  Society, 
^o  cracks  allowed. 
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Figure  1  .-Classification  of  welded  joints. 


SCREEN  FOR  ENCLOSURE  STRENGTH  MODULE 

The  enclosure  strength  can  be  evaluated  for  both  a 
static  internal  pressure  and  an  internal  explosion.  The 
procedure  treats  only  pressure  loads;  thermal  loads  are 
neglected. 

This  task  can  be  used  to  check  the  following  different 
items  of  the  enclosure  strength: 

1.  Penetrations. 

2.  Covers,  including  the  bolts  that  seal  the  cover  to  the 
enclosure  body. 

3.  Body  panels. 

4.  Windows  and  lenses. 

All  penetrations,  including  windows,  are  included  in 
item  1  regardless  of  where  they  are  located  in  the  en- 
closure. The  construction  of  the  reinforcements  are 
checked  to  determine  if  the  reinforcements  contribute  to 
any  reduction  in  strength  of  the  component  containing  the 
reinforcement.    If  there  is  a  reduction  in  strength,  it  is 


considered  when  the  appropriate  component  is  screened 
in  item  3. 

Analysis  of  the  cover  in  item  2  includes  the  bolts  that 
secure  the  cover  to  the  enclosure  body.  The  analysis  of 
the  bolted  cover  is  based  on  a  linear  elastic  analysis. 
However,  the  cover  should  also  be  analyzed  as  a  body 
panel  in  item  3  if  significant  plasticity  is  indicated. 

Body  panels  refer  to  all  panels  in  the  remainder  of  the 
enclosure,  including  the  enclosure  bottom.  The  procedure 
of  item  3  is  set  up  so  that  as  many  panels  as  needed  can 
be  analyzed.  Some  irregular  enclosures  might  require 
several  panels  to  treat  even  one  side  of  the  enclosure. 

The  windows  and  lenses  are  checked  for  strength  to 
react  to  the  internal  pressure.  The  lip  that  supports  the 
window  or  lens  is  checked  also. 

Penetrations 

When  penetrations  for  windows,  cables,  switches,  etc., 
are  placed  through  the  structural  boundary,  reinforcement 
must  be  added  around  the  penetration  if  stiffness  and 


strength  are  not  to  be  reduced.  This  section  of  the 
screening  process  calculates  strength  and  stiffness  factors 
for  body  panels  that  have  circular  penetrations  for  win- 
dows, cables,  switches,  etc.  A  calculated  strength  factor 
less  than  1  for  a  penetration  near  the  center  line  of  a  body 
panel,  is  an  indication  that  the  penetration  has  weakened 
the  structure.  For  this  case,  the  calculated  strength  factors 
are  used  in  item  2  in  the  body  panel  screening  process. 

The  strength  and  stiffness  factors  are  calculated  using 
the  elastic  modulus  and  yield  strength  of  the  plate  and 
reinforcement  and  the  geometry  of  the  reinforcement. 
Users  can  retrieve  the  appropriate  elastic  modulus  and 
yield  strength  from  the  acceptable  materials  data  base  or 
they  can  enter  them.  The  dimensions  and  location  of  the 
penetration  must  be  entered  according  to  the  geometry 
shown  in  figure  2.  The  user  must  refer  to  figure  2  in  this 
report  since  it  was  not  possible  to  include  the  information 
in  the  screening  program.  The  strength  and  stiffness  fac- 
tors are  computed  for  two  orthogonal  planes  that  are  par- 
allel to  the  panel  sides  which  pass  through  the  penetration 
at  the  maximum  opening.  The  strength  factors  STRFA 
and  STRFB  are  used  in  item  2  in  the  body  panel  screening 
task  if  they  meet  the  following  two  conditions:  They  are 
less  than  1  and  the  penetration  lies  near  the  center  line. 

As  an  example,  if  STRFB  of  the  plate  in  figure  2  was 
less  than  1,  it  would  be  used  in  the  analysis  of  the  plate 
because  it  lies  near  the  center  line,  which  is  parallel  to  side 
B;  however,  STRFA  would  not  be  used  because  it  is  too 
far  removed  from  the  center  line,  which  is  parallel  to  side 
A.  The  calculated  values  for  STRFA  and  STRFB  are  writ- 
ten to  the  RESULT.DOC  file  and  are  also  passed  on  to 
the  body  panel  screening  task  if  appropriate. 

The  stiffness  factors  STFFA  and  STFFB  are  calculated 
in  this  procedure  also  but  they  are  not  used  directly  in  any 
other  task  of  the  package.  However,  the  program  writes 
this  information  to  the  RESULT.DOC  file  if  they  are  less 
than  1  so  that  the  user  can  be  aware  of  a  potential 
concern.  A  factor  less  than  1  is  not  cause  to  reject  or  fail 
the  enclosure.  For  more  complex  geometries,  the  in- 
spector must  estimate  an  equivalent  thickness,  t,  or 
calculate  the  area  moments  of  inertia  of  the  before  and 
after  configurations.  Dividing  the  after  value  by  the  before 
value  gives  STFFA  and  STFFB. 

Equations  for  the  strength  reduction  factors  are  based 
upon  the  plastic  sections  modulus.  The  plastic  section 
modulus  is  used  because  the  analysis  of  enclosure  strength 
is  based  upon  yield  line  analysis  methods.  As  for  stiffness, 
the  strength  factors  for  more  complex  reinforcement  must 
be  calculated  or  must  be  estimated  by  using  an 
equivalent  t. 

Bolt  and  Cover  Strength 


Stresses  in  the  bolted  cover  can  be  calculated  for  any 
static  pressure  using  a  previously  developed  series  solution 
algorithm.  This  algorithm  can  be  used  to  calculate  bolt 
and  cover  stresses  for  rectangular  covers  of  arbitrary  size 
and  for  arbitrary  bolt  spacing.  The  following  are 
assumptions  upon  which  the  procedure  is  based: 

1.  Linear  elastic  behavior. 

2.  Cover  bolted  to  a  rigid  structure. 

3.  Bolt  heads  rotate  with  the  cover  plate. 

4.  No  shear  deformations  in  the  bolts  or  cover. 

The  procedure  is  used  to  calculate  the  cover  deflections, 
rotations,  moments,  shears,  and  stresses,  and  the  moments 
and  stresses  of  the  bolts.  The  screening  program  needs 
the  following  information: 

1.  Dimensions  of  the  cover. 

2.  Physical  properties  of  the  material  of  construction. 

3.  Pressure  for  which  the  analysis  is  to  be  made. 

4.  Number  of  bolts  and  their  locations. 

5.  Physical  properties  of  the  bolts. 

6.  Dimensions  of  the  bolts  and  their  engagements. 

The  user  must  enter  the  dimensions  of  the  cover  and 
bolt  locations  as  described  in  figure  3.  The  bolt  locations 
are  given  in  x,y  coordinates  with  the  0,0  location  in  the 
bottom  left  of  the  cover.  Users  can  either  retrieve  the 
required  physical  parameters  of  the  cover  and  bolts  from 
the  data  bases  or  they  can  enter  them. 
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Because  the  behavior  of  the  bolts  and  cover  are  inter- 
dependent, an  elastic  analysis  is  used  to  give  a  coupled 
solution  for  stresses  in  the  bolts  and  cover.  If  cover 
stresses  exceed  the  material  yield  stress,  the  cover  should 
then  be  analyzed  as  a  panel. 


SECTION    A-A 

Figure  2.  Geometry  of  circular  penetrations,  a,  b,  h,  length, 
width,  and  thickness,  respectively,  of  panel  containing 
penetration;  D0,  diameter  of  penetration,  including  reinforcement; 
t,  thickness  of  reinforcement;  H,  buildup  of  reinforcement. 
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KEY 
-$-Bolt  locations 


PLAN  VIEW 


Figure  3. -Geometry  of  bolted  cover.   A,  B,  length  and  width, 
respectively,  of  cover  panel;  A  should  always  be  longer  than  B. 


Users  are  given  the  option  of  having  the  screening  pro- 
gram determine  the  minimum  length  of  thread  engage- 
ment for  the  bolts  used.  If  they  choose  this  option,  they 
must  obtain  bolt  data  from  outside  references  and/or  from 
the  enclosure  drawing.  The  calculated  value  of  the  mini- 
mum length  of  thread  engagement  should  be  less  than  or 
equal  to  the  length  of  thread  engagement  specified  on  the 
drawing  or  measured  from  the  enclosure.  The  calculations 
are  based  on  strength  requirements  only.  If  they  are  in 
conflict  with  the  requirements  of  30  CFR  18.31,  the  users 
must  use  their  judgment  in  passing  or  failing  the  enclosure. 

The  program  uses  the  series  solution  algorithm  to  cal- 
culate the  strength  quantities  for  the  cover  at  the  center 
of  each  side  and  the  middle  of  the  panel,  the  locations 
where  the  maximum  deflections  are  likely  to  occur.  The 
program  also  calculates  the  moments  and  stresses  for  each 
bolt.  The  algorithm  includes  an  infinitely  squared  set  of 
linear  algebraic  equations.  The  deflections  and  stresses 
are  determined  from  the  solutions  to  a  truncated  set  of 
those  equations.  In  order  to  implement  the  algorithm  in 
Microsoft  FORTRAN  and  then  integrate  that  code  as  part 
of  the  screening  program,  it  was  necessary  to  truncate  the 
infinite  series  to  three  terms.  A  Gauss-Seidel  iterative 
procedure  was  used  to  solve  the  resulting  set  of  nine 
simultaneous  linear  equations. 

Stresses  in  the  bolts  and  cover  are  based  upon  elastic 
material  behavior.  For  this  type  of  behavior,  bolt  axial 
stresses  should  not  exceed  the  bolt  yield  strength,  and  the 
combination  of  bending  and  axial  stress  should  not  exceed 
the  bolt  ultimate  strength.  Using  these  criteria  will  assure 
that  bolt  stretch  will  not  occur,  which  could  increase  the 
fiange-to-cover  flame  gap.  Cover  stresses  can  exceed  the 
material  yield  stress.  If  they  do,  the  user  is  encouraged  to 
check  the  cover  as  a  panel  using  the  body  panel  strength 
analysis  procedure. 


Body  and  Panel  Strength 

The  ability  of  an  enclosure  to  withstand  an  internal 
static  pressure  of  150  psig  and  pressure  loading  produced 
by  an  internal  explosion  is  predicted  by  using  the  results  of 
calculations  based  on  limit  or  yield  line  analysis  techniques. 
The  capability  of  an  individual  panel  to  withstand  a  static 
pressure  is  determined  by  performing  a  strength  calcu- 
lation using  the  panel  geometry,  the  panel  material  prop- 
erties, and  the  panel  joint  configuration.  The  calculation 
is  repeated  for  each  body  panel.  If  the  strength  calculation 
is  less  than  150  psig  for  any  body  panel,  the  prediction  is 
that  the  enclosure  may  not  be  able  to  withstand  the  actual 
static  pressure  test.  For  the  internal  explosion  evaluation, 
a  modified  strength  calculation  is  used  to  account  for  dy- 
namic effects.  If  the  modified  strength  calculation  is 
greater  than  70  pet  of  the  static  pressure  calculation  for 
any  individual  panel,  the  prediction  is  that  the  enclosure 
will  not  pass  the  actual  explosion  test  as  specified  in  30 
CFR  18.62.  The  70-pct  factor  is  used  to  account  for  the 
fact  that  permanent  deformations  in  the  panel  cannot 
exceed  0.040  in/ft. 

Internal  Static  Pressure 

An  allowable  static  pressure  is  computed  for  the  indi- 
vidual panels  in  the  enclosure  to  determine  whether  or  not 
they  will  withstand  a  150-psig  static  pressure  as  required  by 
30  CFR  18.  This  calculation  is  based  on  a  limit  analysis 
that  permits  plastic  hinges  to  form  in  the  structure  and 
residual  deformations  to  develop.  If  the  calculated  allow- 
able static  pressure  for  all  the  panels  is  greater  than  150 
psig,  the  enclosure  is  deemed  acceptable  for  the  static 
pressure  loading.  Intermediate  and  final  results  of  the 
calculations  are  written  to  the  RESULT.DOC  file. 

The  allowable  static  pressure,  PA,  is  given  by 

PA  =  RM/AB,  (1) 

where  RM    =  panel  resistance, 

A    =  panel  width, 

and  B    =  panel  length. 

The  panel  resistance  is  calculated  from  the  panel  yield 
strength,  panel  geometry,  joint  efficiency,  and  strength 
factors.  The  panel  yield  strength  can  be  retrieved  from  the 
acceptable  materials  data  base  or  entered.  The  panel 
length,  width,  and  thickness  must  be  entered.  Dimension 
B,  the  length  of  the  panel,  must  be  the  same  or  longer 
than  dimension  A,  the  width  of  the  panel.  The  joint  effi- 
ciencies are  designated  as  EJA  for  the  joint  efficiency  of 
the  joint  parallel  to  side  A,  and  EJB  as  the  joint  efficiency 
of  the  joint  parallel  to  side  B.  The  efficiencies  for  dif- 
ferent types  of  joints  can  be  retrieved  from  data  bases.  If 
the  joints  are  welded,  the  efficiencies  can  be  retrieved 
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according  to  class  from  the  weld  joint  efficiency  data  base. 
If  a  bolted  cover  is  being  analyzed  as  a  panel,  appropriate 
joint  efficiencies  are  calculated.  If  penetrations  are  pre- 
sent, the  strength  reduction  factors  STRFA  and  STRFB 
calculated  by  item  1  of  the  "Screen  for  Enclosure  Strength 
Module"  section  are  used.  If  no  penetrations  are  present, 
or  they  do  not  contribute  to  a  reduction  in  strength  if 
present,  a  value  of  1.0  is  used. 

Internal  Explosion 

To  determine  whether  or  not  a  panel  is  adequate  for 
dynamic  loading,  and  equivalent  static  pressure,  PEQ,  is 
first  calculated  and  then  compared  to  the  previously  de- 
termined allowable  static  pressure,  PA,  from  the  internal 
static  pressure  section.  This  section  should  be  performed 
immediately  after  the  internal  static  pressure  check.  The 
allowable  static  pressure  is  reduced  by  a  factor  of  0.7  to 
limit  permanent  deformations  in  the  enclosure  to  values 
approximately  equal  to  0.04  in/ft.  If  PEQ  is  less  than  0.7 
PA,  the  enclosure  panel  is  considered  adequate  for  the 
dynamic  loading. 

In  order  to  calculate  the  equivalent  static  pressure, 
PEQ,  a  peak  explosion  pressure  and  minimum  pressure 
rise  time  must  be  first  estimated.  A  value  of  100  psig  is 
used  as  a  peak  explosion  pressure  in  the  algorithm.  An 
estimate  for  the  minimum  pressure  rise  time  (TRMIN)  in 
milliseconds  is  estimated  by  the  equation 

(TRMIN)  =  (PMAX/4.77)    *   (^VOLUME  ).    (2) 

The  user  must  enter  the  volume  of  the  enclosure  in 
cubic  meters  or  use  the  volume  if  it  had  been  previously 
determined  during  the  flame  path  screening  section. 

A  dynamic  load  factor,  DLF,  must  then  be  determined 
for  the  panel.  It  is  necessary  to  first  determine  a  funda- 
mental circular  frequency  for  the  rectangular  panel.  The 
circular  frequency  is  a  function  of  the  boundary  conditions 
on  the  panel  and  its  physical  properties.  The  physical 
properties  and  geometry  can  be  either  retrieved  from  the 
data  base  or  entered.  The  length,  B,  of  the  panel  must  be 
greater  than  the  width,  A.  Possible  boundary  configu- 
rations are  shown  in  figure  4.  The  user  must  determine 
the  configuration  number  from  figure  4  for  the  particular 
boundary  conditions  for  the  plate  of  interest.  The  fun- 
damental period,  T,  is  calculated  from  the  circular  fre- 
quency. A  ratio  of  the  minimum  rise  time  (TRMIN)  to 
the  fundamental  period,  T,  is  then  calculated.  The  data 
from  figure  5  are  used  by  the  computer  program  to  cal- 
culate the  DLF  given  the  ratio  TRMIN/T.  PEQ  is  then 
determined  by 


PEQ  =  (100  psig)*DLF. 


(3) 


Window  and  Lens  Strength 

Headlight  lenses  are  required  to  be  at  least  the  equiv- 
alent of  1/2-in  Pyrex  glass.  The  user  is  first  asked  to 
check  all  glass  lenses  to  make  sure  that  the  thickness 
requirement  is  met. 

To  check  the  strength  of  the  windows  and  lenses,  a  pro- 
cedure is  used  that  first  determines  a  bending  stress.  This 
calculated  stress  is  then  compared  to  the  yield  strength  of 
the  material  used  for  the  window  or  lens.  If  the  material 
is  glass,  the  calculated  bending  stress  must  be  less  then 
one-half  the  yield  strength  for  the  window  to  be  accept- 
able. If  the  material  is  plastic,  the  calculated  bending 
stress  must  only  be  less  than  the  yield  strength  of  the 
plastic. 

The  bending  stress  is  calculated  from  the  dimensions  of 
the  window,  an  assumed  applied  pressure  of  150  psig,  and 
a  stress  constant,  K.  The  procedure  is  used  for  both 
rectangular  and  circular  windows.  If  the  window  is  rect- 
angular, the  stress  constant,  K,  is  determined  from  either 
figure  6  or  figure  7. 

Figure  6  is  used  if  the  edges  of  the  window  are  free  and 
figure  7  is  used  if  the  edges  are  clamped.  Dimension  A, 
the  length  of  the  window,  must  be  longer  than  the  width, 
B.  If  the  window  is  circular,  the  stress  constant,  K,  is 
determined  by  the  type  of  support  used  for  the  window. 

A  final  check  done  in  this  section  is  to  determine  if  the 
lip  that  holds  the  window  can  withstand  the  internal  pres- 
sure. An  allowable  shear  force,  VA,  which  produces  a 
plastic  hinge  at  the  built-in  end  of  the  lip,  is  first  cal- 
culated. The  yield  strength  of  the  lip  material  is  needed  by 
the  algorithm.  The  user  can  either  retrieve  it  from  the 
acceptable  materials  data  base  or  enter  it.  The  user  must 
enter  the  geometry  as  defined  in  figure  8.  The  allowable 
shear  force,  VA,  is  compared  to  an  estimated  shear  force, 
VMAX,  due  to  loading  on  the  lip  resulting  from  a  150-psig 
static  pressure.  The  strength  of  the  lip  is  considered 
acceptable  if  VA  >VMAX. 


Configurations  and  No. 
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A  summary  of  intermediate  and  final  results  for 
TRMIN,  T,  PEQ,  and  DLF  are  written  to  the 
RESULT.DOC  file. 


Figure  4. -Boundary  conditions  for  rectangular  plates. 
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TRMIN/T 

Figure  5. -Response  of  one  degree  of  movement  elastic  system. 


Enclosure 
panel 


Figure  8.-Window  and  lip  geometry.  I,  length  of  lip  material; 
tL,  thickness  of  lip;  c,  length  of  window  supported  by  lip;  P, 
uniform  pressure  loading  over  length  of  window. 
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Figure  6. -Stress  factor  for  free  window. 
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Figure  7. -Stress  factor  for  clamped  window. 


SCREEN  FOR  RUGGEDNESS  MODULE 

The  ruggedness  of  the  enclosure  is  checked  by  first 
calculating  the  kinetic  energy  that  can  be  absorbed  by  the 
enclosure.  This  kinetic  energy  is  then  compared  to  a 
proposed  ruggedness  criteria  based  on  kinetic  energy7. 
Because  a  ruggedness  requirement  does  not  now  exist  for 
XP  enclosures,  the  tested  enclosure  should  not  be  rejected 
if  it  does  not  meet  the  proposed  kinetic  energy  ruggedness 
criteria.  The  fact  that  it  does  or  does  not  meet  the 
proposed  criteria  is  recorded  in  the  RESULT.DOC  file. 

The  static  resistance,  RM,  of  the  cover  and  side  panels 
as  calculated  in  the  section  to  screen  the  panels  for 
strength  is  used  to  determine  an  estimate  of  the  kinetic 
energy  that  can  be  absorbed  by  the  enclosure.  The  kinetic 
energy  that  can  be  absorbed  by  a  panel  will  be  equal  to 
the  strain  energy  that  the  panel  absorbs  in  plastic 
deformation.  An  allowable  permanent  deformation  of  one 
plate  thickness,  h,  is  used  to  compute  the  absorbed  strain 
energy.  The  kinetic  energy  that  can  be  absorbed  by  the 
enclosure  cover  or  panel  is  calculated  by 


KE  =  RM  *  h 


(4) 


This  value  is  compared  to  the   appropriate  following 
criteria: 


Surfaces  exposed 
to  roof  fall  .... 

Surfaces  exposed 
to  side  impact    . 

Protected  surfaces 


KE  -  3,700  ft  •  lb. 

KE  =  8,000  ft  -lb. 
KE  =     800  ft -lb. 


7Work  cited  in  footnote  5. 
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EXAMPLE  CALCULATIONS  FOR  BOLTED  COVER 


A  simulation  for  the  cover  and  12  bolts  shown  in  figure 
9  was  performed  using  the  enclosure  strength  task  of  the 
program.  The  input  data  for  this  cover  and  bolt  configu- 
ration are  stored  in  the  program  data  base  and  can  be 
easily  accessed  through  the  menu.  The  results  obtained 
by  running  the  program  on  the  user's  machine  using  the 
data  from  the  program  data  base  can  be  compared  to  the 
results  shown  in  figure  10.  The  dimensions  and  physical 
parameters  used  in  the  simulation  are  given  in  table  6. 

The  following  method  was  used  to  run  the  simulation: 

1.  From  the  main  menu,  the  screen  for  enclosure 
and/or  window  strength  option  was  chosen. 

2.  From  the  next  menu,  the  check  cover  for  strength 
option  was  chosen. 

3.  Cover  dimensions  were  entered.  Dimension  A  was 
entered  as  13.5,  B  as  10.75,  and  the  cover  thickness,  H, 
was  entered  as  0.5. 


4.  In  order  to  use  the  12-bolt  configuration  of  figure  9, 
the  second  option,  use  the  12-bolt  configuration  described 
on  the  following  screen,  was  chosen. 

5.  The  elastic  modulus,  yield  strength,  and  Poisson's 
ratio  for  the  cover  were  entered  (they  could  have  been 
accessed  from  the  acceptable  materials  data  base).  For 
this  simulation,  the  elastic  modulus  was  entered  as  10,000 
kips/in2,  the  yield  strength  was  entered  as  33  kips/in2,  and 
Poisson's  ratio  as  0.3. 

6.  From  the  next  menu,  the  second  option,  use  the 
values  for  1/2-in  class  five  bolts,  was  chosen. 

7.  The  results  of  the  simulation  and  a  record  of  the 
input  data  were  written  to  the  RESULT.DOC  file  in 
drive  B. 


Bolt  3 


Simply  supported 
edges 

Bolt  2 


Bolt 
locations 


Inside  edge  of 
flange 


Simply  supported 
edges 


Figure  9.-Bolted  cover  used  in  sample  simulation. 
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THE  COVER  DIMENSIONS  ARE: 

A:        13.5000000      B:        10.7500000 
THICKNESS:   5.000000E-001 

DIMENSION  OF  AA  -  MATRIX  -  NN  =    3 

PLATE  DIMENSIONS  -  A  =   13.5000,  B  =   10.7500,  H  =    .5000 

INTERNAL  PRESSURE  -  PO  =     100.0000 

PLATE  PROPERTIES 

E  =  10000000.00,  MU  =   .3000,  YP=     33000.00 

BOLT  PROPERTIES 

E  =  30000000.00,  I  =      .003068,  G  =  11540000.00,  KS  =   1.3330 

A  =      .196300,  L  =    .8500,  D  =    .5000 


PLATE  DEFLECTIONS  AND  STRESSES 

X  =    6.7500,  Y  =    5.3750 

DEFLECTION  =   .01807 

SIGMA-X  =    8649.5200,  SIGMA-Y  =   10969.5900 

TAU-XY  =       .0000,  TAU-XZ  =      -.0003,  TAU-YZ  =      -.0043 

X  =    6.7500,  Y  =     .0000 

DEFLECTION  =  .00000 

SIGMA-X  =       .0000,  SIGMA-Y  =       .0000 

TAU-XY  =  -.0003,  TAU-XZ  =       .0000,  TAU-YZ  =    -386.1477 

X  =     .0000,  Y  =    5.3750 

DEFLECTION  =  .00000 

SIGMA-X  =  .0000,  SIGMA-Y  =       .0000 

TAU-XY  =  .0352,  TAU-XZ  =    -133.8358,  TAU-YZ  =       .0000 

BOLT  STRESSES  AND  DEFLECTIONS 

POSITION  OF  BOLT  -  1,  X  =     .4500,  Y  =     .4500 

DEFLECTION  =   .00004 

AXIAL  STRESS  =    1343.1810 

BENDING  STRESSES  -  SIGMA-X  =    1992.4940,  SIGMA-Y  -   -2025.9750 

TAU-XZ  =    -239.2171,  TAU-YZ  =    -243.2367 

POSITION  OF  BOLT  -  2,  X  =    4.6500,  Y  =     .4500 

DEFLECTION  =   .00076 

AXIAL  STRESS  =   26910.0100 

BENDING  STRESSES  -  SIGMA-X  =    4666.1180,  SIGMA-Y  -  -40195.8900 

TAU-XZ  =    -560.2100,  TAU-YZ  =   -4825.8840 

POSITION  OF  BOLT  -  3,  X  =    8.8500,  Y  =     .4500 

DEFLECTION  =   .00076 

AXIAL  STRESS  =   26909.9900 

BENDING  STRESSES  -  SIGMA-X  =   -4666.1190,  SIGMA-Y  -   -40195.8800 

TAU-XZ  =     560.2101,  TAU-YZ  =   -4825.8810 

Figure  10. -Sample  output  for  12-bolt  cover  simulation. 


15 


POSITION  OF  BOLT  -  4,  X  =   13.0500,  Y  =     .4500 

DEFLECTION  =   .00004 

AXIAL  STRESS  =    1343.1860 

BENDING  STRESSES  -  SIGMA-X  =   -1992.4940,  SIGMA-Y  -   -2025.9810 

TAU-XZ  =     239.2170,  TAU-YZ  =    -243.2375 

POSITION  OF  BOLT  -  5,  X  =     .4500,  Y  =    3.7330 

DEFLECTION  =   .00060 

AXIAL  STRESS  =   21024.2900 

BENDING  STRESSES  -  SIGMA-X  =   30986.7400,  SIGMA-Y  -   -4455.9550 

TAU-XZ  =   -3720.2410,  TAU-YZ  =    -534.9780 

POSITION  OF  BOLT  -  6,  X  =   13.0500,  Y  =    3.7330 

DEFLECTION  =   .00060 

AXIAL  STRESS  =   21024.3300 

BENDING  STRESSES  -  SIGMA-X  =  -30986.7300,  SIGMA-Y  -   -4455.9610 

TAU-XZ  =    3720.2390,  TAU-YZ  =    -534.9787 

POSITION  OF  BOLT  -  7,  X  =     .4500,  Y  =    7.0160 

DEFLECTION  =   .00060 

AXIAL  STRESS  =   21032.0000 

BENDING  STRESSES  -  SIGMA-X  =   30998.0200,  SIGMA-Y  -    4454.3310 

TAU-XZ  =   -3721.5940,  TAU-YZ  =     534.7830 

POSITION  OF  BOLT  -  8,  X  =   13.0500,  Y  =    7.0160 

DEFLECTION  =   .00060 

AXIAL  STRESS  =   21032.0400 

BENDING  STRESSES  -  SIGMA-X  =  -30997.9900,  SIGMA-Y  -    4454.3370 

TAU-XZ  =    3721.5920,  TAU-YZ  =     534.7838 

POSITION  OF  BOLT  -  9,  X  =     .4500,  Y  =   10.3000 

DEFLECTION  =   .00004 

AXIAL  STRESS  =    1343.4310 

BENDING  STRESSES  -  SIGMA-X  =    1992.8410,  SIGMA-Y  -    2026.3340 

TAU-XZ  =    -239.2587,  TAU-YZ  =     243.2798 

POSITION  OF  BOLT  -  10,  X  =    4.6500,  Y  =   10.3000 

DEFLECTION  =   .00076 

AXIAL  STRESS  =   26910.2300 

BENDING  STRESSES  -  SIGMA-X  =    4665.8310,  SIGMA-Y  -   40196.4100 

TAU-XZ  =    -560.1755,  TAU-YZ  =    4825.9450 

POSITION  OF  BOLT  -  11,  X  -    8.8500,  Y  =   10.3000 

DEFLECTION  =   .00076 

AXIAL  STRESS  =   26910.2200 

BENDING  STRESSES  -  SIGMA-X  =   -4665.8330,  SIGMA-Y  -   40196.3900 

TAU-XZ  =     560.1758,  TAU-YZ  =    4825.9420 

POSITION  OF  BOLT  -  12,  X  =   13.0500,  Y  =   10.3000 

DEFLECTION  =   .00004 

AXIAL  STRESS  =    1343.4360 

BENDING  STRESSES  -  SIGMA-X  =   -1992.8400,  SIGMA-Y  -    2026.3410 

TAU-XZ  =     239.2586,  TAU-YZ  =     243.2807 

Figure  1 0.-Sample  output  for  1 2-bolt  cover  simulation-Continued. 
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Table  6. -Simulation  dimensions  and  physical  parameters 

Plate: 

Dimension  A  (X  direction)    in 

Dimension  B  (Y  direction)   in 

Thickness  (H)    in 

Elastic  modulus    kips/in2  .  . 

Yield  strength  kips/in2  .  . 

Poisson's  ratio 

Bolts: 

Bolt  1,  long  side: 

X  coordinate 

Y  coordinate 

Bolt  2,  corner: 

X  coordinate 

Y  coordinate 

Bolt  3,  short  side: 

X  coordinate 

Y  coordinate 

Elastic  modulus    kips/in2  .  . 

Area  moment  of  inertia 

about  any  diameter  (I) 

Shear  stress  factor  (KS) 

Cross  sectional  area   in2  .  . 

Thread  engagement 

length  (L)    in 

Diameter in 


13.5 

10.75 

0.5 

10,000 

33 

0.3 


8.85 
0.45 

0.45 
0.45 

0.45 

3.733 

30,000 

0.003068 

1.330 

0.1993 


0.85 
0.5 


COMPARISON  OF  SOLUTION  TECHNIQUES 


The  series  algorithm  was  previously  implemented  using 
VAX  FORTRAN  on  a  32-bit  word  length  VAX  computer. 
The  infinite  series  was  truncated  to  nine  terms.  This  im- 
plementation was  used  to  simulate  the  bolted  cover  of 
figure  9.  The  solutions  obtained  agree  well  with  the  so- 
lutions obtained  by  using  the  PC  implementation  with  the 
series  truncated  to  three  terms.  In  addition,  a  finite- 
element  model  of  the  bolted  cover  was  previously  solved 
using  ANSUL.  The  model  contained  312  elements  and 
208  nodes.    A  comparison  of  the  results  obtained  by  the 


three  solution  techniques  is  presented  in  table  7.  For  bolts 
1  and  3,  the  series  solution  gives  higher  values  for  axial 
stresses  and  lower  values  of  bending  stresses  than  the 
finite-element  solution.  Nevertheless,  the  total  axial  plus 
bending  stresses  for  the  two  methods  are  in  reasonable 
agreement.  For  the  corner  bolt  (bolt  2),  the  solutions  do 
not  agree  well  at  all.  The  difference  may  be  attributed  to 
the  constraint  equations  imposed  on  the  bolt  at  the  corner 
in  the  finite-element  analysis. 


Table  7.-Comparison  of  solution  techniques 

Deflection,  Axial 

Solution  technique  cm  stress, 

106  N/m2 

Bolt  1,  long  side: 

Finite  element   'o.OOSSSS  145.9 

MF  (9-term  series)   .002078  192.5 

PC  (3-term  series)    .001930  185.5 

Bolt  2,  corner: 

Finite  element   1.000640  27.1 

MF  (9-term  series)   .000104  9.7 

PC  (3-term  series)    .000102  9.3 

Bolt  3,  short  side: 

Finite  element    1.002103  114.3 

MF  (9-term  series)   .001631  151.0 

PC  (3-term  series)    .001524 145.0 

MF       Main  frame  computer. 

PC        Personal  computer. 

'includes  rigid  body  motion  as  well  as  elongation  of  bolts  because  of  where  constrain  was  applied  in  the  model. 


Bending  stress  at 
bolt  top,  106N/m2 


30.6  413.9 
34.8  300.0 
32.2  277.2 

49.7  55.9 
15.0  15.3 
13.7  14.0 

321.4  31.3 

232.2  34.3 

213.7  30.7 
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RESULTS  OF  CALCULATIONS  PERFORMED  ON  REPRESENTATIVE  PLATES 


The  PC  implementation  of  the  series  solution  algorithm 
can  be  used  to  estimate  the  deflection  at  the  center  of  a 
simply  supported  flat  plate  subjected  to  a  uniform  pres- 
sure.   The  plate  is  treated  as  a  cover  without  any  bolts. 


Deflections  were  determined  for  10  different  plates  sub- 
jected to  various  uniform  static  pressures.  The  results  of 
the  simulations  are  shown  in  table  8. 


Table  8.-Results  of  simulations  performed  on  various  plates 


Plate 


Dimensions,  in 


B 


H 


Elastic 
modulus,  kips/in2 


Pressure, 
psi 


Deflection, 
in 


1   

2 

3 

4 

5 

6 

7 

8 

9 

10 

A  Length. 

B  Width. 

H  Thickness. 


48.75  30 


48.75  30 


48.75  30 


48.75  30 


0.75 


1.5 


48.75  30  2.5 

30  30.25  1 

20.5  28.69  1 

10  24  1 

10.75  13.5  1 

24  36      .5 


10,000 

10,000 
29,000 

10,000 

10,000 
29,000 
10,700 
10,700 

30,000 

29,000 


50 
100 
150 

50 
150 
250 

104.2 
106 

100 
200 
300 

300 

174.8 

90 

81.5 

150 
180 

125 


0.86425 
1.72850 
2.59275 

.36461 
1.09382 
1.82302 

.26852 
.26654 

.21606 
.43212 
.64819 

.14001 

.21525 

.11223 

.00911 

.27684 
.33220 

.94241 


CONCLUSIONS 


The  U.S.  Bureau  of  Mines  computer  program  described 
in  this  report  can  be  used  as  a  convenient  and  fast  means 
of  evaluating  the  performance  of  XP  enclosures.  The 
program  can  be  used  to  evaluate  components  of  the  en- 
closure as  well  as  the  complete  enclosure.  Useful  in- 
formation concerning  common  materials  of  construction 
and  design  constraints  are  included  in  the  accompanying 
data  base. 


XP  enclosure  designers  and  manufacturers  can  use  the 
program  to  determine  if  the  design  is  adequate  before 
submitting  the  enclosure  for  costly  and  time-consuming 
certification  testing.  MSHA  can  use  the  program  to  screen 
enclosures  before  putting  them  through  the  lengthy  and 
costly  series  of  inspections  and  explosion  tests.  Enclosures 
that  are  clearly  inadequate  can  be  returned  to  the 
manufacturer  for  redesign. 
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APPENDIX  A.-FILES  ON  DISK  1,  XP  ENCLOSURE  SCREENING  PACKAGE 


LIST 

DIR 

0 

3-15-89 

2: 13a 

XP 

EXE 

199236 

7-10-87 

9:40a 

DISP1 

DOC 

709 

6-04-87 

8:57a 

DISP2 

DOC 

880 

8-27-86 

8:21a 

DISP3 

DOC 

940 

9-16-86 

4:34a 

DISP4 

DOC 

786 

6-16-87 

ll:55a 

DI5P5 

DOC 

981 

7-25-86 

l:56p 

DISP6 

DOC 

486 

7-25-86 

2J02p 

DISP7 

DOC 

1166 

7-28-86 

10:42a 

DISP8 

DOC 

1166 

7-28-86 

11:00a 

DISP9 

DOC 

813 

8-04-86 

2:35a 

DISP10 

DOC 

1069 

8-01-86 

ll:06a 

DISP11 

DOC 

737 

8-04-86 

2:48a 

DISP12 

DOC 

1353 

8-19-86 

l:52a 

DISP14 

DOC 

1029 

8-05-86 

10:54a 

DISP15 

DOC 

1048 

8-19-86 

2: 12a 

DISP16 

DOC 

798 

8-19-86 

3:31p 

DISP17 

DOC 

454 

8-20-86 

9:26a 

DISP18 

DOC 

474 

9-16-86 

12:56p 

DISP19 

DOC 

478 

9-16-86 

12:57p 

F00T1 

DOC 

1136 

1-14-86 

4:22a 

WELDEF 

DOC 

491 

1-14-86 

7:53a 

FLANEPR 

DAT 

269 

10-20-85 

2:41a 

WELDJE 

DAT 

120 

10-20-85 

2:49a 

MATSTEEL  DAT 

346 

1-15-86 

10:23a 

STEELDES  DAT 

563 

8-05-86 

9:23a 

MATSLASS  DAT 

461 

6-10-86 

10: 11a 

GLASSDES  DAT 

352 

6-10-86 

10:23a 

MATSEAL 

DAT 

463 

6-10-86 

10:25a 

INPUT 

DAT 

216 

10-19-85 

5:45a 

FLAME1 

DAT 

810 

1-01-80 

12:58a 

MATAL 

DAT 

349 

6-10-86 

9:32a 

ALDES 

DAT 

583 

6-10-86 

12:39p 

SEALDES 

DAT 

176 

6-10-86 

10:30a 

34  File 

(5) 

32944  bytes  free 

■BBBBBBB^B 


■■■ 


19 


APPENDIX  B.-FILES  ON  DISK  2,  XP  ENCLOSURE  SCREENING  PACKAGE 

SOURCE  CODE 


Volume 

in  drive  B  has 

no  label 

Directory  of 

B:\ 

TAB52 

FOR 

1217 

8-19-86 

5:25p 

DISPLAY 

FOR 

1895 

6-16-B7 

l:57p 

XP 

FOR 

16708 

6-16-87 

12:28p 

WELD 

FOR 

1978 

6-16-87 

11:13a 

S 

FOR 

2290 

6-17-87 

2:57p 

CC 

OBJ 

46290 

7-10-87 

9:33a 

BB 

OBJ 

31034 

7-09-87 

3:32p 

W 

FOR 

8562 

6-17-87 

3:09p 

RETURN 

FOR 

331 

6-04-86 

3:52a 

WELDE 

FOR 

362 

6-09-86 

3:56a 

MAT 

FOR 

3360 

8-26-86 

6:58a 

VOLUME 

FOR 

3887 

6-05-86 

l:l6p 

PRINT 

FOR 

464 

6-05-86 

2:09a 

FLAMEP 

FOR 

1630 

8-06-86 

4:44p 

FOOT 

FOR 

1876 

8-26-86 

9:32a 

PENET 

FOR 

10578 

8-06-86 

i:57p 

R 

FOR 

7948 

8-26-86 

7:00a 

DISPLAY 

OBJ 

3055 

6-16-87 

ll:05a 

WELD 

OBJ 

3288 

6-16-87 

11:14a 

FIS4 

OBJ 

1943 

7-09-87 

3:20p 

TAB52 

OBJ 

2356 

8-19-86 

4:5ip 

RETURN 

OBJ 

763 

6-04-86 

4:04a 

WELDE 

OBJ 

1040 

6-09-86 

3:57a 

MAT 

OBJ 

6395 

8-27-86 

9:08a 

VOLUME 

OBJ 

6950 

6-05-86 

l:17a 

PRINT 

OBJ 

1077 

6-05-86 

2:lla 

FLAMEP 

OBJ 

2816 

8-06-86 

4:45p 

FOOT 

OBJ 

3045 

8-26-86 

9:34a 

PENET 

OBJ 

17379 

8-06-86 

2:00p 

SIM 

OBJ 

11328 

6-12-86 

10:57a 

R 

OBJ 

13990 

8-27-86 

9: 16a 

W 

OBJ 

14973 

6-17-87 

3:39p 

TEST1 

OBJ 

3147 

7-09-87 

l2:5lp 

S 

OBJ 

3137 

6-17-87 

3:28p 

XP 

OBJ 

24270 

6-16-87 

11: 17a 

PRT 

1778 

7-24-87 

2:39a 

CC 

FOR 

25693 

7-10-87 

8:08a 

FIS4 

FOR 

651 

7-09-87 

3:l9p 

TEST1 

FOR 

1047 

7-09-87 

12:50p 

APPB 

DOC 

1679 

1-01-80 

3:39a 

BB 

FOR 

18968 

7-09-87 

3:35p 

XPMAIN 

FOR 

16512 

7-24-87 

11:13a 

LIST 

DIR 

0 

3-15-89 

3:02a 

43  Filets)  10240  bytes  -free 


20 


APPENDIX  C.-FORTRAN  SOURCE  CODE  LISTING  OF  XP  SCREENING  PACKAGE 

C     ************************************************************ 
C     *  * 

C    *  XPMAIN  * 

c    *  * 

c   I*********************************************************** 

CHARACTERS  FNAME.ID 

CHARACTERS  pos2(10),CLR(4) ,CR,ESC,ST,P0S1 (9), REV (4), NORM (7) 

CHARACTER*!  TOP (9) 

DIMENSION  AAA (12, 4), PAN (20, 3) 

DATA  CLR/'  ','[', »2V J'/ 


DATA  P0S2/' 
DATA  P0S1/' 
DATA  REV/'  ' 
DATA  NORM/' 


,'[','1V0V?V1VHV  •»'[', 'K'/ 

,'[V7ViVlVHV  VCVKV 

'[','7','m'/ 

,'[','0'.'ni','  ','C, 'A'/ 
DATA  TOP/'   VCV1V5V1VHV   ','C','K'/ 
ESC=CHAR(27) 
CLR(1)=ESC 
T0P(1)=ESC 
T0P(7)=ESC 
P0S2(1)=ESC 
P0S2(8)=ESC 
P0S1(1)=ESC 
P0S1(7)=ESC 
CR=CHAR(13) 
REV(1)=ESC 
N0RM(1)=ESC 
N0RM(5)=ESC 
ST=CR 
NNN=1 

WRITE  (*,*)  CLR 
WRITE(*,1330)  P0S2 

1330  FORMAT  ('   M0A1) 
IFLAGV=0 
VOLUMES. 
IRUG=0. 

IT=0. 

WRITE  (*,*)  CLR 
WRITE  (•,1331)  P0S1 
WRITE  (*,♦)  REV 

1331  FORMAT  ('  ',9AD 

WRITE!*.*)'  XF  ENCLOSURE  SCREENING  PROGRAM   V5.0   08/11/88' 

WRITE!*,*)'  WRITTEN  BY     F.  T.  Duda      ' 

WRITE  (♦,*)  NORM 

Write(*,») 

WRITE(*,*) 

WRITE(*,*) 

WRITE(*,*) 
10   CONTINUE 

DO  3333  1=1,1000 

DO  3333  J= 1,200 
3333  CONTINUE 

WRITE  (♦,*)  CLR 

WRITE  (*,1330)  P0S1 
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NARATIVE  TO  BE  PRINTED  ON  SCREEN 


The  Bureau  of  Mines  expressly  declares  that  there  are 
no  warranties  express  or  implied  which  apply  to  the 
software  contained  herein.  By  acceptance  and  use  of 
said  software!  which  is  conveyed  to  the  user  without 
consideration  by  the  Bureau  of  Mines,  the  user  hereof 
expressly  waives  any  and  all  claims  for  damage  and/or 
suits  for  or  by  reason  of  personal  injury  or  property 
damage?  including  special?  consequential  or  other 
similar  damages  arising  out  of  or  in  any  way 
connected  with  the  use  of  the  software  contained 
herein.' 


WRITEC*,*) 

WRITER,*) 

WRITE**,*) 

WRITE**, *> 

WRITER, *) 

WRITE(*,*> 

WRITE**,*) 

WRITE**,*) 

WRITE**,*) 

WRITE (*!*) 

WRITE**,*) 

WRITE(*,*) 

WRITE**,*) 

WRITE**,*) 

CALL  RETURN *CLR) 

WRITE**, *)CLR,P0S2 

IGO=20 

NUMB=20 

CALL  DISPLAY !IG0, NUMB, CLR) 

IF  (IG0.EQ.1)  THEN 

6100  OPEN *23,FILE='B: RESULT, DOC ,STATUS=' NEW' ) 
ELSE 

6101  0PEN(23,FILE='C:RE5ULT.B0C ,STATUS='NEW' ) 
ENDIF 

6103  CONTINUE 

IG0=21 

NUMB=20 

CALL  DISPLAY (IGO, NUMB, CLR) 
1001  WRITE  **,*)  CLR,P0S2 

WRITE  **,001) 
001  FORMAT  (IX,'  ENTER  THE  IDENTIFICATION  OF  THE  ENCLOSURE  TO  BE 
1SCREENED' ) 

WRITE**,*) 

WRITE**,*) 

WRITE**,*)' 

WRITE**,*) 

WRITE**, +)'  IF 

WRITE**,*)' 

WRITE**,*) 

READ  (*,310)  ID 
310  FORMAT  <A) 

WRITE**,*)  CLR,P0S2 

WRITE**,*)'  THE  IDENTIFICATION  FOR  YOUR  ENCLOSURE  IS:' 

WRITE**,*) 

WRITE!*, 355)  ID 
355  FORMAT (15X,64A) 

WRITE**,*) 

WRITE**,' (A\)')' 

READ**, 310)  ST 

IF  (ST.  EG. 'Y\  OR.  ST.  EG 

GO  TO  25 

ELSE 

GOTO  1001 

ENDIF 


THE  NAME  CAN  BE  UP  TO  64  CHARACTERS  IN  LENGTH' 

HE  ENCLOSURE  DOES  NOT  HAVE  AN  IDENTIFICATION,' 
HIT  RETRN' 


IS  THIS  THE  CORRECT  IDENTIFICATION?  <Y)' 


V 


, OR. ST. EG.'  ')  THEN 
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25   WRITE(23,392)'  THE  ENCLOSURE  BEING  SCREENED  IS  MD 

WRITE(23.*> 
2  WRITE  (*,*)  CLR.P0S2 

WRITEt*,*)'  DO  YOU  WANT  TO:' 
WRITE  (*,*) 

WRITEt*,*)'  (1)      DO  YOU  WANT  TO  PERFORM  A  SINGLE  TASK' 
WRITEt*,*) 

WRITEt*,*)'  OR' 

WRITE(*,*) 

WRITEt*,*)'  (2)      DO  YOU  WANT  TO  BE  DIRECTED  THROUGH  THE' 
WRITE(*,*)'        SCREENING  PROCESS?' 
WRITEt*,*)'         ' 
WRITEt*,*)'  OR' 

WRITE(*,*) 

WRITE(*,*)'  (3)      DO  YOU  WANT  TO  EXIT' 
WRITE  (*,*) 
WRITE  (*,*) 

WRITE  (*,'(A\)')  '  INPUT  YOUR  CHOICE-)' 
READ  (*,'(BN,I6)')I1 
IF  (Il.EQ.l)  THEN 
IFLAGX=1 

ELSEIF  (I1.EQ.2)  THEN 
IFLA6X=2 

ELSEIF  (I1.EQ.3)  THEN 
GOTO  20 
ELSE 

WRITE  (*,*)'  YOU  ONLY  HAVE  TWO  CHOICES,  (1)  OR  (2)' 
WRITE  (*,*)'  TRY  AGAIN' 

GO  TO  2 
ENDIF 
IF  (IFLAGX.EQ.2)  THEN 
WRITE(*,*)  CLR.P0S2 

WRITE(*,*)'  YOU  HAVE  CHOSEN  TO  BE  DIRECTED  THROUGH  THE' 
WRITE(*,*)'  SCREENING  PROCESS' 
WRITE(*,*) 
WRITE(*,») 
CALL  RETURN (CLR) 
WRITE(*,*)  CLR.P0S2 

WRITE(*,*)'  YOU  WILL  BE  DIRECTED  THROUGH  THE  FOLLOWING' 
WRITE(*,*)'  STEPS:' 


SCREEN  FOR  FLAME  PATH  REQUIREMENTS' 

SCREEN  FOR  WELD  QUALITY  EFFICIENCY' 

SCREEN  FOR  ACCEPTABLE  CONSTRUCTION' 

MATERIALS' 

SCREEN  FOR  ENCLOSURE  STRENGTH' 

SCREEN  FOR  ENCLOSURE  RUGGEDNESS' 


WRITE(*,*) 

WRITE(*,*)' 

(1) 

WRITEt*.*)' 

(2) 

WRITEt*,*)' 

(3) 

WRITEt*,*)' 

WRITEt*,*)' 

(4) 

WRITEt*,*)' 

(5) 

WRITEt*,*) 

WRITEt*,*) 

CALL  RETURN (CLR) 

12=1 

GO  TO  3 

ENDIF 

WRITE  (*,*)  CLR,P0S2 

WRITE  (*,*)  CLR,P0S2 
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312  WRITE*.*,  *)'  DO  YOU  WANT  TO:' 
WRITE  (*,*) 
WRITE  (*,*) 

WRITE (*,*)'  (1)  SCREEN  FOR  FLAME  PATH  REQUIREMENTS' 
WRITE(*.*)»  (2)  SCREEN  FOR  WELD  QUALITY  EFFICIENCY' 
WRITE**,*)'  (3)  SCREEN  FOR  ACCEPTABLE  CONSTRUCTION' 
WRITE**,*)'       MATERIALS' 

WRITE(*.*)'  (4)    SCREEN  FOR  ENCLOSURE  AND/OR  WINDOW  STRENGTH' 
WRITE (*,*)'  (5)    SCREEN  FOR  ENCLOSURE  RUGGEDNESS' 
WRITE**,*)'  (6)    END  PROGRAM' 
WRITE  (*,*) 
WRITE!*,*) 

WRITE  (*,'(A\)')  '  INPUT  YOUR  CHOICE — )' 
READ  (*.'(BN.I6)')I2 

3  CONTINUE 

IF  (12.EQ.1)  THEN 
♦MODULE  FOR  FLAME  PATH  REQUIREMENTS********************************** 

WRITE  (*,*)  CLR.P0S2 

WRITE(*,*)'  THIS  SECTION  SCREENS  FOR  FLAME  PATH  REQUIREMENTS' 

WRITE**, *) 

WRITE(*.*) 

CALL  RETURN (CLR) 
♦INSTRUCTIONS  GO  HERE 
390   CONTINUE 

WRITE**,*)  CLR.P0S2 

WRITE  (*,*)  '  DO  YOU  WANT  TO' 

WRITE  (*,♦) 

WRITE  !*.*) 

WRITE (*,*)' (1)  VIEW  FLAME  PATH  REQUIREMENTS  FOR  YOUR  ENCLOSURE' 

WRITE  (*,*) 

WRITE'*,*) '(2)  VIEW  FOOTNOTES  OF' 

WRITE(*,*)'    30  CFR  PART  18  REQUIREMENTS' 

WRITE*.*,*) 

WRITER, *)' (3)  NEITHER' 

WRITE  (*,*) 

WRITE  (*,*) 

WRITE  (*,'(A\)')'  INPUT  YOUR  CHOICE-)' 

READ  (*,'*BN,I6)')  Jl 
IF  (Jl.EQ.l)  THEN 

4  WRITE  (♦,♦)  CLR,P0S2 

311    WRITE(*,*>'  DO  YOU  WANT  TO:' 
WRITE  (♦.♦) 

WRITER,*)'  (1)      ENTER  INTERNAL  VOLUME  MEASURED  FROM' 
WRITE(*,*)'        ENCLOSURE  IN  CUBIC  INCHES' 
HRiTEt*.*)'  OR' 

WRITE**,*)'  (2)      TYPE  IN  ENCLOSURE  DIMENSIONS  KEYED  TO' 
WRITE**, *)'         COMPONENTS.' 
WRITE  (*,*> 
WRITE  (*,*) 

WRITE  (*,'(A\)')  '  INPUT  YOUR  CHOICE--)' 
READ  **,'(BN,I6)')I1 

IF  (Il.EQ.l)  THEN 

WRITE  (*,♦)  CLR.P0S2 

WRITE  (*,*)'  ENTER  VOLUME  IN  CUBIC  INCHES' 

READ  (*,*)  VOLUME 
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WRITE  (*,*)  CLR,  P0S2 

WRITER,  *)'  THE  VOLUME  OF  ENCLOSURE  '.ID,'  IS  '.VOLUME,' 
1  CU.  IN.' 

ELSEIF  (Il.EQ.2)  THEN 
CALL  VOLUM1 (CLR, P0S2, VOLUME) 
ELSE 

WRITER, *)'  YOU  ONLY  HAVE  TWO  CHOICES' 
60  TO  4 
ENDIF 

WRITE  (♦,*) 
WRITEf*,*) 
CALL  RETURN (CLR) 
WRITE  (*,*)CLR,P0S2 
WRITE  (*,*)'  DO  YOU  WANT  TO:' 
WRITE  (*,*) 

WRITE(*,*)'    (1)  CONTINUE  WITH  THE  PROGRAMS  CALCULATED' 
WRITE(*,*)'        VOLUME  OF' 

WRITE(*,*)'  ', VOLUME,'  CUBIC  INCHES' 

WRITE(*,*) 

WRITE (*,*)'    (2)  ENTER  THE  INTERNAL  VOLUME  ' 
WRITE(*,*) 
WRITE(*,*) 

WRITE(*,'(A\)')'  INPUT  YOUR  CHOICE — )' 
READ  (*,'(BN,I6)')IXX 
IF  (IXX.EQ.2)  THEN 

WRITE(*,*)'  ENTER  VOLUME  IN  CUBIC  INCHES' 
READ  (*,*)  VOLUME 
ELSEIF  (IXX.EQ.l)  THEN 
GO  TO  5 
ELSE 
WRITE(*,*)'  YOU  ONLY  HAVE  TWO  CHOICES,  (1)  OR  (2)' 

WRITE  (*,*)'         TRY  AGAIN ' 

GO  TO  6 
ENDIF 

IFLAGV=1 
V=VOLUME 

WRITE(*7*>  CLR.P0S2 
IF  (V.LT.45.)  THEN 

WRITE  (*,*)'  VOLUME  IS  LESS  THAN  45  CU.  IN.  ' 
IFLAG=1 

ELSEIF  (V.GE.45.  .AND.V.LE.124.)  THEN 
WRITE  (*,*)'  VOLUME  IS  BETWEEN  45  AND  124  CU.  IN. 
1  INCLUSIVE' 

IFLAG=2 
ELSE 

WRITE  (*,*)  '  VOLUME  IS  MORE  THAN  124  CU.  IN.' 
IFLAG=3 
ENDIF 
WRITER,  *) 
WRITE(*,*) 
CALL  RETURN (CLR) 
ELSEIF  (J1.EQ.2.0R.J1.EQ.3)  THEN 
GO  TO  11 
ELSE 
WRITE  (*,*)'  YOU  ONLY  HAVE  THREE  CHOICES' 
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GO  TO  3 
ENDIF 
11       CONTINUE 

IF  (Ji.EQ.l)  THEN 
WRITE!*,*)  TOP 
WRITE(*,*)'    FLAME  FATH  REQUIREMENTS  SPECIFIED  IN  PARAGRAPH 
1  18.3  OF  30  CFR  WILL  BE' 
WRITER*) 

WRITE(*,*) 'DISPLAYED  ON  THE  NEXT  SCREEN.  SINCE  THESE 
1  REQUIREMENTS 
l  ARE  VOLUME  DEPENDENT, ' 
WRITE(*,») 

WRITE (*,*)  'THE  REQUIREMENTS  SHOWN  ARE  FOR  YOUR  ENCLOSURE,  MD 
WRITER, *) 'WHICH  HAS  A  VOLUME  OF  ',V,'   CUBIC  INCHES. ' 
WRITE(*,*) 

WRITE(*,*)'    IT  IS  EXPECTED  THAT  YOU  WILL  COMPARE  THE 
1ENCLOSURE  VALUES  EITHER  TAKEN  FROM' 
WRITER, *) 

WRITE(*,*)'THE  ENCLOSURE  DRAWINGS  OR  ACTUAL  MEAURSEMENTS  FROM 
1THE  ENCLOSURE.' 
WRITE(*,*5 

WRITE<*,*)'    AFTER  YOU  HAVE  HAD  AN  OPPORTUNITY  TO  VIEW  THE 
1REQUIREMENTS  AND  COMPARE  YOUR' 
WRITER, *) 

WRITER,*) 'ENCLOSURE  MEASUREMENTS  TO  THE  REQUIREMENTS,  YOU  WILL 
1BE  ASKED  IF  THE  ENCLOSURE' 
WRITE(*,*) 

WRITE(*,*) 'SATISFIES  THE  REQUIREMENTS  OR  NOT.' 
WRITEC+,*) 
WRITE(*,») 
WRITE(*,*) 
CALL  RETURN (CLR) 

CALL  FLAMEP(CLR,IFLAG) 
GO  TO  391 

ELSEIF  (J1.EQ.2)  THEN 
IGO=l 
NUMB=21 

CALL  DISPLAY (160. NUMB, CLR) 
GO  TO  391 

ELSEIF  (J1.EQ.3)  THEN 
IFUFLAGX.EQ.2.AND.  IFLAGV.EQ.0)  THEN 
WRITER, *)  CLR,P0S2 
WRITE!*,*)'  YOU  MUST  CHOOSE  1' 
WRITE  (*,*) 
WRITE(*,*) 
CALL  RETURN (CLR) 
60T0  390 

ELSEIF ( IFLA6X . EQ. 2. AND. IFLAGV. EQ. 1 )  THEN 
WRITE'.*, ♦)  CLR,P0S2 

WRITE  (*,*)  '  SINCE  YOU  ARE  BEING  GUIDED  THROUGH' 
WRITE(*,*)'  THE  SCREENING  PROCESS,  YOU  WILL  ' 
WRITE(*,*)'  BE  DIRECTED  TO  THE  NEXT  STEP  WHICH  IS  TO' 
WRITE(*,*)'  SCREEN  FOR  WELD  QUALITY  EFFICIENCY' 
12=2 
WRITE(*,») 
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WRITE(*,«) 

CALL  RETURN  (OR) 

60  TO  3 

ELSEIF  (IFLA6X.EQ.1)  THEN 

WRITE(*,*)  CLR.P0S2 

WRITE(*,*)'Y0U  ARE  DONE  WITH  THIS  SECTION.  YOU  ONLY' 

WRITE(*,*) 'WANTED  TO  SCREEN  FOR  FLAME  PATH  REQUIREMENTS' 

WRITE(*,*) 

WRITE(*,*) 

WRITER, *) 

CALL  RETURN (CLR) 

GOTO  2 

ENDIF 
ELSE 

WRITE  (*,*)'  YOU  ONLY  HAVE  THREE  CHOICES,  1,  2,3' 
GO  TO  390 
ENDIF 

391  CONTINUE 
IG0=2 
NUMB=21 
WRITE(*,*)  CLR.P0S2 

WRITE(*,*)'  DOES  THE  ENCLOSURE  MEET  THE  FLAME  PATH' 
WRITE(*,*)'  REQUIREMENTS  AS  SPECIFIED  IN  PARAGRAPH' 
WRITE(*,'(A\)')'  18.31  OF  30  CFR  ?  <Y>' 
READ(*,310)  ST 

392  FORMAT (1X,32A,64A) 

IF  (ST.EQ.'Y'.OR.ST.EQ.'  '.OR.ST.EQ.'Y')  THEN 
WRITE(23,*>'    THE  ENCLOSURE  DOES  MEET  THE  FLAME  PATH' 
WRITE(23, *)' REQUIREMENTS  AS  SPECIFIED  IN  PARAGRAPH  18.31' 
WRITE(23,*)'0F  30  CFR.' 
IFLAGN=1 
WRITE(23,*) 
ELSE 

WRITE(23,*)'    THE  ENCLOSURE  DOES  NOT  MEET  THE  FLAME  PATH' 
WRITE(23,*) 'REQUIREMENTS  AS  SPECIFIED  IN  PARAGRAPH  18.31' 
WRITE(23,*)'0F  30  CFR.' 
IFLAGN=0 
WRITE(23,*> 
ENDIF 
GO  TO  390 
21    CONTINUE 

ELSEIF  (I2.EQ.2)  THEN 

CALL  WELDCH(ID,CLR,P0S2,AAA) 

WRITE (*i*(A\)>) 

IF  (IFLAGX.EQ.l)  THEN 

WRITER*,*)  CLR,P0S2 

WRITE(*,*)'  YOU  ARE  DONE  WITH  THIS  SECTION.  YOU  ONLY' 

WRITE (*,*)'  WANTED  TO  SCREEN  FOR  WELD  QUALITY.' 

WRITE(*,*> 

WRITE(*,*) 

CALL  RETURN (CLR) 

GO  TO  2 

ELSE 

12=3 

WRITE(*,*)  CLR,P0S2 
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WRITE**,*)'    SINCE  YOU  ARE  BEING  GUIDED  THROUGH  THE' 
WRITE**,*) 'SCREENING  PROCESS,  YOU  WILL  BE  DIRECTED  TO  THE' 
WRITE(*,*)'NEXT  STEP  WHICH  IS  TO  SCREEN  FOR  ACCEPTABLE' 
WRITE**,*) 'MATERIALS  OF  CONSTRUCTION.' 
WRITE**,*) 
WRITE**,*) 
CALL  RETURN (CLR) 
GOTO  3 
ENDIF 
ELSEIF  (I2.EQ.3)  THEN 
WRITE**,*)  CLR.P0S2 

WRITE**,*)'  THIS  SECTION  SCREENS  FOR  ACCEPTABLE' 
WRITE**,*)'  MATERIALS  OF  CONSTRUCTION' 
WRITE**,*) 
WRITE**,*) 
CALL  RETURN (CLR) 
CALL  MAT(P0S2,CLR,T0P) 

IF  (IFLAGX.EQ.l)  THEN 

WRITE**, *)CLR,P0S2 

WRITE**,*)'    YOU  ARE  DONE  WITH  THIS  SECTION.  YOU  ONLY' 

WRITE**,*) 'WANTED  TO  SCREEN  FOR  ACCEPTABLE  MATERIALS  OF' 

WRITE**,*) 'CONSTRUCTION.' 

WRITE**,*) 

WRITE**,*) 

CALL  RETURN *CLR) 

60  TO  2 

ELSE 

12=4 

WRITE**,*)  CLR,P0S2 

WRITE**,*)'    SINCE  YOU  ARE  BEING  GUIDED  THROUGH  THE' 

WRITE**,*) 'SCREENING  PROCESS,  YOU  WILL  BE  DIRECTED  TO  THE  ' 

WRITE**,*) 'NEXT  STEP  WHICH  IS  TO  SCREEN  THE  ENCLOSURE  FOR' 

WRITE**,*) 'STRENGTH  OF  CONSTRUCTION.' 

WRITE**,*) 

WRITE**,*) 

CALL  RETURN *CLR) 

GOTO  3 

ENDIF 
ELSEIF* 12. EQ. 4)  THEN 
WRITE**,*)  CLR,P0S2 
WRITE**,*)'  THIS  SECTION  SCREENS  FOR' 
WRITE**,*)'  ENCLOSURE  STRENGTH' 
WRITE**,*) 
WRITE**,*) 
CALL  RETURN (CLR) 
WRITE**, *)CLR,P0S2 

WRITE(*,'(A\)')'  DO  YOU  WANT  TO  VIEW  INFORMATION  ?(Y)' 
READ  (*,310)  ST 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 
60  TO  506 
ELSE 

GO  TO  505 
ENDIF 
506   160=4 
NUMB=21 
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CALL  DISPLA/(IGO, NUMB, CLR) 
IG0=5 
NUMB=8 

CALL  BISPLAYdGO, NUMB, CLR) 
505        CONTINUE 

WRITE(*,*)CLR,P0S2 

CALL  STREN ( ID.CLR, P0S2, IFLAGX, TOP, VOLUME, IRU6, PAN, IT) 

CALL  RETURN >CLR) 

GOTO  800 

elseif  (i2.eq.5)  then 

write(*,*)clr,p0s2 

writet*,*)'  this  section  screens  for  ruggedness' 

write(*,*)'  irug:  mrug,'  it:  ',it 

DO  4056  IPP=1,IT 
WRITE(*,*) 

WRITE(*,*>  'RM     THICKNESS  IN  FT.     A  TINES  B' 
4056   WRITE(*,*)PAN(IPP,1),PAN(IPP,2),PAN(IPP,3) 
WRITE(*,*> 
CALL  RETURN (CLR) 
IF(IRUG,EQ.0)  THEN 
WRITE(*,*)CLR,P0S2 

WRITE**, *>'  YOU  MUST  DETERMINE  THE  RESISTANCE  OF  THE  COVER' 
WRITE (*,*)'  AND  PANELS  BEFORE  YOU  CAN  DO  THIS  SECTION.' 
WRITE(*,*) 

WRITE(*,*)'  YOU  SHOULD  CHOOSE  4  FROM  THE  MAIN  MENUE-' 
WRITE(*, *) '  SCREEN  FOR  ENCLOSURE  STRENGTH. . . . ' 

WRITE(*,*)'  AND  THEN  CHOOSE  3  FROM  NEXT  MENUE ' 

WRITE(*,*)'  CHECK  BODY  PANELS  FOR  STRENGTH.5 
WRITE(*,*) 
WRITE(*,*> 
CALL  RETURN (CLR) 
ELSE 
WRITE(*,*)CLR,P0S2 

WRITEt*.*)'  YOU  HAVE  CHECKED  SIT,'  PANELS  FOR  STRENGTH.' 
WRITER*) 
WRITE(*,*)'  YOU  CAN  EITHER  CHECK  ONLY  THOSE  PANELS  FOR  ' 
WRITE(*,*!'  RUGGEDNESS  OR  YOU  CAN  60  TO  THE  SECTION' 
WRITEt*.*)'  THAT  WILL  DETERMINE  THE  PANEL  RESISTANCE' 
WRITE(*,*)'  FOR  ALL  PANELS.' 
WRITE!*,*) 
WRITEt*-*) 

WRITE(*,'(A\)')'  DO  YOU  WANT  TO  CHECK  THOSE  F0R  RUGGE0NES5?<\ 
READ(*,310)STT 
IF  (STT.EQ.'Y'.OR.STT.EQ.'y'.OR.STT.EQ.'  MTHEN 
WRITE(*>*)CLR,P0S2 
DO  4550  IPP=l,IT 
RKE=PAN(IPP,i)*PAN(IPP,2) 

WRITE (23,*) 

WRITE(23,*)'  RKE  FOR  PANEL  MPP,'  IS  ',RKE 

wnte(23,*) 

write(23,*)'  RESULTS  OF  RUGGEDNESS  SCREENING.' 

WRITE!23,») 

IF  (RKE. GT. 3700.)  THEN 

WRITE(23,*) 

WRITE (23,*)'  PANEL  ',IT,'  DOES  NOT  MEET  THE  CRITERIA' 
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WRITE<23,*)'  FOR  SURFACES  EXPOSED  TO  ROOF  FALLS' 
WRITE<23,*> 
ENDIF 

IF(RKE.GT.8000.)  THEN 
WR1TE<23,*) 

WRITEC23,*}'  PANEL  ML'  DOES  NOT  MEET  THE  CRITERIA' 
WRITE (23,*)'  FOR  SURFACES  EXPOSED  TO  SIDE  IMPACT.' 
WRITE<23,*) 
ENDIF 

IF(RKE.GT.800.)  THEN 
WRITE (23,*) 

WRITE (23,*)'  PANEL  ML'   DOES  NOT  MEET  THE  CRITERIA' 
WRITE(23,*)'  FOR  PROTECTED  SURFACES  TO  IMPACT.' 
WRITE (23,*) 
ENDIF 
4550   CONTINUE 
ELSE 

GO  TO  S00 
ENDIF 
ENDIF 
ENDIF 
20    CONTINUE 

WRITE(*,*)'  YOU  ARE  DONE.' 
END 

C  Mm***************************************************** 

C  *  * 

C  *  SUBROUTINE  BODY  (BB.FOR)  * 

C  *  * 

C       M**  *********  *******M*»  ****************************  ****** 

SUBROUTINE  80DY(CLR,POS2, I FLA6P, VOLUME, IRUG.PAN, IT, TOP) 
CHARACTER*!  CLR(4) ,POS2(10),TOP<9) 
DIMENSION  AT (5,3) 
DIMENSION  PAN (20, 3) 
WRITE(*,*)  CLR,POS2 


THIS  SECTION  CHECKS  BODY  PANELS  FOR  AN  INTERNAL' 
STATIC  PRESSURE  OF  150  PSIG  AND  FOR  A  PRESSURE  ' 
LOADING  PRODUCED  BY  AN  INFERNAL  EXPLOSION' 
BECAUSE  NO  DEFORMATION  CRITERIA  IS  SPECIFIED  IN  ' 
30  CFR  PART  18  FOR  A  STATIC  PRESSURE  OF  150PSIG,' 
THE  CALCULATION  IS  BASED  ON  A  LIMIT  OR  YIELD  ' 
LINE' 


WRITE(*,*) 
WRITE(*,*) 
WRITE(*,*> 
WRITE!*,*) 
WRITE(*,*) 
WRITER,*) 
WRITE(*,*) 
WRITE(*,*) 
WRITE(*,*) 
CALL  RETURN (CLR) 
IS=0 
IT=0 

WRITE(*,*>  CLR,P0S2 
160=15 
NUMB=20 

CALL  DISPLAY(IGO,NUMB-CLR) 
200   WRITE(*,*)  CLR,P0S2 

WRITE(*,*)'  THIS  SECTION  CHECKS  BODY  PANEL  FOR  AN  INTERNAL' 

WRITEt*,*)'  STATIC  PRESSURE  OF  150  PSIG.' 

WRITE(*,*) 
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IS=IS+1 

CALL  RETURN (CLR) 
WRITE (*,*!  CLR,P0S2 

WRITE(*.*)'  THIS  SECTION  CHECKS  TO  SEE  IF  BODY  PANEL' 
WRITER,  *)'  HAS  A  PENETRATION' 
WRITER*) 
CALL  RETURN (CLR) 
WRITE(*,*)  CLR.P0S2 

WRITE(*,'(A\)')'  DOES  THE  BODY  PANEL  HAVE  A  PENETRATION^)' 
READ (*, 310) ST 
310   FORtlAT(A) 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  '■)  THEN 
WRITE(*,*)CLR,P0S2 

WRITE(*,*)'  DO  YOU  WANT  TO  HAVE  COMPUTER  CALCULATE  ' 
WRITE(*,'(A\)')'  STIFFNESS  FACTORS?(Y>' 
READ  (*,310)ST1 
IF  (STl.EQ.'Y'.OR.STl.EQ.'Y'.OR.STl.EQ.'  ')  THEN 
IFLA6P=1 
NNN=IS 

CALL  PENET(CLR,P0S2, IFLAQP,STRFA,STRFB,A,B,EP,YP,HH,TOP 
l.NNN) 

WRITE(*,*) 

WRITE(*,*>  STRFA,STRFB 
ELSE 

WRITE(*,*)  CLR.P0S2 

WRITE**,*)'  YOU  MUST  ENTER  STRFA.STRFB' 
WRITE(*,*) 

READ(*,*)  STRFA,STRFB 
END  IF 
ELSE 
STRFA=1. 
STRFB=1. 
ENDIF 
100    WRITE(*,*)CLR,P0S2 

WRITE(*,*)'  THIS  SECTION  CHECKS  FOR  SUPPORTS' 
WRITE(*,*) 

WRITE(*,*)'  WHICH  TYPE  OF  SUPPORT  DOES  THE  PANEL  HAVE' 
WRITE**,*) 

WRITE(*,*)'   (1)  SIMPLE  SUPPORTS' 
WRITE**,*)'   (2)  BEND' 
WRITE**,*)'   (3)  WELDED  JOINTS' 
WRITE**,*!'   (4)  BOLTED  COVER' 
WRITE**,*) 

WRITE**,' (A\)')'  INPUT  YOUR  CHOICE — )' 
READ**, '(BN, 16)') IB 
IF  (IB.EQ.l)  THEN 
EJA=0. 
EJB=0. 

ELSEIF  (IB.EQ.2)  THEN 
EJA=1. 
EJB=1. 

ELSEIF ( IB. EQ. 3)  THEN 
WRITE**, *)CLR,P0S2 

WRITE**,*)'  YOU  MUST  ENTER  WELD  JOINT  EFFICIENCIES' 
WRITE**,*) 
CALL  RETURN (CLR) 
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4000  WRITER 
WRITEC* 

WRITE!* 
WRITE!* 
WRITE!* 
WRITE!* 
WRITE!* 
WRITE (* 
WRITE!* 
WRITE(* 
WRITE!* 
READ  (* 


)CLR,P0S2 

)'  DO  YOU  WANT  TO:' 


(1)  ENTER  THE  WELD  JOINT  EFFICIENCY' 

(2)  VIEW  THE  WELD  JOINT  EFFICIENCY  DATA  BASE' 
WHICH  SHOWS  THE  STATIC  JOINT  EFFICIENCIES' 
SHOWN  BY  CLASS  OF  WELD. ' 


)' 
)' 

)' 
*)' 
*) 
*) 

5 (A\) ' > '  INPUT  YOUR  CHOICE — > ' 
'(BN,I6)')I24 
IF  (I24.EQ.1)  THEN 
WRITE!*, *)CLR,  P0S2 


WRITE!*,*}'  ENTER  THE  WELD  JOINT  EFFICIENCIES;  EJA,  EJB' 
WRITE!*,*)'  EJA  SHOULD  BE  FOR  THE  SHORTER  SIDE  AND  EJB' 
WRITE!*,*)'  SHOULD  BE  FOR  THE  LONGER  SIDE.  YOU  SHOULD' 
WRITE!*,*)'  ENTER  THE  EFFICIENCIES  AS  A  PERCENTAGE.' 
WRITE!*,*)'  ENTER  EJA  FOLLOWED  BY  SPACE  AND  THEN' 
WRITE!*,*)'  EJB.' 

WRITE!*,*)'  AS  AN  EXAMPLE,  FOR  AN  EFFICIENCY  OF  20  V 
WRITE!*,*)'  YOU  SHOULD  ENTER  20.' 
READ!*,*!  EJA, EJB 
EJA=EJA/100. 
EJB=EJB/100. 
ELSE 

WRITE!*, *)CLR,P0S2 
CALL  WELDE(CLR) 
WRITE!*,*) 
CALL  RETURN (CLR) 
GO  TO  4000 
ENDIF 
4500       CONTINUE 

ELSEIFdB.EQ.4)  THEN 
WRITE!*,*)  CLR,P0S2 
WRITE!*, *)'  DO  YOU  WANT  TO:' 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 
WRITE!*,* 

WRITE!*,  '(A\}'P  INPUT  YOUR  CHOICE- — >' 
READ!*, '!BN, 16)') IBB 
WRITE!*,*)  CLR.P0S2 
IFflBB.EQ.D  THEN 

WRITE!*,*) ■  YOU  MUST  ENTER  JOINT  EFFICIENCIES' 
WRITE(*,*> 
WRITE!*,*)'  ENTER  EJA  FOLLOWED  BY  A  SPACE  AND  THEN  EJB' 


(1) 


ENTER  THE  JOINT  EFFICIENCIES  WHICH' 
ARE  THE  RATIO  OF  THE  MOMENTS  THAT' 
CAN  BE  RESTRICTED  BY  THE  BOLTS' 
TO  THE  FULLY  PLASTIC  MOMENTS  OF  THE' 
COVER. ' 


(2)  HAVE  THE  PROGRAM  CALCULATE  EJA  AND' 
EJB  FOR  THE  BOLTED  COVER. ' 
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WRITEC*,*)'  AS  A  PERCENTAGE  LESS  THAN  1007..' 
READ(*,*)EJA,EJB 
EJA=EJA*.01 
EJB=EJB*.01 
ELSE 

WRITEC*.*)1  YOU  HAVE  CHOSEN  TO  HAVE  THE  PROGRAM  ' 
WRITER,*)'  CALCULATE  THE  JOINT  EFFICIENCIES.' 
WRITEC*,*) 

WRITEC*,*)'  THE  PROGRAM  WILL  NEED  THE  FOLLOWING:' 
WRITEC*,*) 
WRITEC*,*)'  COVER  MATERIAL  YIELD  STRESS' 
WRITE(*,*)'  YIELD  LOAD  FOR  THE  BOLT,' 
WRITEC*,*)'  NUMBER  OF  BOLTS  ALONG  SIDE  A,  THE  SHORT  SIDE.' 
WRITEC*,*)'  NUMBER  OF  BOLTS  ALONG  SIDE  B.  THE  LONG  SIDE.' 
WRITEC*,*)'  COVER  THICKNESS.' 
WRITEC*.*)'  LENGTH  EA  FROM  FIGURE  6.' 
WRITEC*.*)'  LENGTH  EB  FROM  FIGURE  6.' 
WRITEC*,*) 
WRITEC*,*) 
CALL  RETURN (CLR) 
WRITEC*. *)CLR,POS2 
IFN=0 

IF(IFLAGP.EQ.l)  THEN 
WRITEC*,*)'  PREVIOUSLY  DETERMINED  VALUES  ARE:' 
WRITEC*,*) 

WRITEC*,*)'  COVER  MATERIAL  YIELD  STRESS:  ',YP 
WRITEC*,*)'  COVER  DIMENSIONS:  A=  ',A,'  B=  ',B 
WRITEC*,*)'  COVER  THICKNESS:  ',H 
WRITEC*,*) 

WRITEC*,' (A\)')!  ARE  THEY  CORRECT? <Y)' 
READ (*, 310) STT 

IFCSTT.EQ.'Y'.OR.STT.EQ.'Y'.OR.STT.EQ.'  ')  THEN 
CONTINUE 
ELSE 
IFN=1 
END  IF 
ELSE 
IFN=1 
ENDIF 

WRITEC*,*)  CLR.P0S2 
IFCIFN.EQ.l)  THEN 
WRITEC*,*)'  DO  YOU  WANT  TO:' 
WRITEC*,*) 
WRITEC*,*) 

WRITEC*,*)'    CD  ENTER  THE  YIELD  STRESS  IN  KSI.' 
WRITEC*,*)'    (2)  RETRIEVE  FROM  DATA  BASE.' 
WRITEC*,*) 
WRITEC*,*) 

WRITEC*,' (A\)')'  INPUT  YOUR  CHOICE- — )' 
READ (*,' CBN, 16)') 124 
IF (124. EG. 1)  THEN 
WRITEC*, *)CLR,P0S2 

WRITEC*,*)'  ENTER  THE  YIELD  STRESS  IN  KSI' 
READ(*,*)YP 
WRITEC*,*) 
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WRITE!*,*) '  ENTER  THE  THICKNESS  OF  THE  COVER  PLATE.1 
READ(*,*)  H 
ELSE 

CALL  RETRI  (CLR.P0S2,  YP,EP,PR,DEN.  IFLAGM,  IFLAGLTOP) 
WRITER,  *>CLR,P0S2 
CALL  RETURN (CLR! 

WRITE!*,*)'  ENTER  THE  THICKNESS  OF  THE  COVER  PLATE.' 
READ(*.*)H 
WRITE(*,*) 

UiR I TE ( * 7  * ) '  ENTER  DIMENSION  A  FOLLOWED  BY  A  SPACE' 
WRITEU,*)'  AND  THEN  DIMENSION  B.' 
HRITE(*i*J 
ENDIF 
ELSE 
ENDIF 

WRITE!*,*)'  ENTER  THE  YIELD  LOAD  FOR  THE  BOLT  IN  KS1.' 
READ(*,*)BYL 
WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  NUMBER  OF  BOLTS  ALONG  SIDE  A.' 
READ  v*,*)EAN 
WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  NUMBER  OF  BOLTS  ALONG  SIDE  B. ' 
READ!*,*)  EBN 
WRITE!*,*} 

WRITE!*,*)'  ENTER  THE  DIMENSION  EA.' 
READ!*,*)EA 
WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  DIMENSION  EB. ' 
READ!*,*)  EB 
WRITE!*,*) 

EJA=(8. *BYL*EA*EAN) / (3. *A* (H**2. ) *YP) 
EJB= (8. *BYL*EB*EBN) / (3. *B* (H**2. ) *YP) 
ENDIF 

WRITE!*, *)CLR,P0S2 

WRITE!*,*)'  THE  JOINT  EFFICIENCIES  ARE:' 
WRITE!*,*)'  EJA  =  ',  EJA,'    EJB  =  ',EJB 
IF(EJA.8T,1.)  THEN 

WRITE!*,*)'  EJA  CANNOT  BE  GREATER  THAN  1.' 
WRITE!*,*)'  YOU  SHOULD  RECALCULATE  OR  REENTER.' 
GO  TO  100 
ELSE 
ENDIF 

IF(EJB.6T.i.)  THEN 

WRITE!*,*)'  EJB  CANNOT  BE  GREATER  THAN  1.' 
WRITE!*,*)'  YOU  SHOULD  RECALCULATE  OR  REENTER.' 
GO  TO  100 
ELSE 
ENDIF 

WRITE!*, *)CLR,P0S2 
WRITEi*,*) 
WRITE!*,*) 

CALL  RETURN (CLR) 

ELSE 

WRITE!*, *>'  YOU  MUST  TYPE  1-4' 

CALL  RETURN (CLR) 
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GOTO  100 

ENDIF 

WRITEt*,*) 

CLR.POS2 

*THIS  SECTION  DETERMINES  C 

IF  (IFLAGP.EQ.l)  THEN 

CALL  RETURN (CLR) 

YP=YP*1000. 

EP=EP*1000. 

WRITEC*,*)' 

VALUES  OF  DIMENSION  A,  DIMENSION  8' 

WRITEt*,*)' 

THICKNESS  AND  YIELD  STRENGTH  OF  THE' 

WRITEt*,*)' 

PANEL  ARE  THE  FOLLOWING. ' 

WRITEt*,*! 

WRITE(*,*)' 

DIMENSION   A:    ',A 

WRITEt*,*)' 

DIMENSION   B:    '.B 

WRITEt*,*)' 

THICKNESS   H:    »,H 

WRITE(*,*)' 
WRITEt*,*) 

YIELD  STRENGTH  IN  KSI  :    '.YP/1000. 

IF2=0 

WRITEt*,' <A\)')'  ARE  THEY  CORRECT?«Y)' 
READ (*, 310)  STT 

IF<STT.E8.'Y!.0R.STT.EQ.'Y',0R.8TT.EQ.'  ')  THEN 
CONTINUE 
ELSE 
IF2=1 
ENDIF 
ELSE 
IF2=1 
ENDIF 
3030       IF(IF2.EQ.l)  THEN 
WRITE (*,*) CLR, P0S2 

WRITEt*,*)'  YOU  MUST  ENTER  DIMENSION  A,  B,  AND  THICKNESS.' 
WRITEt*,*)'  DIMENSION  A  MUST  BE  SHORTER  THAN  DIMENSION  B.' 
WRITE(*,*) 

WRITE (*,*)'  ENTER  A  FOLLOWED  BY  A  SPACE  AND  THEN  B. ' 
READ(*,*)  A,B 
WRITEt*,*) 

WRITEt*,*)'  ENTER  THE  THICKNESS  OF  THE  PANEL.' 
READ(*,*)H 
WRITEt*,*)' 
WRITEt*,*)' 
WRITEt*,*)' 
WRITEt*,*) 

WRITEt*,*)'  DIMENSION   A 
WRITEt*,*)'  DIMENSION   B 
WRITEt*,*)'  THICKNESS   H 
WRITEt*,*) 
IF2=0 

WRITE'.*,'  (A\)')'  ARE  THEY  CORRECT? (Y)' 
READ (*, 310)  STT 

IFtSTT.EQ.'Y'.OR.STT.EQ.'Y'.OR.STT.EQ.'  '<  THEN 
CONTINUE 
ELSE 
IF2=1 

60  TO  3030 
ENDIF 


VALUES  OF  DIMENSION  A 
AND  THICKNESS  OF  THE' 
PANEL  ARE  THE  FOLLOWING 


',A 
',B 
',H 


DIMENSION  B' 
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ELSE 
IF2=1 

ENDIF 

WRITE (*,*!CLR,  P0S2 

WRITE!*, *)'  DO  YOU  WANT  TO:' 

WRITE**,*) 

WRITE(*,*) 

WRITE(*,*)'  (1)  ENTER  THE  YIELD  STRENGTH  IN  KSI.' 

WRITE(*,*)'  (2)  RETRIEVE  FROM  DATA  BASE1 

WRITE**,*) 

WRITE(*,*I 

WRITE (*,'(A\)')'  INPUT  YOUR  CHOICE >' 

READ{*,'{BN,I6)')I23 
WRITE(*,*)  CLR,P0S2 
IF  (I23.EQ.1)  THEN 

WRITE**,*)'  ENTER  THE  YIELD  STRENGTH  IN  KSI' 
READ(*,*)YP 
ELSE 
CALL  RETRI  <CLR,P0S2. YP.EP.PR.DEN,  IFLAGM,  IFLAMLTOP) 

ENDIF 
ELSE 
ENDIF 

IF(IFLAGP.EQ.l)  THEN 
CONTINUE 
ELSE 

YP=YP*i000. 
EF-EP*1000. 
ENDIF 
GX=A/B 

IF  (8X.ST.1.)  THEN 

WRITE**,*)  '  DIMENSION  A  MUST  BE  LESS  THAN  DIMENSION  B.' 
WRITER*) 

WRITES*,*)  '  THE  VALUES  OF  A  AND  B  WILL  BE  SWITCHED.' 
AF=A 
A=B 
B=A 

WRITE**,*) 

WRITE**,*/'  NEW  VALUE  OF  A:  ',A,'  NEW  VALUE  OF  B:   ', 
ELSE 
ENDIF 
GX  =  A/B 
WRITE(*,*) 
WRITE**,*) 
CALL  RETURN (CLR) 
IF  (8X.LE..5)  THEN 

C=1.5 
ELSEIF  (6X.6E..5.ANB.SX.LE..6)  THEN 

C=1.5 
ELSEIF  iGX.GE.,6,AND.3X,LE..7)  THEN 

C=1.29 
ELSEIF  (SX.6E..7.AND.SX.LE..8)  THEN 

C=1.17 
ELSEIF  (SX.8E..8.AND.6X.LE..9)  THEN 
C=1.07 
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ELSEIF  (GX.GE..9.AND.GX.LE.1.0)  THEN 

C=1.0  2 
ELSE 

C=1.5 
ENDIF 
3031   WRITE  (*,*>'  THE  VALUE  FOR  STRFA  USED  IS  ', STRFA 
WRITE  (*.*)'  THE  VALUE  FOR  STRFB  USED  IS  ', STRFB 
WRITE  (*,»)'  THE  VALUE  FOR  EJA   USED  IS  ',EJA 
WRITE  (*,*)'  THE  VALUE  FOR  EJB   USED  IS  !,EJB 
WRITE(*.MA\)')'  ARE  THEY  CORRECT?(Y)' 
READ (*, 310)  STT 

IF(STT.EQ.'Y'.OR.STT.EQ.'Y'.OR.STT.EQ.'  ')  THEN 
CONTINUE 
ELSE 
WRITE(*,*)  'ENTER  STRFA' 
READ  (*,*)  STRFA 
WRITE(*,*)  'STRFA  =  ',  STRFA 
WRITE(*,*)  'ENTER  STRFB' 
READ(*i*)  STRFB 
WRITE(*,*)  'STRFB  =  ', STRFB 
WRITE(*,*)  'ENTER  EJA' 
READ(*,*)  EJA 
WRITE  (*,*)  'EJA  =  ',EJA 
WRITE'*,*)  'ENTER  EJB' 
READ(*,»)  EJB 
WRITE(*,*)'EJB  =  ',EJB 
60  TO  3031 
♦THIS  SECTION  CALCULATES  PA 
PM0=,25*YP*(H**2) 

RM=12. *PMO* (STRFA+EJA+C* (STRFB+EJB) ) 
PA=RM/(A*B) 
IFLAGW=1 
IT=IT+1 
WRITE(23,*) 
WRITE (23,*) 

WRITE<23,*>'  RESULTS  OF  PANEL  SCREENING' 
WRITE (23,*) 

WRITE(23,*)'  STATIC  PRESSURE  SCREENING' 
WRITE(*,*> 

WRITE(23,*)'  PANEL  NUMBER:  MS 
PAN(IS,1)=RM 
PAN(IS,2)=H/12. 
PAN(IS,3)=A*B/144. 
IRUG=1 


WRITE(23,*) 
WRITE (23,*) 
WRITE (23,*) 
WRITE(23,*) 
WRITE(23,*) 
WRITE(23,*) 
l.YP/1000. 

WRITE(23,*) 
WRITE(23,*) 
WRITE(23,*) 
WRITE(23,*) 


DIMENSIONS  OF  PANEL: ' 

A=  ',A,'  B=  ',B,'   THICKNESS=   ',H 

THE  VALUE  OF  A/B  USED  IN  CALCULATION  IS  ',GX 

THE  VALUE  OF  C  USED  IS  ',C 

THE  VALUE  OF  YIELD  STRENGTH  IN  KSI  USED  IS  ' 

THE  VALUE  FOR  STRFA  USED  IS  ', STRFA 

THE  VALUE  FOR  STRFB  USED  IS  ', STRFB 

THE  VALUE  FOR  EJA  USED  IS  \EJA 

THE  VALUE  FOR  EJB  USED  IS  ',EJB 
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THE  VALUE  DETERMINED  FOR  RM  IS  ',RM 
THE  VALUE  DETERMINED  FOR  PA  IS  \PA 


WRITE (23,* 
WRITE(23,* 
WRITE (23,*) 
IF (PA. LE. 150.)  THEN 

WRITE(*,*)'  BODY  PANEL  SHOULD  BE  REJECTED  ' 
WRITEi*,*)'  PA=   ',PA 

WRITE(23,*)'  THE  BODY  PANEL  SHOULD  BE  REJECTED 
WRITE(23,*)'  THE  CALCULATED  STATIC  PRESSURE  PA' 
WRITE(23,*)'  IS  LESS  THAN  150  PSIG.' 
WRITE(23,*)'  PA=   ',PA 
ELSE 

WRITE (*,*)'  BODY  PANEL  OK   ' 
WR1TE(*,*)'  PA=   ',PA 
WRITE (23,*)'  BODY  PANEL  OK  ' 
WRITE(23,*)'  PA=   ',PA 
ENDIF 
♦CHECK  ANOTHER  PANEL 

CALL  RETURN (CCR) 


400 


2001 


WRITE ( 
WRITE ( 
WRITE ( 
WRITE( 
WRITE ( 
WRITE( 
WRITE( 
WRITE( 
WRITE ( 
WRITE ( 
WRITE( 
WRITE ( 
WRITE( 
WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 


,*)  CLR,P0S2 

,*)'  DO  YOU  WANT  TO 

,*! 


(11 


',*)' 
>,*)' 
',*)' 
',*)' 
-,*)' 
',*)' 
■,*)' 
>,*)' 
>,*) 
>,*)' 
>,*) 
>,*)' 
',*) 
',*) 

,'(A\)')' 
READ(*,'(BN,I6)') 
IF  (IG.EQ.l)  THEN 
GO  TO  200 

ELSEIF  (IG.EQ.2)  THEN 
GO  TO  2001 

ELSEIF  (IG.EQ.3)  THEN 
RETURN 
ELSE 

WRITE(*,*)' 
WRITE(*,*)' 
WRITE(*,*) 
WRITE(*,*) 
CALL  RETURN (CLR) 
GO  TO  400 
ENDIF 
WRITE(*,*)CLR,P0S2 
CONTINUE 
160=16 
NUMB=20 


CHECK  ANOTHER  PANEL  FOR  STATIC 
PRESSURE. ' 

NOTE:  YOU  SHOULD  ONLY  CHOOSE  THIS' 
OPTION  IF  YOU  ARE  NOT  PLANNING' 
TC  CHECK  THE  PREVIOUS  PANEL' 
FOR  DYNAMIC  PRESSURE.  YOU  WILL' 
NEED  THE  PA  VALUE  TO  CHECK  FOR' 
DYNAMIC  PRESSURE, ' 


(2)  CHECK  PANEL  FOR  DYNAMIC  PRESSURE' 


(3)  RETURN  TO  MAIN  MENUE' 


INPUT  YOUR  CHOICE — )' 
IG 


YOU  ONLY  HAVE  3  CHOICES' 
TRY  AGAIN' 
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CALL  DISPLAY (160, NUMB, CLR) 
WRITE(*,*)CLR,P0S2 
************************ 

WRITE(*,*)'  YOU  MUST  ENTER  THE  NUMBER  AS  AN  INTEGER  FROM' 
WRITER, *)'  FIGURE  4   —ENTITLED' 
WRITE(*,*)'  BOUNDARY  CONDITIONS  FOR  RECTANGULAR  PLATES' 
WRITE(*,*)'  (INTEGER  FROM  1  TO  5)' 

WRITE<*,*)'  WHICH  CORRESPONDS  TO  THE  BOUNDARY  CONDITIONS' 
WRITE(*,*)'  FOR  YOUR  PANEL' 
READ(*,*)IR 
5002   6X=A/B 

WRITER,  *)'  THE  RATIO  OF  B/A  USED  IS' 

WRITE(*,*)  'B/A  =  :M./GX 

AT(1,1)=19.74 

AT (1,2) =32. 08 

AT(1,3)=71.56 

AT(2,1)=11.68 

AT(2,2)=13.71 

AT(2,3)=18.8 

AT(3,1)=28.95 

AT (3, 2) =56. 35 

AT(3,3)=145.5 

AT(4,1)=24.02 

AT (4, 2) =26. 73 

AT (4, 3) =37. 66 

AT (5,1) =35. 99 

AT (5, 2) =60. 77 

AT(5,3)=147.80 

GX=1./GX 

IF  (GX.GE.1.0.AND.GX.LT.1.5)  THEN 

IC=1 
ELSEIF  (GX.GE.1.5.AND.GX.LT.2.5)  THEN 

IC=2 
ELSEIF  (GX.GE.2.5)  THEN 

IC=3 

ELSE 

WRITE(*,*)  '  A  MUST  BE  LESS  THAN  B' 

WRITE (*,*)  '  REENTER  A  FOLLOWED  BY  B' 

READ  (*,*)  A,B 

SO  TO  5002 

ENDIF 

GX=1./GX 

IF  (IC.EQ.l)  THEN 

XX=AT(IR,IC) 

ELSEIF  (IC.EQ.2)  THEN 

XX=AT(IR,IC) 

ELSE 

XX=AT(IR,IC) 

ENDIF 

WRITEt*,*)'  THE  VALUE  OBTAINED  FOR  THE' 
WRITE(*,*)'  FUNDAMENTAL  CIRCULAR  FREQUENCY  FOR' 
WRITE(*,*)'  THE  RECTANGULAR  PLATES' 
WRITE(*,») 

WRITE(*,*)1  WHICH  IS  REFERENCED  AS  XX' 
WRITE(*,*)'  FOR  THE  BOUNDARY  CONDI TON' 
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WRITEC*,*)'  NUMBER  MR,'  IS   ',XX 
WRITE(*,*) 

WRITE(*,*)'  ENTER  THE  ELASTIC  MODULUS  IN  KSI' 
READ  (*,*)  EP 
EP=EP*1000. 

WRITE  (*,*)  '  ENTER  POISSONS  RATIO  AS  A  DECIMAL' 
READ  (*,*)  PR 

WRITE  (*,*)  '  ENTER  THE  DENSITY  IN  LB/IN3' 
READ(*,*)  DEN 
3000   CONTINUE 

WRITE(*,*>'  DETERMINED  VALUES  ARE:' 

WRITE**,*) 

WRITE(*,*)'  YIELD  STRESS  IN  KSI  =     ',YP/1000. 

WRITE(*,*)'  ELASTIC  MODULUS  IN  KSI=  '.EP/1000. 

WRITE(*,*)'  POISSONS  RATIO=   ',PR 

WRITE(*,*)'  DENSITY=        ',DEN 

WRITE(*,*)'  A=   ',A,'   B=   SB 

WRITE(*,*)'  THICKNESS=   ',H 

WRITE  (*,*)'  B/A  =  ',  l./GX 

WRITE  (*,*)'  XX  =  ',  XX 

WRITE(*,*) 

WRITE(*,'(A\)')'  ARE  THEY  CORRECT?(Y)' 

READ (*, 310) STT 

IF(STT,EQ.'y'.OR.STT.EQ.'Y'.OR.STT.EG.'  ')  THEN 

CONTINUE 

ELSE 

WRITE(*,*)'  ENTER  THE  YIELD  STRESS  IN  KSI' 

READ(*,*)YP 

YP=YP*1000. 

WRITE(*,*) 

WRITE(*,*)'  ENTER  THE  ELASTIC  MODULUS  IN  KSI' 

READ(*,*)EP 

EP=EP*i000. 

WRITE(*,*) 

WRITE(*,*)'  ENTER  POISSONS  RATIO  AS  A  DECIMAL' 

READ(*,*)PR 

WRITE(*,*) 

WRITE(*,*)'  ENTER  DENSITY  IN  LB/IN3' 

READ (*,*) DEN 

WRITE!*,*)'  ENTER  A/B' 

READ(*,*)GX 

WRITE!*,*)'  ENTER  XX' 

READ(*,*)  XX 

60  TO  3000 

ENDIF 

D=EP* !H**3. ) / ( 12. * ( 1 . -PR) **2. ) 

W=(XX)/((B**2.)*((DEN*H/D)**.5)) 

T=2000,*3.14/W 

WRITE  (*,*)  'D  =  ',  D  ,  'W  =  ',  W,  'T  =  ',  T 

WRITE(*,*)'  A  REASONABLE  UPPER  LIMIT  FOR  THE  EXPLOSION' 

WRITE(*,*)'  PRESSURE,  WITHOUT  PRESSURE  PILING,  IS' 

WRITE!*,*)'  CONSIDERED  TO  BE  100  PSIG.' 

WRITE(*,*) 

WRITE(*,'(A\)')'  DO  YOU  WANT  TO  USE  100  PSIG?<Y)' 

READ(*,310)STT 
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IF  (STT.EQ.'Y'.OR.STT.EQ.'y'.OR.STT.EQ.'  ')  THEN 
PMAX=100. 
ELSE 

WRITEf*,*)'  YOU  MUST  ENTER  PMAX  IN  PS 16.' 
WRITE(*,*/'  ENTER  PMAX' 
READ(*,*)  PMAX 
ENDIF 
1050   CONTINUE 

WRITER.*)'  THE  CURRENT  VOLUME  OF  YOUR  ENCLOSURE  IS:' 
WRITEt*,*)'  VOLUME  IN  CU  INCHES:    ', VOLUME 
WRITE(*,*) 

WRITE(*,'(A\)')'  IS  THAT  CORRECTLY)' 
READ(*,310)STT 

IF  (STT.EQ.'Y'.OR.STT.EQ.'y'.OR.STT.EQ.'  '!  THEN 
CONTINUE 
ELSE 

WRITER, *)'  YOU  MUST  ENTER  THE  VOLUME  IN  CUBIC  INCHES.' 
WRITE<*,*)'  ENTER  VOLUME' 
READ <*,*) VOLUME 
GO  TO  4050 
ENDIF 

V=V0LUME/1728. 

TMIN=(PMAX/4.72)*((V**.333)) 
TX=TMIN/T 
**************** 

WRITE  (*,*)CLR,P0S2 

CALL  FIG4(TX,DLF) 

WRITE(*,*) 

WRITE(*,*) 

CALL  RETURN (CLR) 

WRITE  (*,*)  '  DLF  =  ',  DLF 

WRITE(*,*)'  VOLUME=  ',  VOLUME,'   TMIN=   ',TMIN,'  T=   M 

WRITE!*,*)'  PMAX=   ',  PMAX 

WRITE(*,*)'  (TMIN/T)=   ',TX 

WRITE(*,*)'  D=  ',D 

WRITE(*,*) 

WRITER, ♦) 

WRITE(*,*) 

WRITE(*,*) 

PEQ=PMAX*DLF 

WRITEf*,*)'  PEQ=   ',PEQ 

WRITE(*,*) 

WRITE(23,*) 

WRITE(23,*) 

WRITE(23,*) 

WRITE(23,*)'  RESULTS  OF  DYNAMIC  SCREENING' 

WRITE(23,*) 

WRITE(23,*)'  PANEL  NUMBER:  SIS 

WRITE(23,*) 

WRITE(23,*)'  VOLUME=  ',  VOLUME,'   TMIN=   ',TMIN,'  T=   ',T 

WRITE(23,*)'  PMAX=   ',  PMAX 

WRITE(23,*)'  BOUNDARY  CONDITION  NUMBER  FROM  FIGURE  4  IS  MR 

WRITE(23,*)'  THE  VALUE  OF  A/B  USED  IS  ',  GX 

WRITE (23,*)'  THE  VALUE  FROM  TABLE  5.2  ' 

WRITE(23,*)'  FROM  REFERENCE  3-ESTIMATING  RECTANGULAR' 
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WRITE<23,*)'  XP  ENCLOSURE  PERFORMANCE,  FINAL  ENGINEERING' 

WRITE(23,*)'  REPORT ' 

WRITE (23,*)'  IS  ',  XX 

WRITE  (23,*)'  PA  FROM  STATIC  SCREENING  =  ',  PA 

WRITE(23,*)'  (TMIN/T)=   MX 

WRITE(23,*)'  D=  ',D 

WRITE  (23,*)  '  DLF  =  ',  DLF 

WRITE(23,*)'  DIMENSIONS  OF  PANEL  ARE:' 

WR1TE(23,*)'  A=   *,A,'  B=   ',B,'  THICKNESS=  ',H 

WRITE(23,*)'  YIELD  STRENGTH  USED  IN  KSI  IS   ',YP/1000. 

WRITE(23,*)'  ELASTIC  MODULUS  USED  IN  KSI  IS  ',EP/1000. 

WRITE (23,*)'  DENSITY  IS  ',  DEN 

WRITE(23,*)'  POISSONS  RATIO  IS  ',  PR 

WRITE(23,*)'  VOLUME  USED  IS        ', VOLUME 
PTEST=PA*.7 
IRPTEST.LT. PEQ)  THEN 

WRITE(23,*)'  THE  PANEL  FAILED  THE  DYNAMIC  TEST' 

WRITE(23,*)'  PEQ  IS  GREATER  THAN  .7*PA' 
ELSE 

WRITE (23,*) 
WRITE(23,*) 

WRITE(23,*)'  THE  PANEL  IS  OK' 

WRITE(23,*)'  PEQ  IS  LESS  THAN  .7  ♦  PA' 
ENDIF 

WRITE (23,*) 

WRITE(23,*)'  PEQ=  SPEG,'   .7*PA=   ',PTEST 

5001   WRITE(*,«)  '  YOU  ARE  FINISHED' 
RETURN 

C    *  * 

C    *         SUBROUTINE  COVER  (CC.FOR)  * 

C    *  * 

C    444******************************************************** 
$LARGE 

SUBROUTINE  COVER (ID, TOP, IFLAGM, CLR, P0S2) 

CHARACTER*64  FNAME,ID,FLM2(33) 

REAL  MU,MX,MY,MXY,IB,KS,LB,K1,K2,K11,K12,K13,K,KSV,MXT,MYT 

DIMENSION  AA(25,25),EPSILON(25,25),X(10),Y(10),W(10),WX(10), 
WXX(10),WXXX(10),WY(10),WYY(10),WYYY(10),WXY(10), 
WXXY ( 10) , WYYX ( 10) , SIGMAX ( 10) , SIGMAY ( 10) , TAUXY ( 10) , 
TAUXZ(10),TAUYZ(10),XB(40),YB(40),C(25,25),K(200), 
CSV (25, 25) , KSV (200) , WB (40) , THETABX (40) , 
THETACY(40) ,SI6MABX(40) ,SIGMACY(40) ,TAUBXZ(40) , 
TAUBYZ(40),SIGMAAX(40),WORK(200),FU2),THICK(6) 

CHARACTER*1  POS2(10),CLR(4),CR,ESC,ST,POS1(9),REV(4),NORM(7) 

CHARACTER*1  TEXT(80),TEXTX(18,80),TOP(9) 

DIMENSION  6(11*3) 

WRITE  (*,*)  CLR.P0S2 

WRITE  (*,'(A\)')'  DO  YOU  WANT  TO  VIEW  INFORMATION? (Y)' 

READ  (*,310)  ST 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 

GO  TO  506 

ELSE 

60  TO  505 

ENDIF 
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506  WRITE!*,*)  CLR,P0S2 
IG0=8 

NUMB=20 

CALL  BISPLAY!IGO,NUMB,CLR) 
IFLA8A=0 
IFLAGB=0 
600  WRITE**.*)  CLR,P0S2 
CONTINUE 
160=9 
NUi1B=20 

CALL  DISPLAY (160. NUMB, CLR) 
505  CONTINUE 

WRITE!*,*)  CLR,P0S2 

IF  (IFLAGA.EQ.0)  THEN 

WRITE(*,*)'  YOU  WILL  BE  ASKED  TO  ENTER  THE  DIMENSIONS  OF  THE' 

WRITE!*,*)'  COVER.' 

WRITE(*,*) 

WRITE(*,' (A\)')'  DO  YOU  WANT  TO  VIEW  PREVIOUS  SCREEN?(Y>' 

READ!*, 310)  ST 

IF  (ST.EQ.'Y'.OR  .ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 

60  TO  600 

ENDIF 

WRITE!*,*)  CLR,P0S2 

WRITE!*,*)'  ENTER  THE  DIMENSIONS  A  FOLLOWED  BY  B' 
WRITE!*, ♦)'  OF  PLATE  IN  INCHES:' 
WRITE!*,*) 
READ!*,*)  A,B 
WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  THICKNESS  OF  COVER  IN  INCHES:' 
WRITE!*,*) 
READ!*,*)  H 
WRITE!*,*) 
IFLAGA=1 
ENDIF 

WRITE!*,*) CLR, P0S2 

WRITE!*,' (A\)')'  DO  YOU  WANT  TO  ANALYZE  A  COVER  WITH  B0LT3?<Y>' 
READ!*, 310) ST 

IF  (ST.EQ.'Y'.OR. ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 
60  TO  507 
ELSE 
NB=0 
GO  TO  508 

ENDIF 

507  CONTINUE 

IF(IFLAGB.EQ.0)  THEN 

WRITE!*,*)  CLR,P0S2 

WRITE (*,*)'     YOU  WILL  BE  ASKED  TO  INPUT  BOLT' 

WRITE (*,♦)'  POSITIONS  AND  BOLT  SPECIFICATIONS' 

WRITE!*,*)'  FOR  COVER  ANALYSIS.' 

WRITE!*,*) 

WRITE (♦,*)'  THE  BOLT  POSITIONS  ARE  DESCRIBED  ON  THE  PREVIOUS' 

WRITE!*, ♦)'  SCREEN.' 

WRITE!*,*) 

WRITE!*,' (A\)')'  DO  YOU  WANT  TO  VIEW  IT  AGAIN?<Y>' 

READ  (*,310)  ST 


43 


IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.1  ')  THEN 
GO  TO  600 
ENDIF 
651   WRITE(*,*)CLR,P0S2 

WRITE(*,*)*  DO  YOU  WANT  TO:' 

WRITER,*) 

WRITE**,*)1    (1)  ENTER  YOUR  OWN  BOLT  CONFIGURATION.' 

WRITE!*,*)'    (2)  USE  THE  12  BOLT  CONFIGURATION' 

WRITE!*.*)'       DESCRIBED  ON  THE  FOLLOWING  SCREEN.' 

WRITE!*,*)'       (IF  YOU  CHOOSE  THIS  OPTION,  DIMENSIONS' 

WRITE!*,*)'       A,  B,  AND  THICKNESS  OF  COVER  WILL  BE' 

WRITE!*,*)'       CHANGED  TO  A=13.5,  B=10.75,  THICKNESS=.5' 

WRITE!*,*) 

WRITE!*,*) 

WRITE !*,'(A\)')'  INPUT  YOUR  CHOICE — >' 

READ!*, '(BN,  16)') 124 

IF  (I24.EQ.1)  THEN 

GO  TO  650 

ELSEIF  (I24.EQ.2)  THEN 

A=13.5 

B=10.75 

H=.5 

WRITE (*,*)  CLR.TOP 

150=11 

NUMB=20 

CALL  DISPLAY ( 160, NUMB, CLR) 

N8=12 

XB!1)=.45 

YB(1)=.45 

XB(2)=4.65 

YB(2)=.45 

XB(3)=8.85 

YB(3)=.45 

XB(4)=13.05 

Y8<4)=.45 

XB!5>=.45 

YBi5) =3.733 

XB(6!=13.05 

YB(6)=3.733 

XB(7)=.45 

YB(7)=7.016 

XB(8)=13.05 

YB(Si=7.016 

XB<9)=.45 

YB(9)=10.3 

XB<10)=4.65 

YB!10)=10.3 

XB(11)=B.S5 

YB(ll)=i0.3 

XB(12)=13.05 

YB!12)=10.3 

GO  TO  610 
ELSE 

WRITE!*-,*)'  YOU  ONLY  HAVE  TWO  CHOICES' 
30  TO  651 
ENDIF 
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650   WRITE (*,*)  CLR.P0S2 

WRITE!*,*)'  ENTER  THE  NUMBER  OF  BOLTS:  ' 

READ  <*,*)  NB 

WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  POSITIONS  OF  BOLTS  AS  DESCRIBED  ON  ' 

WRITER, *)'  PREVIOUS  SCREEN.  ENTER  THEM  AS  X,Y  COORDINATES' 

WRITE!*,*)'  SEPARATED  BY  A  SPACE  BETWEEN  X  AND  Y  AND  ONE  BOLT' 

WRITE!*,*)'  COORDINATE  PER  LINE  FOLLOWED  BY  A  RETURN.' 

WRITE!*,*) 

DO  601  IDD=1,NB 

READ!*,*)  XB(IDB),YB!IDD) 
601   CONTINUE 

IFLAGB=1 

ENDIF 

WRITE!*,*) 

WRITE!*,*) 

CALL  RETURN (CLR) 
WRITE!*,*) 
508  CONTINUE 

WRITE!*,*)  CLR.P0S2 
CALL  RETURN !CLR) 
610  WRITE!*,*)  CLR,P0S2 

WRITE!*,*)'     YOU  MUST  ENTER  THE  ELASTIC  MODULUS,  YIELD' 
WRITE!*,*)'  STRENGTH  AND  POISSONS  RATIO  OF  THE  COVER 
MATERIAL  OR' 

WRITE!*,*)'  RETRIEVE  IT  FROM  THE  ACCEPTABLE  MATERIALS' 
WRITE!*,*)'  DATA  BASE.' 
WRITE!*,*) 
WRITE!*,*) 
CALL  RETURN (CLR) 
510  WRITE!*,*) CLR, P0S2 

WRITE!*,*)'  DO  YOU  WANT  TO:' 

WRITE!*,*) 

WRITE!*,*) 

WRITE!*,*)'  (1)  ENTER  THE  ELASTIC  MODULUS  , YIELD' 

WRITE!*,*)'     STRENGTH  AND  POISSONS  RATIO  ' 

WRITE!*,*)'  (2)  RETRIEVE  THOSE  VALUES  FROM' 

WRITE!*,*) '     ACCEPTABLE  MATERIALS  DATA  BASE' 

WRITE!*,*) 

WRITE!*,*) 

WRITE (*,'(A\)')'  INPUT  YOUR  CHOICE — >' 

READ  (*,'(BN,I6)')I222 

IF  (I222.EQ.1)  THEN 

WRITE!*,*)  CLR,P0S2 

WRITE!*,*)'  ENTER  THE  ELASTIC  MODULUS  OF  COVER  IN  K5I' 

WRITE!*,*) 

WRITE!*,*)'  FOR  EXAMPLE,  THE  ELASTIC  MODULUS  OF  A-36  STEEL' 

WRITE!*, ♦)'  WOULD  BE  ENTERED  AS  29000.' 

WRITE!*,*) 

READ  !*,*)  EP 

WRITE!*,*)  CLR,  P0S2 

WRITE!*,*)'  ENTER  THE  VIELD  STRENGTH  OF  THE  COVER  IN  KSI' 

WRITE!*, ♦) 

READ(*,*)YP 

WRITE!*,*) 
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WRITE!*,*)'  ENTER  POISSONS  RATIO  OF  COVER  IN  DECIMAL  FORMAT:' 

WRITE**,*) 

READ(*,»)  MU 

WRITE!*,*)  CLR,P0S2 

ELSE  IF  (I222.EQ.2)  THEN 

CALL  RETRI (CLR, P0S2, YP, EP, PR, DEN, IFLAGM, IFLAMI , TOP) 

MU=PR 

WRITE(*,*) 

CALL  RETURN (CLR) 

ELSE 

WRITE!*,*)'  YOU  ONLY  HAVE  TWO  CHOICES' 

WRITE!*,*) 

CALL  RETURN (CLR) 

60  TO  510 

ENDIF 

EP=EP*1000. 

YP=YP*1000. 

WRITE!*,*) CLR,  P0S2 

WRITE!*,*) '  ENCLOSURE  IDENTIFICATION:   ' 
WRITE!*, 315) ID 
315  FORMAT (1X.64A) 

WRITE!*,  *) ' ' 

WRITE!*,*) 
WRITE!*,*) 
WRITE!*, 9030)  A,B,H 
WRITE!*,*) 

WRITE!*,*)'  COVER  PLATE  PROPERTIES' 
WRITE!*, 9050)  EP,MU,YP 
WRITE!*,*) 
CALL  RETURN (CLR) 
IF!NB.EQ.0)THEN 
BO  TO  316 
ELSE 

CONTINUE 
ENDIF 
317  WRITE!*,*)  CLR,P0S2 

WRITE!*,*)'     YOU  MUST  ENTER  BOLT  PROPERTIES  AND  DIMENSIONS' 
WRITE (*,♦)'  AS  DESCRIBED  ON  NEXT  SCREEN  OR  USE  VALUES  FROM' 
WRITE!*,*) '  DATA  BASE.' 
WRITE!*,*) 
WRITE!*,*) 

WRITE!*,' (A\)')'  DO  YOU  WANT  TO  VIEW  INFORMATION?(Y)' 
READ (♦, 310) ST 
319  CONTINUE 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 

WRITE!*,*)  CLR.TOP 

160=10 

NUMB=20 

CALL  DISPLAY (160, NUMB, CLR) 

WRITE(*-*)CLR,TOP 

ELSE 

CONTINUE 

ENDIF 
318   WRITE(*,*)CLR,P0S2 

WRITE (♦,*)'  DO  YOU  WANT  TO:' 
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(1)  ENTER  THE  BOLT  PROPERTIES  AND  DIMENSIONS' 

(2)  USE  THE  VALUES  FOR  1/2  INCH  CLASS  FIVE  BOLTS' 
LISTED  ON  THE  PREVIOUS  SCREEN. ' 
VIEW  PREVIOUS  SCREEN. ' 


(3/ 


,'(A\)')' 
MBN,I6)')I23 
.EQ.l)  THEN 
CLR.TOP 


INPUT  YOUR  CHOICE >' 


)' 
)' 

)' 
EB 
) 
)' 


ENTER  THE  ELASTIC  MODULUS  OF  THE  BOLT  IN  KSI' 
FOR  EXAMPLE,  THE  ELASTIC  MODULUS  OF  STEEL' 
WOULD  BE  ENTERED  AS  30000. ' 


ENTER  POISSONS  RATIO  IN  DECIMAL  FORMAT* 


PRR 


(2 


♦  U+PRR): 

) 
)' 

)DB 
) 

)' 
)' 


ENTER  THE  DIAMETER  OF  THE  BOLT  IN  INCHES. ' 


WRITER, *) 

WRITE'! 

WRITE ( 

WRITE! 

WRITE! 

WRITE! 

WRITE! 

WRITE! 

READ!* 

IF  (12 

WRITE! 

WRITE! 

WRITE! 

WRITE! 

READ!* 

WRITE! 

WRITE! 

READ  ( 

EB=EB* 

GB=EB/ 

WRITE! 

WRITE! 

READ  ( 

WRITE! 

WRITE! 

WRITE! 

WRITE! 

WRITE! 

READ!* 

WRITE! 

WRITE! 

WRITE! 

WRITE! 

READ!* 

WRITE! 

WRITE! 

WRITE! 

WRITE! 

READ!*,*) 

WRITE!*,*) 

CALL  RETURN (CLR) 

AB=3.14*(DB/2.)**: 

ELSEIF  (I23.EQ.2) 

EB=30000000. 

IB=. 003068 

KS=1.333 

AB=.1963 

GB=11J 

LB= 

DB=.5 

ELSEIF  (I23.EQ.3)  THEN 

ST='Y' 

GO  TO  319 

ELSE 

WRITE!*,*)'  YOU  ONLY  HAVE  3  CHOICES.' 


ENTER  THE  AREA  MOMENT  OF  INERTIA  ABOUT  ANY' 
DIAMETER.  FOR  EXAMPLE,  FOR  A  1/2  INCH  BOLT,  YOU' 
WOULD  ENTER  .003068.  FOR  OTHER  BOLTS,  THIS' 
INFORMATION  IS  FOUND  IN  MACHINERYS  HANDBOOK.' 


IB 


'     ENTER  THE  SHEAR  SHAPE  FACTOR  KS.  <  KS  =4/3' 
'  FOR  A  CIRCULAR  CROSS  SECTION.  ENTER  1.333  IF' 
'  YOUR  BOLTS  HAVE  A  CIRCULAR  CROSS  SECTION. ' 
KS 

'     ENTER  THE  THREAD  ENGAGEMENT  OF  THE  BOLT.' 
'  IN  INCHES.  A  TYPICAL  VALUE  FOR  A  1/2  INCH  BOLT' 
'  WOULD  BE  .B5  INCH.' 
LB 


THEN 


,85 
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GO  TO  318 
ENDIF 
316  NN=3 

CALL  RETURN (CLR) 

WRITE (*,*)'     STRESSES  IN  THE  COVER  OR  BOLTED  COVER  SHOULD  ' 

WRITE!*,*)'  BE  CALCULATED  FOR  A  150  PSI6  STATIC  PRESSURE' 

WRITE!*,*)'  USING  THE  ALGORITHM  IN  THIS  SECTION.' 

WRITE!*,*) 

WRITE!*,*)'    HOWEVER,  YOU  WILL  BE  PROVIDED  AN  OPPORTUNITY' 

WRITE!*,*)'  TO  ENTER  ANY  PRESSURE  THAT  YOU  WANT  TO  USE.' 

WRITE!*,*) 

WRITE!*,*)'    DO  YOU  WANT  TO  USE  SOME  OTHER  PRESSURE  BESIDES' 

WRITE!*,' (A\)')'  THE  150  PSIG  NORMALLY  USED?!Y>' 

READ!*, 310)  ST 

IF  !ST,EQ,'Y'.OR.ST.EQ.'Y\OR.ST.EQ.'  ')  THEN 

WRITE!*,*)  CLR,POS2 

WRITE!*,*)'    ENTER  THE  STATIC  PRESSURE  IN  PSIG.' 

READ!*,*)  P0 

WRITE!*,*) 

ELSE 

P0=150. 

ENDIF 

IF  (NB.EQ.0)  THEN 

GO  TO  556 

ELSE 

CONTINUE 

ENDIF 

WRITE!*,*)  CLR-.P0S2 

WRITE (*,'(A\)')'  DO  YOU  WANT  TO  CHECK  BOLT  THREAD  ENGAGEMENT? <Y>' 

READ(*,310)ST 

IF'ST.EQ.'Y'.OR.  ST.EQ.'Y'.OR.ST.EQ,'  ')  THEN 

NT1=1 

WRITE!*,*)  CLR,TOP 

WRITE!*,*)'  THREAD  ENGAGEMENT' 

WRITE!*,*) 

WRITE!*,*)'     THE  MINIMUM  LENGTH  OF  THREAD  ENGAGEMENT,' 

WRITE!*,*) 'REQUIRED  TO  DEVELOP  THE  FULL  STRENGTH  OF  THE  BOLTS,' 

WRITE!*,*) 'IS  CALCULATED  IN  THIS  SECTION.' 

WRITE!*,*) 

WRITE!*,*) 

CALL  RETURN (CLR) 

IG0=13 

NUMB=20 

WRITE!*,*) CLR, P0S2 

CALL  DISPLAY !IGO, NUMB, CLR) 

WRITE!*,*)  CLR, TOP 

WRFE!*,*)'  ENTER  THE  BOLT  TENSILE  STRESS  AREA  IN  SQUARE  INCHES.' 

READ(*,*)AT 

WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  MAXIMUM  MINOR  DIAMETER  OF  THE  INTERNAL  ' 

WRITE!*,*)'  THREAD  IN  INCHES.' 

READ  (*,♦)  BKNMAX 

WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  MINIMUM  PITCH  DIAMETER  OF  THE  EXTERNAL' 

WRITE!*,*)'  THREAD  IN  INCHES.' 


4S 


READ(*,*)ESMIN 

WRITE(*, 

WRITE!*, 

READ!*,* 

WRITE(*. 

WRITE!*, 

WRITE!*, 

WRITE!*, 

WRITE(*t 

WRITER, 

WRITE(*, 

WRITEt*, 

WRITE!*, 


)'  ENTER  THE  NUMBER  OF  THREADS  PER  INCH.' 
BT 
) 

)CLR,POS2 

)'    THE  MINIMUM  LENGTH  OF  THREAD  ENGAGEMENT  DEPENDS' 
)'  ON  WHETHER  THE  TAPPED  MATERIAL  (THE  FLANGE)  IS  AS' 
)'  STRONG  AS  OR  WEAKER  THAN  THE  BOLT.  THE  TENSILE' 


'  STRENGTH  SOULD  BE  USED  FOR  COMPARISON. ' 


)'     IS  THE  TAPPED  MATERIAL  AS  STRONG  AS  THE  BOLT-" 
WRITE!*,' (A\)')'  PLEASE  ANSWER  YES  OR  NO.W 
READ  <*,310)ST 

TLE=2. * ( AT) / (3. 14*BKNMAX*  i . 5+. 57735*BT* (ESMIN-BKNMA)! ) ) ) 
WRITE!*, *>CLR,P0S2 

IFtST.EQ. 'Y' .OR.ST.EQ. ' V .OR.ST.EQ- '  ')  THEN 
CONTINUE 
ELSE 

WRITE!*,*)  CLR,P0S2 

WRITE!*,*)'  ENTER  THE  TENSILE  STRENGTH  OF  THE  BOLT  MATERIAL' 
READ!*,*)TSB 
WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  TENSILE  STRENGTH  OF  THE  TAPPED  MATERIAL' 
READ!*,*)TSM 
WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  MINIMUM  MAJOR  DIAMETRER  OF  THE' 
WRITE (*,*)'  EXTERNAL  THREAD  IN  INCHES.' 
READ!*,*)  DSMIN 
WRITE!*,*) 

WRITE!*,*)  '  ENTER  THE  MAXIMUM  PITCH  DIAMETER  OF  THE' 
WRITE!*,*)  '  INTERNAL  THREAD  IN  INCHES.' 
READ!*,*)  ENMAX 
WRITE!*,*) 

AS=3. 14*TLE*BKNMAX*  < ( 1 . / (2. *BT) ) +. 57735* (ESMIN-BKNMAX ) ) 
AN=3.14*TLE*DSMIN  *((1./(2.*BT))+.57735*(DSMIN-ENMAX  )) 
TJ=(AS*TSB)/(AN*TSM) 
TLE=TLE*TJ 
ENDIF 

WRITE!*,*)  CLR,P0S2 

WRITE!*,*)'  THE  MINIMUM  THREAD  ENGAGEMENT  LENGTH  IS:' 
WRITE!*, *)TLE 
WRITE!*,*) 
WRITE!*,*) 
CALL  RETURN (CLR) 
ELSE 


CONTINUE 

ENDIF 

CONTINUE 

300 

FORMAT  ('  INPUT  FILE  NAME-'\) 

310 

FORMAT (A) 

556 

CONTINUE 

WRITE!*,*) CLR, P0S2 

WRITE!*, *) '  WORKING 
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H2  =  H  *  H 

H3  =  H  *  H2 

PI  =  3.141593 

PI2  =  PI  *  PI 

PI4  =  PI2  *  PI2 

PI6  =  PI2  *  PI4 

D  =  (EP  *  H3)/(12 

0*(1 

.0-MU*MU)) 

IF  (NB.EQ.0)  THEN 

GO  TO  555 

ELBE 

CONTINUE 

ENDIF 

Kll  =  EB  *  IB  /  (, 

>.0  * 

LB) 

K12  =  12.0  *  EB  * 

IB  * 

KS 

+ 

4.0 

*  AB 

*  GB 

*  LB 

K13  =  12.0  *  EB  * 

IB  * 

KS 

+ 

AB  * 

GB  * 

LB  * 

LB 

Kl  =  Kll  *  K12  /  K13 

K2  =  AB  *  EB  /  (2 

0  * 

LB) 

c 

c 

CALCULATE  AA  MATRIX 

c 

555 

CONTINUE 

*  LB 


DO  1010  M=1,NN 

DO  1000  N=1,NN 

EPSILON(M,N)  =  0.0 

AA(M,N)  =  0.0 
1000  CONTINUE 
1010  CONTINUE 

DO  1030  H=1,NN,2 

DO  1020  N=1,NN,2 

EPSILON(M,N)  =1.0 
1020  CONTINUE 
1030  CONTINUE 

II  =  0 

DO  1500  MP=1,NN 
DO  1400  NQ=1,NN 
II  =  II  +  1 
JJ  =  0 

AP  =  MP*PI/A 
BQ  =  NQ+PI/B 
DO  1300  M=1,NN 
DO  1200  N=1,NN 
JJ  =  JJ  +  1 

IF  (M  .EQ.  MP  .AND.  N  .EG.  NQ)  THEN 

C(IIiJJ)  =  <A*B*D*PI4/4.0)  *  ((MP*MP)/(A*A)+(NQ*NQ)/(B*B))**I 
ELSE 

C(II,JJ)  =  0 
END  IF 
AM  =  M*PI/A 
BN  =  N*PI/B 
DO  1100  1=1 jNB 

CI  =  AM  *  COS(AH»XB(D)  *  SIN(BN*YB(D) 
C2  =  AP  *  COS(AP*XB(D)  *  SIN(BQ*YB(D) 
C3  =  BN  ♦  SIN(AM*XB(D)  *  COS(BN*YB(D) 
C4  =  BQ  *  SIN(AP*XB(D)  *  COS(BQ*YB(D) 
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C5  =  SIN(AM*XB(I))  ♦  SIN(BN*YB(D) 
C6  =  SIN(AP*XB<D)  *  SIN(BQ+YB(D) 

C(ILJJ)  =  C(ILJJ)  +  2.0*K1  *  (C1*C2  +  C3*C4)  +  2.0*K2*C5*C6 
1100   CONTINUE 
1200  CONTINUE 
1300  CONTINUE 

K(II)  =  4.0  *  P0  *  A  ♦  B  *  EPSILON(MP,NQ)  /  (PI2  *  MP  *  NQ) 
1400  CONTINUE 
1500  CONTINUE 
C     BO  6665  M1=1,NN*NN 
C     WRITE  (2,6666)  (C(M1,N1) ,N1=1,NN*NN) 
C6665  CONTINUE 
C     WRITE  (2,6667)  <K(M2),M2=1,NN*NN) 

6666  FORMAT  (9E10.4) 

6667  FORMAT  (9E10.4) 

WRITE  (*,*)  'GOING  INTO  SIMULTANEOUS  SOLUTION' 
C     GO  TO  6668 

CALL  SIMUL  (C,K) 
WRITE  (*,*)'  SOLUTION  OBTAINED' 
1=0 

BO  1700  M=1,NN 
DO  1650  N=1,NN 
1=1+1 

AA(M,N)  =  K(I) 
1650  CONTINUE 
1700  CONTINUE 
C     GO  TO  6668 
C 

C   LOCATION  FOR  STRESS  OUTPUT 
C 

NS  =  3 
X(l)  =  A/2.0 
Yd)  =  B/2.0 
X(2)  =  A/2.0 
Y(2)  =  0.0 
X(3)  =  0.0 
Y(3)  =  B/2.0 
C 

C   COMPUTE  PLATE  STRESSES 
C 

DO  5000  1=1, NS 
W(I)  =  0.0 
WX(I)  =  0.0 
WXX(I)  =  0.0 
WXXX(I)  =  0.0 
WYU)  =  0.0 
WYY(I)  =  0.0 
WYYY(I)  =  0.0 
WXY!I)  =  0.0 
WXXY(I)  =  0.0 
WYYX(I)  =0.0 
DO  4950  M=1,NN 
AM  =  M  *  PI  /  A 
AM2  =  AM  *  AM 
AM3  =  AM  ♦  AM2 
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DO  4900  N=i,NN 

BN  =  N  *  PI  /  B 

BN2  =  BN  *  BN 

BN3  =  BN  *  BN2 

BA  =  SIN<AM*X(I)) 

SB  =  SINiBN*Y(D) 

CA  =  COS(AH*X(D) 

CB  =  COS(BN*Y'I)> 

C     PLATE  DEFLECTIONS 

W(I)  =  WH)  +  AA(M,N)  *  SA  * 

SB 

Mil)   =  WX  ( I )  + 

AA(M,N)  *  A!" 

*  CA 

*   S3 

WXX(I)  =  WXX(I)  - 

AA(M,N)  * 

AM2  * 

SA  * 

SB 

WXXX(I)  =  WXXX(I) 

-  AA(M,N) 

*  AM3 

*  CA 

*  SB 

WY(I)  =  WY(I)  + 

AA(M,N)  *  Bh 

*  SA 

*  CB 

WYY(I)  =  WYY(I)  - 

AA(M,N)  * 

BN2  * 

SA  * 

SB 

WYYY(I)  =  WYYYU) 

-  AA(M,N) 

*  BN3 

*  SA 

*  CB 

WXY«I)  =  WXY(I)  4 

AA(M,N)  * 

AM  *  BN  *  CA  *  CB 

WXXYU)  =  WXXY(I) 

-  AA(M,N) 

*  AM2 

*  BN 

*  SA  * 

CB 

WYYX(I)  =  WYYX(I) 

-  AA(M,N) 

♦  AM  i 

»  BN2 

*  CA  * 

SB 

4900   CONTINUE 

4950  CONTINUE 

MX  =  -D  *  (WXX(I)  + 

MU  *  WYY(D) 

MY  =  -D  *  (WYY(I)  + 

MU  *  WXX(D) 

MXY  =  D  *  (1  -  MU)  * 

WXYU) 

QX  =  -D  *  (WXXX(I)  + 

WYYX(D) 

OY  =  -D  *  (WXXY(I)  + 

WYYY(D) 

SIGMAX(I)  =  6.0  *  MX 

/  H2 

SIGMAYU5  =  6,0  *  MY  /  H2 

TAUXY(I)  =  6.0  *  MXY 

/  H2 

TAUXZ(I)  =  1.5  *  QX 

/  H 

TAUYZ(I)  =  1.5  *  QY 

/  H 

5000  CONTINUE 

C   COMPUTE  BOLT  STRESSES 
C 

DO  5300  1=1, NB 

WB'I)  =  0.0 
THETABX(I)  =0.0 
THETACY  (I)  =0.0 
DO  5200  M=1,NN 
DO  5100  N=1,NN 

AM  =  M  *  PI  /  A 

BN  =  N  *  PI  /  B 
]    BOLT  DEFLECTIONS 

WB(I)  =  NB1I)  +  AA(i1,N)  *  SIN(AM*XB(D)  *  SIN(BN*YBU)) 

THETABX(I)  =  THETABX(I)  +  AM  *  AA(M,N)  *  COS(AM*XB(D)  * 

*  SIN(BN»YB(D) 

THETACY(I)  =  THETACY(I)  +  BN  *  AA(M.N)  *  SIN(AM*XB<D)  * 

*  COS(BN*YB(D) 
5100  CONTINUE 

5200  CONTINUE 

SIGMAAX(I)  =  EB  *  WBU)  /  LB 

XM1  =  2.0  *  EB  *  IB  /  LB 

XM2  =  6.0  *  EB  *  IB  *  KS  +  2.0  •  AB  *  GB  *  LB  *  LB 

XM3  =  12.0  *  EB  *  IB  *  XS  +  AB  *  GB  *  LB  *  LB 
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MYT  =  XM1  ♦  XM2  ♦  THETABX(I)  /  XM3 

MXT  =  -XM1  *  XM2  ♦  IHETACY(I)  /  XM3 
XVI  =  XM1  /  LB 

XV2  =  3.0  *  AB  *  GB  *  LB  *  LB 
XV3  =  XM3 

VXT  =  -XVI  *  XV2  * 

VYT  =  -XVI  *  XV2  ♦ 
SIGMABX(I)  =  32.0  * 


SIGMACY(I)  =  32.0  + 


THETABX(I)  /  XV3 
THETACY(I)  /  XV3 
MYT  /  (PI 
MXT  /  (PI 


TAUBXZ(I) 
TAUBYZ(I) 
5300  CONTINUE 


16.0 
16.0 


VXT  /  (3.0 
VYT  /  (3.0 


DB*DB*DB) 

DB*BB*DB) 


BB*DB) 
BB*BB) 


1  A:  ',A,' 
'  THICKNESS: 
1)  THEN 


C  OUTPUT 
6668  CONTINUE 
WRITE (23, 
WRITE (23, 
WRITE (23, 
WRITE (23, 
WRITE (23, 
WRITE (23, 
WRITE (23, 
WRITE(23, 
WRITE(23, 
WRITE(23, 
IF  (NT1.EQ 
WRITE(23, 
WRITE(23, 
WRITE (23, 
WRITE(23, 
WRITE(23, 

IF  (TLE.LT.LB)  THEN 

WRITE (23 

ELSE 

WRITE (23 

WRITE (23 

WRITE (23 

WRITE(23 

WRITE (23 

WRITE (23 

WRITE(23 

WRITE (23 

WRITE (23 

WRITE(23 

WRITE (23 

WRITE (23 

WRITE (23 

ENDIF 

ELSE 

WRITE (23 

WRITE (23 

ENDIF 

WRITE  (2,*) 
WRITE(2,315)ID 
WRITE  (2,*) 


RESULTS  OF  COVER  SCREENING' 


THE  COVER  DIMENSIONS  ARE: 


B:  ',B 


',H 


THE  MINIMUM  CALCULATED  BOLT  THREAD  ENGAGEMENT  IS: ' 


H_E 


*)'  THE  BOLT  THREAD  ENGAGEMENT  IS  CORRECT' 

»)'  THE  BOLT  THREAD  ENGAGEMENT  IS  NOT  CORRECT' 

*) 

SOME  ALLOWANCE  CAN  BE  MADE  IF  AXIAL  BOLT  ' 
STRESS  ARE  LOW.  ADJUST  THE  MINIMUM  CALCULATED' 
BOLT  THREAD  ENGAGEMENT  BY  THE  FOLLOWING  EQUATION: ' 


*)' 

*)' 

*)' 

*) 

*)' 

*)' 

*)' 

♦  )' 

*)' 

*)' 

») 


*»' 

+  ) 


L(  ADJUSTED  )  =  (L  OLD)  *  1.5  *  K' 
WHERE  K  IS  : ' 

CALCULATED  BOLT  AXIAL  STRESS  FROM  ' 

ALGORITHM  AT  150  PSIG. ' 
K  =  ' 

BOLTED  TENSILE  STRENGTH' 


YOU  CHOSE  NOT  TO  CHECK  MINIMUM  THREAD  ENGAGEMENT. 


ENCLOSURE  IDENTIFICATION: 


53 


WRITE  (23,*) 

WRITE  (23,*) 

WRITE  (23,9020)  NN 

WRITE  (23,9030)  A,B,H 

WRITE  (23,9040)  P0 

WRITE  (23,*) 

WRITE  (23,*)  'PLATE  PROPERTIES' 

WRITE  (23,9050)  EJMIU.YP 

WRITE  (23,*) 

IF  (NB.EQ.0)  THEN 

WRITE (23,*) 'THE  COVER  DOES  NOT  HAVE  ANY  BOLTS,' 

GO  TO  557 

ELSE 

CONTINUE 

END  IF 

WRITE  (23,*)  'BOLT  PROPERTIES' 

WRITE  (23,9052)  EB,IB,GB,KS 

WRITE  (23,9054)  AB,LB,DB 

WRITE  (23,*) 
C     GO  TO  6669 
C     WRITE  (2,*)  'AA  -  MATRIX' 
C     DO  7100  11=1,  NN 
C     WRITE  (2,*)  M 
C     WRITE  (2,9060)  (AA(M,N) ,N=1,NN) 

WRITE  (23,*) 
7100  CONTINUE 
C     GO  TO  6669 

557  CONTINUE 

WRITE  (23,*)  'PLATE  DEFLECTIONS  AND  STRESSES' 
DO  7200  1=1, NS 
WRITE  (23,*) 

WRITE  (23,9070)  X(I),Y(I) 
WRITE  (23,9075)  W(I) 
WRITE  (23,9080)  SISMAX (I) ,SIGMAY(I) 
WRITE  (23,9090)  TAUXY(I) ,TAUXZ(I),TAUYZ(I) 
7200  CONTINUE 
WRITE  (23,*) 
WRITE  (23,*) 
IF  (NB.EQ.0)  THEN 
GO  TO  558 
ENDIF 

WRITE  (23,*)  'BOLT  STRESSES  AND  DEFLECTIONS' 
DO  7300  1=1, NB 

WRITE  (23,*) 

WRITE  (23,9100)  I,XB(I),YB(I) 

WRITE  (23,9075)  WB(I) 

WRITE  (23,9110)  SIGMAAX(I) 

WRITE  (23,9120)  SIGMABX(I),SIGMACY(I) 

WRITE  (23,9130)  TAUBXZ ( I ) , TAUBYZ (I) 
7300  CONTINUE 
6669  CONTINUE 

558  CONTINUE 
C   FORMATS 

313    FORMATUX,'     THIS  OPTION  IS  NOT  FULLY  SUPPORTED-TRY  AGAIN' 
9010  FORMAT  (60A1) 
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9015  FORMAT  (1X.64A1) 

9020  FORMAT  ('  DIMENSION  OF  AA  -  MATRIX  -  NN  =  ',14) 

9030  FORMAT ('  PLATE  DIMENSIONS  -  A  =  '.F9.4, 

*  ',  B  =  ',F9.4,',  H  =  '.F8.4) 
9040  FORMAT  ('  INTERNAL  PRESSURE  -  P0  =  '-F12.4) 

9050  FORMAT  ('  E  =  '.F12.2,',  MU  =  ',F6.4,',  YP=  \F12.2) 

9052  FORMAT  ('  E  =  '.F12.2,',  I  =  SF12.6,',  6  =  *,F12.2»*i  KS  =  ', 

*  F8.4) 

9054  FORMAT  ('  A  =  ',F12.6,',  L  =  ',F8.4,',  D  =  ',F8.4) 

9060  FORMAT  (5E14.4) 

9070  FORMAT  ('  X  =  ',F9.4,',  Y  =  ',F9.4) 

9075  FORMAT  ('  DEFLECTION  =  ',F8.5) 

9080  FORMAT  ('  SIGMA-X  =  ',F12.4,',  SIGMA-Y  =  '.F12.4) 

9090  FORMAT  ('  TAU-XY  =  '.F12.4,1,  TAU-XZ  =  ',F12.4, 

*  ',  TAU-YZ  =  ',F12.4) 

9100  FORMAT  ('  POSITION  OF  BOLT  -  ',12,' ,  X  =  ',F9.4, ' ,  Y  =  SF9.4) 

9110  FORMAT  ('  AXIAL  STRESS  =  ',F12.4> 

9120  FORMAT  ('  BENDING  STRESSES  -  SIGMA-X  =  ',F12.4,',  SIGMA-Y  -  ', 

*  F12.4) 

9130  FORMAT  ('  TAU-XZ  =  '.F12.4,',  TAU-YZ  =  '.F12.4J 
WRITE(*,*)  '  OUTPUT  FILE  IS  PRINTED' 
CALL  RETURN (CLR) 
160=14 
NUMB=20 

CALL  DISPLAY (ISO, NUMB, CLR) 
IF  (NB.EQ.0)  THEN 
CONTINUE 
ELSE 

WRITE(*,*)'  YOU  MUST  ENTER  THE  YIELD  STRENGTH  OF  THE  BOLT' 
WRITE(*,*)'  IN  KSI.  FOR  EXAMPLE  THE  YIELD  STRENGTH  OF  A-36' 
WRITE(*,*)'  STEEL  WOULD  BE  ENTERED  AS  36.' 
READ  (*,*)  BYS 
BYS=BYS*1000. 
WRITE(*,*)  CLR.P0S2 

WRITEf*,*)'  YOU  MUST  ENTER  THE  BOLT  ULTIMATE  STRENGTH  IN  KSI.' 
READ  (*,*)  BUS 
BUS=BUS*1000. 
WRITEf*,*)  CLR,P0S2 
IF11=0 
WRITE (23,*) 
WRITE(23,*) 

WRITE(23,*)'  THE  YIELD  STRENGTH  OF  THE  BOLT  IS:  \6YS 
WRITE(23,*) 

WRITE(23,*)'  THE  BOLT  ULTIMATE  STRENGTH  IS:  '.BUS 
WRITE(23,*) 
DO  4000  1=1, NB 

IF<  SIGMAAX(I).6T.BYS)  THEN 

WRITE(23,*)'  PREDICTED  BOLT  AXIAL  STRESS  IS  GREATER  THAN' 

WRITE(23,*)'  BOLT  YIELD  STRENGTH  AT  BOLT  NUMBER  ',1 

WRITE(23,») 

WRITE(23,*)'  THE  COVER  SHOIXD  BE  REJECTED.' 

WRITE(23,») 

WRITE(*,*)CLR,P0S2 

WRITE(*,*)'  PREDICTED  BOLT  AXIAL  STRESS  IS  GREATER  THAN' 

WRITE(*,*)'  BOLT  YIELD  STRENGTH  AT  BOLT  NUMBER  ',1 
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WRITE!*,*) 

WRITE!*. *>'  THE  COVER  SHOULD  BE  REJECTED.' 

WRITE!*,*) 

CALL  RETURN (CLR) 

ELSE 

IF11=1+IF11 
ENDIF 

4000  CONTINUE 

IF  (IF11.EQ.NB)  THEN 

WRITE(23,*)'  ALL  OF  THE  BOLT  AXIAL  STRESSES  ARE  LESS  THAN  THE' 
WRITE(23,*)'  BOLT  YIELD  STRENGTHS.' 
WRITE!*,*) CLR,  F0S2 

WRITE!*,*)'  ALL  OF  THE  BOLT  AXIAL  STRESSES  ARE  LESS  THAN  THE' 
WRITE(*,*)'  BOLT  YIELD  STRENGTHS.' 
WRITE!*,*) 
CALL  RETURN (CLR) 

XBX=SIGMAAX ( I ) +S IQMABX ( I ) +SIGMACY ( I ) 
DO  4011  1=1, NB 
IF(XBX.GT.BUS)  THEN 

WRITE!23,*)'  PREDICTED  COMBINATION  OF  BENDING  AND  AXIAL  STRESS' 
WRITEC23,*)'  EXCEEDS  THE  BOLT  ULTIMATE  STRENGTH  AT  BOLT.,....' 
WRITE(23,*)I 
WRITE<23,*) 

WRITE (23,*)'  THE  COVER  SHOULD  BE  REJECTED.' 
WRITE(23,*) 
WRITEC*,*)  CLR,P0S2 
WRITE!*,*)'  PREDICTED  COMBINATION  OF  BENDING  AND  AXIAL  STRESS' 

WRITE(*,*)'  EXCEEDS  THE  BOLT  ULTIMATE  STRENGTH  AT  BOLT ' 

WRITE(*,*)I 

WRITE(*,*) 

WRITE!*,*)'  THE  COVER  SHOULD  BE  REJECTED.' 

WRITE!*,*) 

CALL  RETURN (CLR) 

ELSE 

CONTINUE 

ENDIF 

4001  CONTINUE 

ELSE 

ENDIF 

ENDIF 

DO  4003  1=1,  NS 

SS=SIGMAX(I)+SIGMAY(I) 

IF(SS.GT.YP)  THEN 

WRITE(23,*) 

WRITE(23,*)'  PREDICTED  COVER  STRESSES  EXCEED  THE  COVER' 

WRITE(23,*)'  YIELD  STRENGTH.  THE  COVER  SHOULD  BE  ANALYZED' 

WRITE!23,*)'  IN  THE  NEXT  SECTION  AS  A  PANEL.' 

WRITE!*,*) 

WRITE!*,*)'  PREDICTED  COVER  STRESSES  EXCEED  THE  COVER' 

WRITE!*,*)'  YIELD  STRENGTH.  THE  COVER  SHOULD  BE  ANALYZED' 

WRITE!*,*)'  IN  THE  NEXT  SECTION  AS  A  PANEL.' 

WRITE!*,*) 

CALL  RETURN (CLR) 

ELSE 

ENDIF 
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4003  CONTINUE 
RETURN 
END 

C  ********  ********************************** ♦ ************** 

c  *  * 

C  *  SUBROUTINE  DISPLAY        (DISPLAY. FOR)  * 

C  *  * 

c  ***********♦********+************************  ************ 

SUBROUTINE  DISPLAY ( I GO, NUMB. CLR) 

CHARACTER*!  TEXT (80), CLR (4) 

IF  (IGO.EQ.l)  THEN 

OPEN  (21, FILE=' FOOTl.DOC') 

ELSEIF  (IGO.EQ. 2)  THEN 

OPEN  (21,  FILE=' DISP3.DOC) 

ELSEIF  (IG0.EQ.3)  THEN 

OPEN  (21,  FILE=' DISP4.DOC) 

ELSEIF  (IGO.EQ.4)  THEN 

OPEN  (21,  FILE=' DISP5.DOC) 

ELSEIF  (IG0.EQ.5)  THEN 

OPEN  (21,  FILE=' DISP6.DOC) 

ELSEIF  (IG0.EQ.6)  THEN 

OPEN  (21,  FILE=' DISP7.DOC) 

ELSEIF  (IGO.EQ. 7)  THEN 

OPEN  (21, FILE=' DISP8.DOC') 

ELSEIF  (IG0.EQ.8)  THEN 

OPEN  (21,  FILE=' DISP9.DOC) 

ELSEIF  (I60.EQ.9)  THEN 

OPEN  (21,  FILE=' DISF10.DOC) 

ELSEIF (IGO.EQ. 10)  THEN 

0PEN(21,FILE='DISP11.D0C') 

ELSEIF(IGO.EQ.ll)  THEN 

0PEN(21,FILE='DISP12.D0C') 

ELSEIF (IGO.EQ. 12)  THEN 

0PEN(21,FILE='DISP13.D0C) 

ELSEIF (IGO.EQ. 13)  THEN 

OPEN  (21,  FILE=' DISP14.DOC) 

ELSEIF (IGO.EQ. 14)  THEN 

0PEN(21,FILE='DISP15.D0C) 

ELSEIF (IGO.EQ. 15)  THEN 

OPEN  (21,  FILE='  DISP16.DOC) 

ELSEIF (IGO.EQ. 16)  THEN 

OPEN  (21,  FILE=' DISP17.DOC) 

ELSEIF (IGO.EQ. 17)  THEN 

0PEN(21,FILE='DISP18.D0C') 

ELSEIF (IGO.EQ. 18)  THEN 

0PEN(21,FILE='DISP19.D0C) 

ELSEIF (IGO.EQ. 20)  THEN 

OPEN (21, FILE=' DISP1.DOC') 

ELSEIF (IGO.EQ. 21)  THEN 

OPEN  (21,  FILE=' DISP2.DOC) 

ELSEIF (IGO.EQ. 22)  THEN 

OPEN (21, FILE=' DISP20.DOC') 

ELSE 

CONTINUE 

ENDIF 
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WRITE(*,*)CLR 

DO  3100  IX=1,NUMB 

READ  (21,3101) (TEXT (JX),JX=1, 80) 

WRITE(*,3102)(TEXT(JX),JX=1,79) 

3100  CONTINUE 

3101  FORMAT (80A) 

3102  FORMAT (1X.79A1) 
WRITE;*,*)'    » 
WRITE(*,*)'    ' 
REWIND  (21) 
CLOSE  (21) 

CALL  RETURN (CLR) 

RETURN 

END 

c    *  * 

C     *         SUBROUTINE  FIS4   (FI64.F0R)  * 

C     *  * 

C    $**M*t***t*****iM*m*m****«*m**m**********t$m**** 
SUBROUTINE  FI64(X,V) 

WRITE (*,*)'  RATIO  OF  TR/T  USED  TO  OBTAIN  DLF',X 
READ  (*,*)  X 
IF  (X.LE.1.0)  THEN 
V=1.0*  C0S(1.57*X) 
V=V+1. 
ELSEIF  ((X.LE.2.0KAND. (X.6T.1.0))  THEN 
XX=X-1.0 

V=.22*SIN(3.14*XX) 
V=V+1. 
ELSEIF  ((X.LE.3.0).AND.(X.GT.2.0))  THEN 
XX=X-2.0 

V=.14*SIN(3.14*XX) 
V=V+1, 
ELSEIF  ((X.LE.4.0J.AND. (X.GT.3.0))  THEN 
XX=X-3.0 

V=.1*SIN(3.14*XX) 
V=V+1. 
ELSE 

V=1.0 
ENDIF 

WRITE(*,*>  'VALUE  OF  DLF  RETURNED',  V 
RETURN 
END 

C   »»**»***♦**»♦»********♦♦*********♦*************»********+*** 

c    *  * 

C    *         SUBROUTINE  FLAMEP   (FLAMEP.FOR)  * 

C    *  * 

SUBROUTINE  FLAMEP (CLR, I FLAG) 
CHARACTER*64  FLM2(33) 
CHARACTER* 1  CLR (4) 
OPEN  (10,FILE='FLAMEPR,DAT') 
OPEN  dl,FILE=' FLAME1.DAT') 
WRITE  (*,*)  CLR 
ILOOP=0 
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READ  (10,'(A64)')FLM2(1) 
WRITE  (4,'(1X,A64)')FLM2(1) 
WRITE  (4,4) 
WRITE  (4,4) 
DO  2100  1=1,5 

2100  CALL  PRINT1 ( I , ILOOP. IFLAG) 
2104   WRITE  (*,*) 

WRITE  (4,4) 

WRITE  (4,4)'  MORE  DATA  ON  NEXT  SCREEN ' 

WRITE  (4,4) 
1350   WRITE  (*,'(A\)')  '  TO  VIEW  NEXT  SCREEN,  TYPE  RETURN' 

READ  (*,310)  ST 
310   FORMAT (A) 

IF  (ST.EQ.'  ')  THEN 

WRITE  (4,4) 

ELSE 

WRITE  (4,4)  '  YOU  MUST  TYPE  RETURN  TO  CONTINUE' 

GO  TO  1350 

ENDIF 

WRITE  (4,4)  CLR 

DO  2101  1=1,6 

2101  CALL  PRINT1 ( I , ILOOP, IFLAG) 
WRITE  (4,4) 

WRITE  (4,4) 
2103   WRITE  (4,4)  '  IF  YOU  WANT  TO  REVIEW  PREVIOUS  SCREEN,  TYPE  1' 
WRITE(4,4) 

WRITE  (4,MA\)')  '  OTHERWISE,  TYPE  RETURN' 
READ  (4,310)  ST 
IF  (ST.EQ. '1')  THEN 
REWIND  10 
REWIND  11 
ILOOP=0 

WRITE  (4,4)  CLR 
READ  (10,'(A64)'5FLM2(1) 
WRITE  (4,'(1X,A64)MFLM2(1) 
WRITE  (4,4) 
WRITE  (*,*) 
DO  2102  1=1,5 

2102  CALL  PRINT1 ( I , ILOOP, IFLAG) 
GO  TO  2104 

ELSEIF  (ST.EQ.'  ')  THEN 

REWIND  10 

REWIND  11 

CLOSE (10) 

CLOSE  (11) 

WRITE (*,4)  CLR 

RETURN 

ELSE 

WRITE(4,4)  '  YOU  CAN  ONLY  TYPE  A  1  OR  A  RETURN' 

GO  TO  2103 

ENDIF 

END 
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C     *  * 

C     *  SUBROUTINE  FOOT   (FOOT, FOR)  * 

C     *  » 

SUBROUTINE  FOOT (160, NUMB, CLR) 
CHARACTERS  TEXT (80), CLR (4) 
IF  (IGO.EQ.l)  THEN 
OPEN  (21,  FILE=' F00T1.DOC) 
ELSEIF  (I60.EQ.2)  THEN 
0PEN(21,FILE='TW0.D0C) 
ELSEIF  (I80.EQ.3)  THEN 
OPEN  (21,  FILE=' THREE.  DOC) 
ELSEIF  (IGO.EQ. 4)  THEN 
OPEN  (21,  FILE='  FOUR.  DOC) 
ELSEIF  (IG0.EQ.5)  THEN 
OPEN  (21,  FILE='  FIVE,  DOC) 
ELSEIF  (IGO.EQ. 6)  THEN 
0PEN(21,FILE='SIX.D0C) 
ELSEIF  (I60.EQ.7)  THEN 
OPEN  (21, FILE=' SEVEN. DOC) 
ELSEIF  (ISO. EG. 8)  THEN 


OPEN  (21, FILE 


ELSEIF  UG0.EQ.9)  THEN 
OPEN  (21,FILE='NINE.D0C) 


'EI8HT.DOC') 


ELSEIF ( 180,  EQ. 
0PEN(21,FILE=' 
ELSEIF ( ISO. EQ. 
0PEN(21,FILE=' 
ELSEIF ( ISO. EQ. 
0PEN(21,FILE=1 
ELSEIF (IGO.EQ. 
QPEN(21,FILE=' 
ELSEIF (IGO.EQ. 
0PEN(21,FILE=' 
ELSEIF (IGO.EQ. 
0PEN(21,FILE=' 
ELSEIF ( ISO. EQ. 
0PEN(21,FILE=' 
ELSEIF (IGO.EQ. 
0PEN(21,FILE=' 
ELSEIF (IGO.EQ. 
0PEN(21,FILE=' 
ELSEIFdGO.EQ. 
0PEN(21,FILE=' 
ELSEIFdGO.EQ. 
OPEN  (21,  FILE---' 
ELSEIFdGO.EQ. 
0PEN(21,FILE=' 

else 
continue 

ENDIF 

WRITE(*,*)CLR 

DO  3100  IX=1,NUMB 

READ  (21,3101) (TEXT(JX),JX=1, 80) 


10)  THEN 
TEN.  DOC) 

11)  THEN 
ELEVEN.DOC) 

12)  THEN 
TWELVE.  DOC) 

13)  THEN 
THIRT.DOC) 

14)  THEN 
F0URT.DOC) 

15)  THEN 
FIFT.DOC) 

16)  THEN 
SIST.DOC) 

17)  THEN 
SEVENT.DOC) 

18)  THEN 
EIGHTT.DOC) 

20)  THEN 
TWENT.DOC) 

21)  THEN 
TWENT1.DOC) 

22)  THEN 
TWENT2.DOC) 
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WRITE(*,3102)(TEXT(JX),JX=1,79) 

3100  CONTINUE 

3101  FORMAT (80A) 

3102  F0RMAT(1X,79A1) 
WRITER.*)'    ' 
WRITE(*,*>'    ' 
REWIND  (21) 
CLOSE  (21) 

CALL  RETURN (CLR) 

RETURN 

END 

0        ****************  *******  *******************  **************** 
C        *  * 

C     *        SUBROUTINE  MAT    (MAT. FOR)  * 

C     *  * 

C     ********************************************************** 

SUBROUTINE  MAT(P0S2,CLR,T0P) 
CHARACTERS  POS2(10),CLR(4),TEXT(80),TOP(9) 

3101  FORMAT (80A) 

3102  FORMAT (1X,79A1) 
3120   CONTINUE 

3113   WRITE  (*,*)  CLR.P0S2 

WRITE  (*,*)'  DO  YOU  WANT  TO' 

WRITE  <*,*) 

WRITE  (*,*) 

WRITE  (*,*)'     (1)  CHECK  ACCEPTABLE  STEELS' 

WRITE  (*,*) 

WRITE  (*,*)'     (2)  CHECK  ACCEPTABLE  ALUMINUMS' 

WRITE  (*,*) 

WRITE(*,*)  '     (3)  CHECK  ACCEPTABLE  GLASSES  AND  PLASTICS' 

WRITE  (*,*) 

WRITE  (*,*)'    (4)  CHECK  ACCEPTABLE  ADHESIVES  ' 

WRITE  (*,*) 

WRITE  (*,♦>'     (5)  RETURN  TO  MAIN  MENUE' 

WRITE(*,*) 

WRITE(*,'(A\)')'  INPUT  YOUR  CHOICE—)' 

READ  (*,'(BN,I6)')  J2 

WRITE  (*,*)  CLR 

IFLAGM=1 

IF  (J2.EQ.1)  THEN 

OPEN ( 13, FILE=* MATSTEEL. DAT' ) 

OPEN ( 14, FILE=' STEELDES. DAT' ) 
DO  3111  JX=1.6 

READ  (13,3101) (TEXT(IX).IX=1, 80) 
WRITE  (*,3102)(TEXT(IX),IX=1,79> 

3111  CONTINUE 

DO  3112  JX=1,18 

READ  (14,3101) (TEXT(IXi,IX=l, 80) 

WRITE  (*,3102)(TEXT(IX),IX=1,79) 

3112  CONTINUE 
WRITE(*,*) 
REWIND  13 
REWIND  14 

CALL  RETURN (CLR) 

CALL  STEELR(CLR,P0S2,YP,EP,PR,DEN,IFLASM,T0P  ) 
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CLOSE (13) 
CLOSE (14) 
ELSEIF  (J2.EQ.2)  THEN 

OPEN ( 15, FILE*' MATAL.DAT') 

OPEN ( 16, FILE=' ALDES.DAT') 

DO  3114  JX=1,6 

READ  (15,3101) (TEXT(IX),IX=1, 80) 

WRITE  (*,3102)(TEXT(IX),IX=1,79) 

3114  CONTINUE 

DO  3115  JX=1,16 

READ  (16,3101) (TEXT(IX),IX=1, 30) 

WRITE  (*,3102)(TEXT(IX),IX=1,79) 

3115  CONTINUE 
WRITE(*,*) 
REWIND  15 
REWIND  16 

CALL  RETURN(CLR) 

CALL  flLUHR(CLR,P0S2iYPiEPiPR.DEN,  IFLASMJOP) 
CLOSE (15 J 
CLOSE (16) 
ELSEIF  (J2.EQ.3)  THEN 
OPEN ( 17, FILE=' MAT6LASS, DAT' ) 
OPEN(ia,FILE='GLASSDES.DAP) 
DO  3116  JX=1,6 

READ  (17,3101) (TEXT(IX),IX=1, 80) 
WRITE  (*,3102)(TEXT(IX),IX=1,79) 

3116  CONTINUE 

DO  3117  JX=1 , 12 

READ  tI8:318in(TEXT(IXJtIX=l,BB) 

WRITE  (*,3i02MT£XiUX.i,IX-i,7?l 

3117  CONTINUE 
WRITE(*,») 
REWIND  17 
REWIND  18 

CALL  RETURN (CLR) 
:        CALL  ALUP1R  (CLR, P0S2,  YP,  EP ,  PR,  DEN,  IFLAGM,  TOP) 
CLOSE (17) 
CLOSE (18) 
ELSEIF  (J2.EQ.4)  THEN 

0PEN(19,FILE='MATSEAL,DAT') 
OPEN (20, FILE=' SEALDES. DAT' ) 
DO  3118  JX=1,6 

READ  (19,3101) (TEXT(IX),IX=1, 80) 
WRITE  (*,3102)(TEXT(IX),IX=1,79) 

3118  CONTINUE 

DO  3119  JX=i,6 

READ  (20,3101)(TEXT(IX),IX=1,80) 

WRITE  (*,3102)(TEXT(IX),IX=1,79) 

3119  CONTINUE 
WRITE(*,*) 
REWIND  19 
REWIND  20 

CALL  RETURN (CLR) 
CLOSE (19) 
D_OSE!20) 
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ELSEIF  (J2.EQ.5)  THEN 
CALL  RETURN (CLR) 
RETURN 
ELSE 
WRITE(*,*)'  YOU  ONLY  HAVE  FIVE  CHOICES?' 

WRITE(*,*)'   TRY  AGAIN.., ' 

GO  TO  3120 
ENDIF 

GO  TO  3120 
END 
C    MM«m*****M^m#*«m*m***m********M**t**»*+*t** 

c        *  * 

C    *       SUBROUTINE  RETRI   (R.FOR)  * 

C    *  * 

c   *♦♦*♦♦♦***»*♦*♦*#*«**##****♦#♦♦*♦*»#******♦*♦**♦♦*♦***♦***** 

SUBROUTINE  RETRI (CLR,P0S2,YP,EP,PR,DEN, IFLAGM, IFLAM.TOP) 

CHARACTERS  CLR  (4) , P0S2  ( 10) , TOP (9) 

0PEN(14,FILE='STEELDES.DAT' ) 

OPEN ( 16, FILE=' ALDES.DAT') 

WRITE(*,*)  CLR,P0S2 

WRITE(*.*)'   WHAT  IS  THE  MATERIAL?  ' 

WRITE(*,*) 

WRITE**,*)'      (1)  STEEL' 

WRITE(*,*)'      (2)  ALUMINUM' 

WRITE(*,*) 

WRITE(*,'(A\)')'  INPUT  YOUR  CHOICE >' 

READ  (*,'(BN,I6)')IQ 

WRITE (*.*) CLR, P0S2 

IF  (IQ.EQ.l)  THEN 

IFLAGK=0 

CALL  STEELR (CLR, P0S2, YP, EP, PR, DEN, IFLAGM, TOP) 

ELSEIF  (IQ.EQ.2)  THEN 

CALL  ALUMR (CLR, P0S2, YP, EP, PR, DEN, IFLAGM, TOP) 
ENDIF 
CLOSE (14) 
CLOSE (16) 
RETURN 
END 

SUBROUTINE  STEELR ( CLR, P0S2,YP,EP, PR, DEN, IFLAGM, TOP) 
CHARACTERS  CLR(4) ,POS2(10) ,T0P(9) , AB 
INTEGER  IT(7),ISS(7) 
CHARACTER*1  TEXTX(18,80> .TEXT (25) 
WRITE(*,*)  CLR, TOP 

WRITE(*,*)'  THE  FOLLOWING  IS  A  LIST  OF  ACCEPTABLE  STEELS  THAT' 
WRITE(*,*)'  ARE  NORMALLY  USED  FOR  XP  ENCLOSURE  CONSTRUCTION.' 
WRITE(*,*) 
WRITE(*,*) 
DO  4110  JX=1,18 
4110  READ(14,3101)  (TEXTX(JX, IX), IX=1,80> 
DO  4111  JX=1,7 

ISS(JX)=1 

IF(JX.EG.l)  THEN 

IS=1 

ISS(JX)=0 

ELSEIF  (JX.EQ.2)  THEN 
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IS=3 

ISS(JX)=0 

ELSEIF  (JX.EQ.3)  THEN 
I S-5 

ELSEIF  (JX.EQ.4)  THEN 
IS=8 

ELSEIF  (JX.EQ.5)  THEN 
13=11 

ELSEIF  (JX.EQ.6)  THEN 
IS=14 
ELSE 
IS=17 
END1F 
IT(JX)=IS 
DO  4  IXX=1,80 
4     TEXTX(JX,IXX)=TEXTX(IS,IXX) 

4111  WRITE  (*,4112)'  (',JX,')',<TEXTX(JX,IY),IY=1,25) 

4112  FORMAT  (IX, 1A,I5,1A,3X,25A1) 
REWIND  14 

WRITE!*, *> 

WRITE!*,*) 

WRITE!*,*)'  IF  THE  MATERIAL  USED  IN  YOUR  ENCLOSURE  IS  NOT' 

WRITE!*,*)'  LISTED,  ENTER  AN  8.  OTHERWISE,  ENTER  THE  NUMBER' 

WRITE(*,*)'  OF  THE  MATERIAL.' 

WRITE(*,*) 

WRITE!*,*) '  ENTER  THE  NUMBER  FOR  THE  MATERIAL' 

WRITE(*,*J 

WRITE(*,'(A\)')  '  INPUT  YOUR  CHOICE >' 

READ  (*,'<BN,I6)')JX 

WRITE!*,*)  CLR,P0S2 

IF(IFLAGM.EQ.l)  THEN 

GO  TO  6 

ELSE 

CONTINUE 

ENDIF 

IF  (JX.EQ.8)  THEN 

WRITE!*,*)  CLR,  P0S2 

WRITE!*,*)'  YOU  MUST  ENTER  YP,EP,PR,  AND  DEN  FOR  YOUR  MATERIAL' 

WRITE!*,*) 

WRITE!*,*) 'ENTER  YIELD  STRENGTH  IN  KSI.  FOR  EXAMPLE,  YIELD  STRE 
1NGTH' 

WRITE!*,*) 'FOR  2024-T351  WOULD  BE  ENTERED  AS  42.' 

READ  (*,*)  YP 

WRITE!*,*) 'ENTER  ELASTIC  MODULUS  IN  KSI.  FOR  EXAMPLE, 
1ELASTIC 

WRITE (*,*)' MODULUS  FOR  2024-T351  WOULD  BE  ENTERED  AS  10730.' 

READ  !*,*)  EP 

WRITE!*, *i 'ENTER  POISSONS  RATIO  AS  A  DECIMAL' 

READ  (*,*)  PR 

WRITE!*, ♦) 'ENTER  DENSITY  IN  LB/IN3.  FOR  EXAMPLE,  DENSITY  FOR' 

WRITE!*, *)'2024-T351  WOULD  BE  ENTERED  AS  .100' 

READ  !*,*)  DEN 

ELSE 

ITT=IT(JX)-1 

DO  5  ITTT=1,ITT 
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5  R£AD(14,3102)AB 

READ  (14,3103)  YP,EP, PR-DEN 

ENDIF 

WRITE  (*!*) 

WRITE**, ♦>  CLR.P0S2 

WRITE(*,*)'         YIELD         ELASTIC      POISSONS 
1       DENSITY' 

WRITE(*,*)'        STRENGTH       MODULUS      RATIO 
1  ' 

WRITE(*,*)'         <KSI)         (KSI) 
1      (LB/IN3)' 

WRITE(*,*) 

WRITE  (*,*)YP,EP,PR-DEN 

3101  FORMAT (80A) 

3102  FORMAT (80A) 

3103  FORMAT(26X,F2.0,10X,F5.0,10X,F3.2,8X,F4.3> 
REWIND  14 

GO  TO  7 

6  CONTINUE 

IF  (JX.EQ,8)  THEN 

WRITE(23,*)'  THE  STEEL  USED  IN  THE  ENCLOSURE  IS  NOT  LISTED' 
WRITE(23,*)'  IN  THE  DATA  BASE  OF  ACCEPTABLE  CONSTRUCTION' 
WRITE(23,*)'  MATERIALS.' 
ELSE 

ITT=IT(JX)-1 
DO  8  ITTT=1,ITT 
8     READ  (14, 3102) AB 

WRITE(23,*>'  THE  STEEL  USED  IN  THE  ENCLOSURE  WAS  FOUND  IN' 
WRITE(23,*)'  THE  DATA  BASE  OF  ACCEPTABLE  CONSTRUCTION' 
WRITE(23,*)'  MATERIALS.  THE  STEEL  USED  IS:' 
IF(ISS(JX).EQ.0)  THEN 
READ(14,3i04)(TEXT(JXX),JXX=l,25) 
WRITE(23,*HTEXT(JXX),JXX=1,25) 
ELSE 

READ(14,3104)(TEXT(JXX),JXX=1,25) 
WRITE(23,*)(TEXT(JXX),JXX=1,25) 
READ(14,3104)(TEXT(JXX),JXX=1,25) 
WRITE(23,*)(TEXT(JXX),JXX=1,25) 
WRITE (23,*) 
ENDIF 
ENDIF 
REWIND  14 

7  CONTINUE 

3104  FORMAT (25A) 
RETURN 

END 
SUBROUTINE  ALUMR(CLR,P0S2,YP,EP,PR,DEN, IFLAGM,TOP) 
CHARACTERS  CLR(4) ,POS2(10),TOP(9) -AB 
INTEGER  IT(7),ISS(7) 
CHARACTER*1  TEXTX(16,80) .TEXT (25) 
WRITE)*, *)  CLR,TOP 

WRITE (*,*)'  THE  FOLLOWING  IS  A  LIST  OF  ACCEPTABLE  ALUMINUMS  THAT' 
WRITE(*,*)'  ARE  NORMALLY  USED  FOR  XP  ENCLOSURE  CONSTRUCTION.' 
WRITE(*,*) 
WRITE(*,*) 
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DO  4110  JX=i,16 

4110  READ (16, 3131)  (TEXTX(JX,IX),IX=1,80) 
DO  4111  JX=t,7 

ISS(JX)=1 
IF(JX.EQ.l)  THEN 
IS=1 

iss<JX)=a 

ELSEIF  (JX.EQ.2)  THEN 
IS=3 

ISS(JX)=0 

ELSEIF  (JX.EQ.3)  THEN 
IS=5 

ISS(JX)=0 

ELSEIF  (JX.EQ.4)  THEN 
IS=7 

ISS(JX)=0 

ELSEIF  (JX.EQ.5)  THEN 
IS=9 

ELSEIF  (JX.EQ.6)  THEN 
IS=12 
ELSE 
IS=15 
ENDIF 
IT(JX)=IS 
DO  4  IXX=1,80 
4     TEXTX(JX,IXX)=TEXTX(IS,IXX) 

4111  WRITE  (*,4112)'  (',JX,')',(TEXTX(JX,IY),IY=1,25» 

4112  FORMAT  (1X,1A,I5,1A,3X,25A1) 
REWIND  16 

WRITEt*,*) 

WRITEt*,*} 

WRITE(*,*J' 

WRITE(*,*)' 

WRITE**.*)' 

WRITEt*,*) 

WRITEt*,*)' 

WRITEt*,*) 

WRITEt*,' (A\)')  '  INPUT  YOUR  CHOICE >' 

READ  (*,'(BN,I6)')JX 

WRITEt*, *)CLR,P0S2 

IF(IFLASH.EQ.l)  THEN 

GO  TO  6 

ELSE 

CONTINUE 

ENDIF 

IF  (JX.EQ.8)  THEN 


IF  THE  MATERIAL  USED  IN  YOUR  ENCLOSURE  IS  NOT' 
LISTED,  ENTER  AN  8.  OTHERWISE,  ENTER  THE  NUMBER' 
OF  THE  MATERIAL. ' 

ENTER  THE  NUMBER  FOR  THE  MATERIAL' 


WRITE (* 
WRITE(* 
WRITEt* 
WRITEt* 
READ  (* 
WRITEt* 
READ  (* 
WRITEt* 


*)  CLR,  P0S2 

*)'  YOU  MUST  ENTER  YP,EP,PR,  AND  DEN  FOR  YOUR  MATERIAL' 

*) 

*) 'ENTER  YP  ' 

*)YP 

♦  )'  ENTER  EP  ' 

*)  EP 

*)'  ENTER  PR  ' 


READ'*,*)  PR 
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WRITER.*)'  ENTER  DEN  ' 

READ  (*,*)  DEN 

ELSE 

ITT=IT(JX)-1 

DO  5  ITTT=1, ITT 

5  READ (16, 3102) AB 

READ  (16,3103)  YF,EP,PR,DEN 

END  IF 

WRITE(*,*)  CLR,P0S2 

WRITE(*,*)'  YP  ','  EP  ',  '  PR  ','  DEN' 

WRITE  (*,*)YP,EP,PR,DEN 

3101  FORMAT (80A) 

3102  FORMAT (80A) 

3103  FORMAT(26X,F2.0,10X,F5.0,10X,F3.2,8X,F4.3) 
WRITE(*,*) 

WRITE(*,*) 
CALL  RETURN (CLR) 
REWIND  16 
60  TO  7 

6  CONTINUE 

IF  (JX.EQ.8)  THEN 

WRITE(23,*)'  THE  ALUMINUM  USED  IN  THE  ENCLOSURE  IS' 
WRITE(23,*)'  NOT  LISTED  IN  THE  DATA  BASE  OF  ACCEPTABLE' 
WRITE<23,*)'  CONSTRUCTION  MATERIALS.' 
ELSE 

ITT=IT(JX)-1 
DO  8  ITTT=1,ITT 
8     READ (16, 3102) AB 

WRITE(23,*)'  THE  ALUMINUM  USED  IN  THE  ENCLOSURE  WAS' 
KRITE(23,*)'  FOUND  IN  THE  DATA  BASE  OF  ACCEPTABLE  ' 
WRITE(23,*)'  CONSTRUCTION  MATERIALS.  THE  ALUMINUM' 
WRITE(23,*)'  IS:' 
IFdSS(JX).EQ.O)  THEN 

READ(16,3104)(TEXT(JXX),JXX=1,25) 

WRITE(23,»)  (TEXT(JXX),JXX=1,25) 
ELSE 

READ(16,3104)  (TEXT(JXX), JXX=1,25) 

WRITE f 23, *) (TEXT (JXX),JXX=1, 25) 

READ(16,3104)  (TEXT(JXX) , JXX=1,25) 

WRITE(23,*)(TEXT(JXX),JXX=1,25) 

WRITE(23,*) 
ENDIF 
ENDIF 
REWIND  16 

7  CONTINUE 

3104  FORMAT  (25A) 

RETURN 
END 

C      .fc**¥*****i*.**Y***********1i*********  *********************  **** 
C     *  * 

C    *         SUBROUTINE  PENET   (PENET, FOR)  * 

C    *  * 

C    ************************************************************ 

SUBROUTINE  PENET (CLR, P0S2, IFLAGP,STRFA,STRFB,A,B,EP,YP,HH,TOP 
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l.NNN) 
CHARACTER*!  POS2U0)  ,CLR(4) ,  TOP (9)  -  ST 
IFLAGP=0 

WRITE**,*)  CLR,P0S2 
57   WRITE**, *)  '  THIS  SECTION  CHECKS  PENETRATION  ' 
WRITE  (*,*) 
WRITE(*t*) 

WRITE(*,t)  '  IS  THE  REINFORCEMENT  MADE  FROM  THE  SAME' 
WRITE(*, ' (A\J ' ) '  MATERIAL  AS  THE  PLATE?<Y)' 
READ  (*,3i0)  ST 
310   FORMAT (A) 

IF  (ST.EQ.'Y'.OR.ST.EQ.1  '.OR.ST.EQ.'Y')  THEN 
WRITE**,*)  CLR.P0S2 
CALL  RETURN (CLR) 
53   WRITE**, *)  CLR,POS2 

WRITE**, *)'  DO  YOU  WANT  TO:' 
WRITE (*,*) 
WRITE**, *) 

WRITE**, *?'  (i)  ENTER  THE  ELASTIC  MODULUS  AND' 
MRIIEt*.*)1     YIELD  STRENGTH  OF  THE  MATERIAL' 
WRITE**-*)'  (2)  RETRIEVE  THOSE  VALUES  FROM  THE' 
WRITER*)'     ACCEPTABLE  MATERIALS  DATA  BASE' 
WRITE**,*) 
WRITE(*,*i 

WRITE**, '(Ay5)'  INPUT  YOUR  CHOICE — >' 
READ  (*,"BN:I6)')I222 
IF  CI222.EQ,1)  THEN 

WRITE**,*)  CLR.P0S2 

WRITE**,*)'  ENTER  THE  ELASTIC  MODULUS  IN  KSI' 

WRITE**, *> 

WRITER,*)'  FOR  EXAMPLE,  THE  ELASTIC  MODULUS  OF  A-36  STEEL' 

WRITE**,*)'  WOULD  BE  ENTERED  AS  29000.' 

WRITE**,*) 

READ  (*.*/  ER 

EP=ER 

WRITE**,*)  CLR,P052 

WRITE'*,*) 

WRITE**, *)'  ENTER  THE  YIELD  STREN6TH  IN  KSI' 

WRITE**,*) 

WRITE?*,*)'  FOR  A-36,  YOU  WOULD  ENTER  36.' 

WRITE**,*) 

READ**,*)  YR 

yp=VR 

WRITE  (*,*) 

ELSEIF(I222.EQ.2)  THEN 

CALL  PETRI (CLR, P0S2, YP, EP, PR, DEN, IFLA6M, IFLAMI , TOP) 

WRITE**,*) 

CALL  RE  TURN (CLR) 

YR=YP 

ER=EP 

ELSE 

WRITE-'*,*)'  YOU  ONLY  HAVE  TWO  CHOICES' 

CALL  RETURN (CLR) 

SO  TO  53 

END  IF 
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ELSE 

WRITE  I*.*)  CLR,P0S2 

WRITE!*,*) '  MATERIAL  OF  REINFORCEMENT  IS  DIFFERENT  THAN  ' 
WRITE!*,*)'  THE  MATERIAL  OF  THE  PLATE.' 
WRITER,*) 

WRITE!*,*)'  YOU  MUST  ENTER  THE  ELASTIC  MODULUS  AND  YIELD  ' 
WRITE!*,*)'  STRENGTH  OF  EACH  SEPARATELY  OR  RETRIEVE  THEM  FROM  ' 
WRITE!*,*)'  ACCEPTABLE  MATERIALS  DATA  BASE,' 
WRITE!*,*) 
CALL  RETURN (CLR) 
54   WRITE!*,*)  CLR.P0S2 

WRITE!*,*)'  DO  YOU  WANT  TO;' 
WRITE!*-,*) 
WRITE!*,*) 

WRITE!*,*)'   (1)  ENTER  THE  ELASTIC  MODULUS  AND  YIELD' 
WRITE!*,*)'     STRENGTH  OF  REINFORCEMENT  AND  PLATE?' 
WRITE!*,*)'   (2)  RETRIEVE  THOSE  VALUES  FROM  THE' 
WRITE!*,*)'      ACCEPTABLE  MATERIALS  DATA  BASE?' 
WRITE!*,*) 
WRITE!*,*) 

WRITE!*,' (A\)')'  INPUT  YOUR  CHOICE — >' 
READ  {*,'(BN,I6)'5I22 
IF  (I22.EQ.1)  THEN 
WRITE!*,*)  CLR.P0S2 

WRITE!*, ♦)'  ENTER  THE  ELASTIC  MODULUS  OF  THE  REINFORCEMENT' 
WRITE!*,*)'  IN  KSI.  FOR  EXAMPLE,  THE  ELASTIC  MODULUS' 
WRITE!*,*)'  FOR  A-36  STEEL  WOULD  BE  ENTERED  AS  29000.' 
READ!*,*)  ER 
WRITE!*, *)CLR,P0S2 
WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  YIELD  STRENGTH  OF  THE  REINFORCEMENT' 
WRITE!*,*)'  IN  KSI.  FOR  A-36  STEEL  IT  WOULD  BE  36.' 
READ!*,*)  YR 
WRITE!*,*)  CLR,P052 

WRITE!*,*)'  ENTER  THE  ELASTIC  MODULUS  OF  THE  REINFORCEMENT' 
READ  (*,*)  EP 
WRITE!*,*)  CLR,P0S2 

WRITE!*,*)'  ENTER  THE  YIELD  STRENGTH  OF  THE  PLATE' 
READ!*,*)  YP 
WRITE!*, *)CLR,P0S2 

ELSEIF  (I22.EQ.2)  THEN 

WRITE!*,*)  CLR.P0S2 

WRITE!*,*)'  YOU  ARE  TO  RETRIEVE  THE  ELASTIC  MODULUS' 

WRITE!*,*)'  AND  YIELD  STRENGTH  OF  THE  REINFORCEMENT' 

WRITE!*,*)'  MATERIAL.' 

WRITE!*,*) 

WRITE!*,*) 

CALL  RETURN (CLR) 

CALL  RETRI (CLR, P0S2, YP, EP, PR, DEN, IFLAGM, IFLAMI , TOP) 

ER=EP 

YR=YP 

WRITE!*,*) 

CALL  RETURN (CLR) 

WRITE!*, *)CLR,P0S2 

WRITE!*,*)'  YOU  ARE  TO  RETRIEVE  THE  ELASTIC  MODULUS' 
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WRITE!*,*)'  AND  YIELD  STRENGTH  OF  THE  PLATE' 

WRITE (*,*) 

WRITE!*,*) 

CALL  RETURN (CLR) 

CALL  RETRI(CLR,PQS2,YP,EP,PR,DEN!IFLAGM,IFLAMI,TGP) 

WRIT£(*, *) 

CALL  RETURN  (CLR) 

ENDIF 
ENDIF 

WRITE!*,*)  CLR, TOP 

WRITE!*,*)'  THE  ELASTIC  HODULOUS  IN  KSI  OF  THE  REINFORCEMENT  IS:1 
HRITE(*i*)ER 
WRITE!*, *) 

WRITE!*,*)'  THE  YIELD  STRENGTH  IN  KSI  OF  THE  REINFORCEMENT  IS:' 
WRITE!*, *)YR 
WRITE!*,*) 

WRITE!*,*)'  THE  ELASTIC  MODULOUS  IN  KSI  OF  THE  PLATE  IS:' 
WRITE!*, *)EP 
WRITE!*,*) 

WRITE!*,*)'  THE  YIELD  STRENGTH  IN  KSI  OF  THE  PLATE  IS:' 
WRITE!*, *)YP 
WRITE!*,*) 
WRITE!*,*) 
WRITE!*,*) 
WRITE!*,*) 

WRITE!*,' !A\)')'  ARE  THESE  VALUES  CORRECT? W 
READ  (*,310)  ST 

IF  (ST.EQ.'Y'.OR.  ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 
GO  TO  56 
ELSE 

WRITE!*,*)  CLR.P0S2 

WRITE!*,*)  '  YOU  SHOULD  REENTER  THE  VALUES' 
WRITE!*,*! 
CALL  RETURN (CLR) 
GO  TO  57 
ENDIF 
56   CONTINUE 

WRITE!*,*!  CLRtTOP 

WRITE!*,*)'     YOU  MUST  ENTER  IN  INCHES  THE  FOLLOWING  ' 

WRITE!*, *)'  DIMENSIONS  WHICH  ARE  REFERENCED  TO  FIGURE  5.1:' 

WRITE!*,*) 

WRITE!*,*) 

WRITE!*,*)  '  ENTER  DIMENSION  A,  WIDTH  OF  PLATE  CONTAINING' 

WRITE!*,*)  '  PENETRATION' 

READ(*,*)A 

WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  DIMENSION  6,  LENGTH  OF  PLATE' 

WRITE!*,*)5  CONTAINING  PENETRATION' 

READ!*,*)  B 

WRITE!*,*) 

WRITE!*,*)'  ENTER  THE  DIMENSION  DO,  OUTSIDE  DIAMETER  OF' 

WRITE!*,*)'  REINFORCEMENT' 

READ i*-*) DO 

WRITE!*, ♦) 

WRITE!*,*)'  ENTER  THE  DIEMNSION  T,  WIDTH  OF  CIRCUMFERENTIAL' 
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WRITE(*.*>'  REINFORCEMENT' 
READ  **,*)  T 
WRITE**,*) 

WRITE(*,*)'  ENTER  THE  DIMENSION  H,  HEIGHT  OF  REINFORCEMENT' 
READ  (*,*)  H 
WRITE(*-*) 

WRITE**,*)'  ENTER  THE  DIMENSION  H,  THICKNESS  OF  PLATE' 
READ(*,*)HH 
WRITE(*,*! 

STFFA=1 . - (DO/A)  +  (ER/EF) * (2. * (T/A) * ( *H/HH) **3. ) ) 
STFFB=1 . - (DO/B) + (ER/EP) * (2. * (T/B) * ( (H/HH) **3. ) ) 
STRFA=i . - (DO/A) + ( YR/YP) * (2. * (T/A) + ( (H/HH) **2. ) ) 
STRFB=1.-(00/B)  +  (YR/YP)*(2.»(T/B)*((H/HH)**2J) 
IFLAGE=0 

WRITE!*,*)  CLR.P0S2 

WRITE (*,*)  '  THE  CALCULATED  VALUES  ARE:' 
511   WRITE(*,*) 
WRITE!*,*) 

WRITE(*,*)'  STIFFNESS  FACTORS: ' 
WRITE**,*) 

WRITE(*,*)'  STFFA=  ',STFFA,'    STFFB=  '.STFFB 
WRITE (*,*) 
WRITE(*,*) 

WRITE**,*)3  STRENGTH  REDUCTION  FACTORS:' 
WRITER,  *) 

WRITE(*,*)'  STRFA=  SSTRFA,'    STRFB=  ',STRFB 
WRITE**,*) 
WRITE**,*) 
CALL  RETURN (CLR) 
IF  (IFLAGE.EQ.l)  GO  TO  510 
IFLAGP=1 
SGA=STRFA 
SGB=STRFB 

IF  (STRFA.LT.1..0R.STRFB.LL1.)  THEN 
WRITE(*,*)  CLR,P0S2 

WRITE**,*)'  IS  THE  PENETRATION  NEAR  A  CENTER' 
WRITE**,' (A\)')'  LINE  OF  THE  PANEL? <Y>' 
READ  (*,310)  ST 
IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  »)  THEN 
If (STRFA.LT. 1.)  THEN 
WRITE**, *!CLR,P0S2 

WRITE**,*)'     IS  THE  CENTER  OF  THE  PENETRATION  CLOSE  TO' 
WRITE**,*!'  TO  THE  CENTER  LINE  THAT  IS  PARALLEL  TO' 
WRITE**,' (A\)')'  SIDE  A  ?<Y>' 
READ(*,310)ST 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 
CONTINUE 
ELSE 
STRFA=1, 
END  IF 
ELSE 
ENDIF 
IF*STRFB.LT.l.)  THEN 
WRITE**, *)CLR,P0S2 
WRITE**,*)'    IS  THE  CENTER  OF  THE  PENETRATION  CLOSE  TO' 
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WRITER*)'  THE  CENTER  LINE  THAT  IS  PARALLEL  TO' 

WRITE(*,'(A\)')'  SIDE  B  ?<Y>' 

READ(*t310)ST 
IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 

CONTINUE 

ELSE 

STRFB=1. 

ENDIF 
ELSE 
ENDIF 

ELSE 

STRFA=1. 

STRFB=i. 

ENDIF 

ELSE 

STRFA=1. 

STRFB=1. 

ENDIF 
CONTINUE 

WRITE (*,*)CLR,POS2 

WRITE'*,*!'  THE  STIFFNESS  AND  STRENGTH  FACTORS  ARE:' 
WRITE(*,*) 
WRITE(*t*) 
IFLAGE=1 
60  TO  511 
CONTINUE 
WRITE(23,*> 


WRITE(23, 
WRITE (23, 
WRITE (23, 
WRITE (23, 
WRITE (23. 
WRITE (23, 
NRITE(23i 
WRITE (23, 
WRITE (23, 
WRITE  (23, 
WRITE (23, 
WRITE (23, 
WRITE(23, 
WRITE (23- 
WRITE (23, 
WRITE(23, 
WRITE (23, 
WRITE (23, 
IF'STFFA. 

«te;23, 

WRITE i 23, 
WRITE (23, 
WRITE (23, 
WRITE (23, 
WRITE '23, 
ELSE 

WRITE (23, 
WRITE (23, 


) '  THE  PLATE  THAT  HAS  THE  PENETRATION  IS: ' 

)NNN 

) 

)'  ITS  DIMENSIONS  ARE:' 

)'  A:  ',A,'   B:  ',B 

)'  THICKNESS:  ',HH 

)'  OUTSIDE  DIAMETER  OF  REINFORCEMENT,  DO:  ',D0 

)'  WIDTH  OF  CIRCUMFERENTIAL  REINFORCEMENT,  T:  ',T 

)'  HEIGHT  OF  REINFORCEMENT,  H:  ',H 


) '  THE  RESULTS  OF  THE  CHECK  FOR  PENETRATION1 

) '  REDUCTION  ARE  THE  FOLLOWING: ' 

) 

T ,  1  ■ )  THEN 

)■■   THE  STIFFNESS  FACTOR  IK   THE  X  DIRECTION  IS 

)'  LESS  THAN  1.  ALTHOUGH  THIS  IS  NOT  CAUSE' 

»»  FOR  REJECTION,  THERE  IS  A  CONCERN.' 

j 

)'  THE  STIFFNESS  FACTOR  SFFA  IS." 

JSTFFA 


THE  STIFFNESS  FACTOR  IN  THE  X  DIRECTION 
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WRITE (23, *>STFFA 

ENDIF 

IFtSTFFB.LT.l.)  THEN 

WRITE(23,*! 

WRITE(23,*>'  THE  STIFFNESS  FACTOR  IN  THE  Y  DIRECTION  IS' 

WRITE(23,*)'  LESS  THAN  1.  ALTHOUGH  THIS  IS  NOT  CAUSE' 

WRITE(23,*)!  FOR  REJECTION,  THERE  IS  A  CONCERN.' 

WRITE  (23,0 

WRITE (23,*)'  THE  STIFFNESS  FACTOR  STFFB  IS:' 

WRITE (23, *! STFFB 

ELSE 

WRITE (23,*! 

WRITE(23,*)'  THE  STIFFNESS  FACTOR  IN  THE  Y  DIRECTION  IS' 

WRITE (23,*) STFFB 

ENDIF 

WRITE(23,*) 

WRITE(23,*)'  THE  CALCULATED  STRENGTH  REDUCTION  FACTORS  ARE:' 

WRITE (23,*) 

WRITEC23,*)'  STRFA=  ',SGA,'      STRFB=  ',SGB 

WRITE(23,*) 

WRITE (23,*) 

WRITE(23,*)'  AFTER  CONSIDERING  THE  LOCATION  OF  THE' 

WRITE(23,*)'  PENETRATION.  THE  STRENGTH  REDUCTION' 

WRITE(23,*)'  FACTORS  THAT  SHOULD  BE  USED  IN  FUTURE' 

WRITE(23,*)'  CALCULATIONS  ARE:' 

WRITE(23,*)'  STRFA=  '.STRFft,'    STRFB=  '.STRFB 

NNN=NNN+1 

RETURN 

END 
C    M**t**m****t*m***MM****************m*m*MiM*m** 
C    *  * 

C    *  SUBROUTINE  PRINT1  (PRINT.FOR)  * 

C    *  * 

C      *******  +  ***^*  +  ***  +  *  +  ****  +  *******  +  +  **********  +  'P******  +  ****'!*  + 

SUBROUTINE  PRINT1 ( I , ILOOP, IFLAG) 
CHARACTER*64  FLM2(33) 
DIMENSION  6(5,3) 
DO  30  1C0UNT=1,3 
IS=ICOUNT+ILOOP 
READ  (li,'(A64)')FLM2(IS) 
IF  (ICOUNT.EQ.l  )  THEN 

READ  (10,,(F5.0.F6.0iF6.0)')6(I.l)iG(I»2).6(I.3) 
WRITE  (*,' (IX, A64,F6.4)')FLM2(IS), 6(1, IFLAG) 
ELSE 

WRITE  (*,'(1X,A64)')FLM2(IS) 
ENDIF 
30    CONTINUE 

IL00P=IL00P+3 

RETURN 

END 

c   *M^Mm************************************************** 

c    *  * 

C    *       SUBROUTINE  RETURN   (RETURN. FOR)  * 

C    *  * 

C      ♦*****♦**  M^MM**********************************  ********* 
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SUBROUTINE  RETURN (CLR) 

CHARACTER*!  CLR (4), ST 
1340      WRITE(*,'(A\)')'  TO  CONTINUE,  HIT  RETURN' 

READ  (*,350)  ST 

IF  iST.EQ.'   '!  THEN 

WRITE  (*,*)  CLR 

ELSE 

WRITE  (*,*)'  YOU  MUST  HIT  RETURN  TO  CONTINUE' 

GO  TO  1340 

END  IF 
350        FORMAT (A) 

END 

C         *  * 

C  *  SUBROUTINE  STREN  (S.FOR)  * 

C  *  * 

C  fr*********** ************************************************ 

SUBROUTINE  STREN (ID, CLR, P0S2, IFLA6X, TOP, VOLUME, IRUG.PAN, IT) 
CHARACTER*!  CLR (4) , P0S2 ( 10) , TOP (9) , ST 
CHARACTER*64  ID 
DIMENSION  PAN (20, 3) 
310      FORMAT (A) 

WRITE  (♦,*)  CLR,P0S2 
CALL  RETURN (CLR) 
51        WRITE(*,»)  CLR,P0S2 

WRITE  (*,*)'  DO  YOU  WANT  TO:' 

WRITE(*,*> 

WRITE(*,*) 

WRITE**,*)'     (1)     CHECK  PENETRATIONS  FOR  STRENGTH' 

WRITEi*,*)'     (2)     CHECK  COVER  FOR  STRENGTH' 

WRITE(*,*)'     (33     CHECK  BODY  PANELS  FOR  STRENGTH' 

WRITE(*,*)'     (4!     CHECK  WINDOWS  AND  LENSES  FOR  STRENGTH' 

WRITE(*,*)'     '5)     END  THIS  SECTION' 

WRITE(*,*) 

WRITE(*-.*) 

WRITE!*,' (A\)'P   INPUT  YOUR  CHOICE — >' 

READ(*,'(BN,I6)')I22 

IF(I22.EQ.!)  THEN 

WRITE (*,*)  CLR.P0S2 

WRITE!*,' (A\)')'  DO  YOU  WANT  TO  VIEW  INFORMATION? (Y)' 

READ  (*,310)  ST 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'   ')  THEN 

GO  TO  506 

ELSE 

GO  TO  505 

ENDIF 
506  IG0=6 

NUMB=20 

CALL  DISPLAY (IGO, NUMB, CLR) 

IG0=7 

NUMB-20 

CALL  DISPLAY (IGO, NUMB, CLR) 
505  CONTINUE 

CALL  PENFT<CLR,P0S2;IFLAGP,STRFA,STRF8,A,B,EP,YP,HH,T0P 
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l,f 

ELSEIF(I22.EQ.2)  THEN 
11=5 

write;*,*)  clr,pos2 

write!*-*)  '   this  section  screens  the  cover  and  bolts' 

write(*,*)'  if  present  for  strength' 

WRITE(*,*) 

CALL  RETURN (CLR) 

CALL  COVER ( ID, TOP, IFLAGM, CLR, P0S2) 

ELSEIF  (I22.EQ.3)  THEN 

IFLAGV=0 

IFLAGP=0 

IFLAGQ=0 

IFLA6Z=0 

CALL  BODY (CLR, P0S2, I FLAGP, VOLUME, IRUG,PAN, IT,TOP) 

ELSEIF ( 122. EQ. 4)  THEN 

WRITE!*, *)CLR,P0S2 

WRITE!*,*)'  YOU  HAVE  CHOSEN  TO  CHECK  WINDOWS.' 

WRITE(*,*) 

WRITE(*,*) 

CALL  RETURN (CLR) 

CALL  WINDOW (ID) 

ELSEIF  (122. EG. 5)  THEN 

GO  TO  50 

ELSE 

WRITEt*!*)'  YOU  ONLY  HAVE  FIVE  CHOICES' 

CALL  RETURN (CLR) 

GO  TO  51 

ENDIF 
GO  TO  51 
50     WRITE!*, *)CLR,P0S2 
RETURN 
END 

c    *  * 

C    *        SUBROUTINE  TAB52   (TAB52.F0R)  * 

C    *  * 

SUBROUTINE  TAB52 ( IR, GX, XX, CLR, P0S2) 
DIMENSION  AT (5, 3) 
CHARACTERS  CLR(4),POS2(10) 
OPEN  (21,  FILE=' TAB52.DOC) 
DO  3103  IX=i,5 

3103  READ(21,3104)(AT(IX,JX),JX=1,3) 

3104  FORMAT (F5.0,F5.0,F5.0) 
CLOSE (21) 

201  CONTINUE 
GX=1./GX 
IF  (GX.GE. 1.0.AND.GX.LT. 1,5)  THEN 

IC=1 
ELSEIF (GX.GE. 1.5. AND.6X.LT. 2. 5)  THEN 

IC=2 
ELSEIF (GX.GE. 2. 5)  THEN 

IC=3 
ELSE 
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WRITE(*,*)CLR,PGS2 

WRITE(*,*)'  A  MUST  BE  SHORTER  THAN  B. ' 
WRITE**,*) 

WRITE'*,*)5  YOU  MUST  ENTER  A  FOLLOWED  BY  B.' 
READ  (*,*)  A.B 
GO  TO  201 
ENDIF 
200   CONTINUE 

IF  (IC.EQ.l)  THEN 

XX=AT(IR:IC) 
ELSEIF  (IC.EB.2)  THEN 

XX=AT(IR,IC) 
ELSEIF (IC.EQ.; 3)  THEN 

XX=AT(IR,IC) 
ELSEIF (IC.EQ. 4)  THEN 

XX=ATQR,IC) 
ELSEIF (IC.EQ. 5)  THEN 

XX=AT(IR,IC) 
ELSE 

WRITE(*,*>'  YOU  MUST  ENTER  A  NUMBER  1-5  CORRESPONDING' 
WRITE(*.*)'  TO  A  BOUNDARY  CONDITION  IN  TABLE  5,2' 
WRITE**,*) 

WRITE**,*)'  ENTER  A  NUMBER  1-5.' 
READ  (*,*)IC 
GO  TO  200 
ENDIF 
RETURN 
END 
C    ************ ***************************#******»************* 
C    *  * 

C    *         SUBROUTINE  SIMUL   (TESTi.FOR)  * 

C    *  * 

SUBROUTINE  SIMUL  (AA,X) 
INTEGER  FLAG 

DIMENSION  AA(9,9),X(50),A<9,10) 
N=9 

ITMAX=  200 
EPS=. 0000901 
DO  1000  M1=1,N 
1000  WRITE  (*,6666)  (AA(M,N1) ,N1=1,N) 
6666  FORMAT  (9E10.4) 
DO  42  11=1, N 
DO  41  JJ=1,N 

41  A(II,JJ)=AA(II,JJ) 

42  CONTINUE 

DO  40  11=1, N 
40  A(II,10)=X(II) 

NP1=N+1 

DO  3  1=1,  N 

ASTAR=A(I,I) 

DO  3  J=1,NP1 
3   A(I,J)=A(I,J)/ASTAR 

DO  9  ITER=1,ITMAX 

FLAG=1 
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DO  7  1=1, N 

XSTAR=X(I) 

X(I)=A(I,NP1) 

DO  5  J=1,N 

IF  (I.  EQ.  J)  GO  TO  5 

X(I)  =  X(I)  -A(I,J)*X(J) 
CONTINUE 
IF  (ABS(XSTAR-Xd))  .LE.  EPS  )  60  TO  7 

FLAG  =0 
CONTINUE 

IF  (FLA6.NE.1  )  GO  TO  9 
WRITE(*,*)'  METHOD  DID  CONVERGE' 
WRITE(*,*)  ITER,  (X(I),I=1,N) 
GO  TO  1 
CONTINUE 

WRITE  (*,*)  '  METHOD  DID  NOT  CONVERGE)' 
CONTINUE 

WRITE(*,*)  ITER,  (X(I),I=1,N) 
RETURN 
END 

*  * 

*  SUBROUTINE  V0LUM1   (VOLUME. FOR)         * 

*  * 

SUBROUTINE  V0LUM1 (CLR, P0S2, VOLUME) 
DIMENSION  F( 12), THICK (6) 
CHARACTER*!  POS2(10),CLR(4) 
WRITE  (*,*)  CLR,P0S2 


WRITE(*, 
WRITE(*, 
WRITE(*, 
WRITE(*, 
WRITE(», 
WRITE(*. 
WRITE(*, 
WRITE(», 
WRITER, 
WRITE(», 
WRITE(», 
WRITE (*, 


ENTER  THE  OUTSIDE  DIMENSIONS  AND  THICKNESS' 

OF  EACH  FACE. ' 

THE  PROGRAM  WILL  SUBTRACT  2  TIMES  THE  THICKNESS' 

OF  THE  FRONT  FACE,  2  TIMES  THE  THICKNESS  OF  THE' 

LEFT  FACE  AND  2  TIMES  THE  THICKNESS  OF  THE' 

TOP  FACE  IN  ORDER  TO  CALCULATE  THE  INTERNAL' 

VOLUME' 


IF  THIS  IS  NOT  CORRECT  FOR  YOUR  ENCLOSURE, ' 

YOU  MUST' 

ENTER  THE  INTERNAL  VOLUME  AS  YOU  MEASURED  IT' 

FROM  THE  ENCLOSURE' 
CALL  RETURN (CLR) 
WRITE  (*,*)  CLR,P0S2 

WRITE  (*,*)'  ENTER  THE  DIMENSIONS  OF  THE  FRONT  FACE' 
WRITE  (*,*)'  LENGTH  FIRST,  FOLLOWED  BY  SPACE  AND  THEN  WIDTH' 
READ  (♦,*)  F(1),F(2) 

WRITE  (*,*)'  ENTER  THE  THICKNESS  OF  THE  FRONT  FACE' 
READ  (*,*)  THICK (1) 
WRITE  (*,*)CLR,P0S2 

WRITE  (*,*)'  ENTER  THE  DIMENSIONS  OF  THE  BACK  FACE' 
WRITER, *)'  LENGTH  FIRST,  FOLLOWED  BY  SPACE  AND  THEN  WIDTH' 
READ  (*,*)  F(3),  F(4) 

WRITER,*)'  ENTER  THE  THICKNESS  OF  THE  BACK  FACE' 
READ  (*,*)  THICK (2) 
WRITE  !*,*)  CLR,  P0S2 
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WRITE !*,*)'  ENTER  THE  DIMENSIONS  OF  THE  LEFT  FACE' 

WRITE!*,*)'  LEN8TH  FIRST,  FOLLOWED  BY  SPACE  AND  THEN  WIDTH' 

READ  !*,*)  F(5),F!6) 

WRITE!*,*)'  ENTER  THE  THICKNESS  OF  THE  LEFT  FACE' 

READ!*,*)  THICK (3) 

WRITE  (*,*)  CLR.P0S2 

WRITE  (*.*)'  ENTER  THE  DIMENSIONS  OF  THE  RIGHT  FACE' 

WRITE  (*,*)'  LENGTH  FIRST,  FOLLOWED  BY  SPACE  AND  THEN  WIDTH' 

READ  !*■*)  F(7),F(8) 

WRITE!*,*)'  ENTER  THE  THICKNESS  OF  THE  RIGHT  FACE' 

READ  (*,*)  THICK (45 

WRITE!*, *)CLR,  P0S2 

WRITE!*, *)'  ENTER  THE  DIMENSIONS  OF  THE  TOP  FACE' 

WRITER,*}'  LENGTH  FIRST,  FOLLOWED  BY  SPACE  AND  THEN  WIDTH' 

READ  (*,*)  F(9),F(10) 

WRITE  (*,*)'  ENTER  THE  THICKNESS  OF  THE  TOP  FACE' 

READ!*,*)  THICK (5) 

WRITE  (*,*)  CLR,  P0S2 

WRITE  (*,*)'  ENTER  THE  DIMENSIONS  OF  THE  BOTTOM  FACE' 

WRITE!*, *)'  LENGTH  FIRST,  FOLLOWED  BY  SPACE  AND  THEN  WIDTH' 

READ  (*,*)  F(li),F(12) 

WRITE  (♦,*)  '  ENTER  THE  THICKNESS  OF  THE  BOTTOM  FACE' 

READ  (*,*)  THICK (61 

WRITE(*,*)CLR,P0S2 

WRITE  !*,*)'  FRONT  FACE: 

WRITE!*,*)  ' 

WRITE!*,*)  ' 

WRITE  (*,*)'  BACK  FACE: 

5 


LENGTH 

WIDTH 

THICKNESS- 
LENGTH 

WIDTH 

THICKNESS- 
LENGTH 

WIDTH 

THICKNESS- 
LENGTH 

WIDTH- 

THICKNESS- 
LENGTH 

WIDTH 

THICKNESS- 
LENGTH 

WIDTH 

THICKNESS- 


WRITE?*,*) 

WRITE!*,*)  ' 

WRITE!*,*)  '  LEFT  FACE: 

WRITE!*,*)  ' 

WRITE!*,*)  ' 

WRITE!*,*)  '  RIGHT  FACE." 

WRITE!*,*)  ' 

WRITE!*,*)  ' 

WRITE!*,*)  '  TOP  FACE: 

WRITE!*,*)  ' 

WRITE!*,*)  ' 

WRITE!*,*)  '  BOTTOM  FACE: 

WRITE!*,*)  ' 

WRITE!*,*)  ' 

CALL  RETURN !CLR) 

WRITE  (*,*)  CLR,P0S2 

IF((F(1).EQ.F(9)).0R.(F(2).EQ.F(9)))  THEN 

VOLUME= (F !2) -2. *THICK !3> ) * (F ( 1 ) -2. *THICK (5) ) * (F  < 10) -2. *THICK  ( 1 ) ) 

ELSEIF  !(F(1).EQ.F(10)).OR.(F(2).EQ.F(10)))  THEN 

V0LUME=(F(2)-2.*THICK(3))*(F(1)-2.*THICK!5!)*(F(9)-2.*THICK(1)) 

ENDIF 

WRITE  !*,*)'  VOLUME=  '.VOLUME 

RETURN 

END 


'tF(l) 

',F<2) 

', THICK (1) 

',F(3) 

',F(4) 

'•THICK (2) 

',F!5) 

',F(6) 

', THICK (3) 

',F(7> 

',F(8) 

'.THICK (4) 

',F(9) 

',F!10) 

', THICK (5) 

',F(1D 

',F(12) 

', THICK (6) 
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C     M**MMMMMm****M*MMMMmMMMMt***************t 

C    *  * 

C    *        SUBROUTINE  WINDOW   (W.FOR)  * 

C    *  * 

SUBROUTINE  WINDOW (ID) 

CHARACTER*!  CLR 14) ,POS2(10) ,ESC,CR,ST 

CHARACTERS  ID 

DATA  CLR/'   VCV2VJ'/ 

data  POS2/5  vtvivfvsviVHV  vevkv 

ESC=CHAR(27) 

CLR(1)=ESC 

FOS2(l)=ESC 

P0S2(8)=ESC 

CR=CHAR(13) 

ST=CR 
NN=1 
2222  WRITE(*,*)CLR,P0S2 

WRITE(*t*)»  THIS  SECTION  CHECKS  THE  STRENGTH  OF  WINDOWS' 

WRITE(*,*)'  AND  LENS  FOR  A  STATIC  PRESSURE  OF  150  PSIG. ' 

WRITE'*,*) 

WRITE(*,*>'  IN  ADDITION,  THE  LIP  THAT  HOLDS  THE  WINDOW  OR  ' 

WRITE**,*)'  LENSE  IS  ALSO  EVALUATED.' 

WRITE(*,*> 

WRITE(*,*) 

WRITE(*,*) 

CALL  RETURN(CLR) 

WRITE(*,*)CLR,P0S2 

WRITE(*,*)'  30  CFR  PART  18  REQUIRES  THAT  ALL  HEADLIGHT  LENSES  ' 

WRITE (*,*)'  BE  THE  EQUIVALENT  OF  1/2  INCH  THICK  PYREX  GLASS.' 

WRITE(*,*) 

WRITE(*,*)'  DO  ALL  GLASS  HEADLIGHT  LENSES  MEET  THIS' 

WRITE(*, ' (A\) ' ) '  REQUIREMENT? <Y> ' 

READ(*,310)ST 

WRITE (*,*) CLR, P0S2 

WRITE(23,*) 

WRITEC23,*)'  RESULTS  CHECK  FOR  STRENGTH  OF  WINDOWS  AND  LENSES.' 

WRITE(23,») 

WRITE (23,*) 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 

WRITE(23,*)'  HEADLIGHT  LENSES  COMPLY  WITH  30  CFR' 

WRITE (23,*)'  REQUIREMENTS.' 

ELSE 

WRITE (23,*)'  HEADLIGHT  LENSES  DO  NOT  COMPLY  WITH  30  CFR' 

WRITE (23,*)'  REQUIREMENTS.' 

ENDIF 

WRITE(*,*)CLR,P0S2 

WRITE(*,'(A\)')'  DO  YOU  WANT  TO  VIEW  INFORMATION? <Y>' 

READ (*, 310) ST 
310   FORMAT (A) 

WRITE(*,*)  CLR.PQS2 

IF  tST.EQ.'Y'.CR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 

180=17 

NUMB=20 

CALL  DISPLAY (IGO. NUMB, CLR) 
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IG0=18 
NUMB=20 

CALL  DISPLAY UGO, NUMB, CLR) 
END  IF 

WRITE (»,*)CLR,P0S2 
WRITE(*t*) '  IS  THE  WINDOW:' 
WRITE(*1*) 
WRITE?*,*) 

WRITE?*,*)'    (15  RECTANGULAR' 
WRITE?*,*)'    (2)  CIRCULAR' 
WRITE?*,*) 
WRITE?*,*) 

WRITE?*, '(AM')'  INPUT  YOUR  CHOICE — )' 
READ(*,'?BN,I6)')IB 
2035  WRITE?*,*)  CLR,P0S2 
IF(IB.EQ.l)  THEN 
WRITE?*,*)'  ENTER  THE  DIMENSIONS  OF  THE  RECTANGULAR  WINDOW  IN' 
WRITE?*,*)'  INCHES.  A  SHOULD  BE  ENTERED  FOLLOWED  BY  A  SPACE  AND' 
WRITE?*,*) '  THEN  ENTER  B.  DIMENSION  A  MUST  BE  LARGER  THAN  ' 
WRITE?*,*)'  DIMENSION  B. ' 
READ?*,*)  A,B 
G=A/B 

IF(B.LT.l.)  THEN 
GO  TO  2035 
ELSE 
ENDIF 
PB=150.*B 

WRITE?*, *)CLR,P0S2 

WRITE?*,*)'  ENTER  THE  PLATE  THICKNESS  IN  INCHES.' 
READ?*,*)T 
WRITE?*, *)CLR,P0S2 

WRITE?*,*)'  ENTER  THE  CONSTANT  Kl  FOR  THE  CURVE  SHOWN  IN' 
WRITE?*,*)'  FIGURE  5.   USE  CURVE  A  WHEN  EDGES  ARE' 
WRITE?*,*)'  SIMPLY  SUPPORTED;  USE  CURVE  B  WHEN  EDGES  ARE' 
WRITE?*,*)'  CLAMPED.' 
WRITE?*,*) 

WRITE?*,*)'  THE  VALUE  FOR  RATIO  A/B  IS' 
WRITE?*, *)G 
WRITE!*,*) 

WRITE?*,' (A\)')'  ENTER  VALUE  FROM  GRAPH. >' 

READ?*,*)GK1 

SI6MAX=GK1*((B/T)**2.)*150. 

IBX=1 

WRITE?*, *)CLR,P0S2 

WRITE?*,*)'  VALUE  CALCULATED  FOR  SIGMAX  IS  '.SIGMAX 

ELSEIF  (IB.EQ.2)  THEN 

IBX=2 

WRITE''*,*)  CLR,P0S2 

WRITE?*,*)'  ARE  THE  EDGES:' 

WRITE?*,*) 

WRITE?*,*)'  (U  SIMPLY  SUPPORTED' 

WRITE?*,*)'  (2)  CLAMPED' 

WRITE?*,*) 

WRITE?*,*) 

WRITE?*, '(A\)')'  ENTER  YOUR  CHOICE )' 


so 


READ (*, 310) ST 
WRITER, *)  CLR,P0S2 

IF(ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EG.'  ')  THEN 

GK2=.3025 

ELSE 

GK2=. 1875 

END  IF 

WRITE!*,*)'  ENTER  THE  UNSUPPORTED  DIAMETER  OF  THE  WINDOW 
WRITE(*,*)'  IN  INCHES.' 
WRITE(*,*) 
READ(*,*)D 
PB=150.*D/2. 
WRITE (*,*)CLR,P0S2 

WRITE(*,*)'  ENTER  THE  THICKNESS  OF  THE  WINDOW  IN  INCHES.' 
READ(*,*)T 
WRITE(*.»)  CLR.P0S2 
SI6MX=GK2* ( (D/T) **2. ) *150. 

WRITE(*,*>'  VALUE  CALCULATED  FOR  SIGMAX  IS  '.SIGMX 
ENDIF 
WRITE(*.*) 
WRITE(*,*) 
CALL  RETURN (CLR) 
IWW=0 
2323     WRITE(*,*)  CLR.P0S2 

WRITE(*,*)'  IS  THE  WINDOW  MATERIAL:' 
WRITER,  *) 

WRITE(*,*)'    (1)  SODA  LIME  GLASS' 
WRITE(*,*)'    (2)  BOROSILICATE  GLASS' 
WRITER*)'    (3)  LEXAN  101  PLASTIC 
WRITEt*.*)'    (4)  MERLON  3113  PLASTIC 
WRI TE ( * , * ) '    (5)  NONE  OF  THE  ABOVE' 
WRITE(*,*) 

WRITE(*, ' (A\) ' ) '  ENTER  YOUR  CHOICE ) ' 

READ(*,*)IW 
WRITER  *)CLR,P0S2 
IF(IW.EQ.l)  THEN 
SIGMAX=6600. 
IWW=1 

ELSEIF(IW.EQ.2)  THEN 
SIGMAX=6300. 
IWW=1 

ELSEIF(IW.EQ.3>  THEN 
SIGMAX=9000. 
ELSEIF(IW.EQ.4)  THEN 
SIGMAX=9200. 
ELSEIF(IW.EQ.5)  THEN 

WRITE(*,*)'  ENTER  THE  YIELD  STRENGTH  IN  KSI  OF  THE' 

WRITE (*,*)'  WINDOW  MATERIAL.' 

READ(*.*)  SIGMAX 

SIGMAX=SI6MAX*1000. 

ELSE 

GO  TO  2323 

ENDIF 

IF(IWW.EQ.l)  THEN 

STEST=.5*SIGMAX 
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ELSEIFdfcW.EQ.0)  THEN 

STEST=SIGMAX 

END  IF 

WRITE (23,*) 

WRITE (23,*)  '  WINDOW  NUMBER   ',NN,'   IS  BEING  SCREENED.' 

WRITE(23,*) 

WRITE(23,*>'  THE  DIMENSIONS  OF  THE  WINDOW  ARE:' 

WRITE(23,*) 

IF  (IBX.EQ.l)  THEN 

WRITE(23,*)'  THE  WINDOW  IS  RECTANGULAR' 

WRITE(23,*)'  DIMENSION  A   ',A 

WRITE(23,*)'  DIMENSION  B   ',B 

WRITE(23,*)'  PLATE  THICKNESS   ',T 

WRITE (23,*)'  RATIO  A/B  IS' .6 

WRITE (23,*)'  VALUE  READ  FROM  GRAPH  »,6K1 

WRITE (23,*) 

WRITEC23,*)'  CALCULATED  VALUES  ARE  ' 

WRITE(23,*)'  SIGMAX=   '.SIGMAX 

WRITE (23,*! 

WIAX=. 5*150. *B 

ELSEIF  (IBX.EQ.2)  THEN 

WRITEC23,*)'  WINDOW  IS  CIRCULAR' 

WRITE(23,*)'  UNSUPPORTED  DIAMETER  OF  WINDOW  IN  INCHES  ',D 

WRITE(23,*)'  THICKNESS  OF  WINDOW  ',T 

WRITE(23,*)'  VALUE  OF  6K2  ',GK2 

WRITE(23,*) 

WRITE (23,*)'  CALCULATED  , VALUES  ARE  ' 

WRITE(23,*)'  SIGMAX=   '.SISHX 

VHAX=150.*.5*D/2. 

ENDIF 

WRITE (23,*) 

WRITE(23,*)'  YIELD  STRENGTH  OF  WINDOW  MATERIAL  IS:   ' 

WRITE (23, *)SIGMAX 

IF(STEST.GT.SISHX)  THEN 

WRITE(23,*)'  WINDOW  IS  ACCEPTABLE' 

ELSE 

WRITE(23,*)'  WINDOW  IS  UNACCEPTABLE  SINCE  SIGMAX  IS  ' 

WRITE (23,*)'  GREATER  THAN  YIELD  STRENGTH' 

ENDIF 

WRITE(*,*)CLR,P0S2 

WR!TE(*,*)'  THIS  SECTION  CHECKS  THE  LIP  FOR  STRENGTH' 
WRITE(23,*) 

WRITE(23,*)'  RESULTS  OF  CHECK  FOR  LIP  STRENGTH' 
WRITE(23,*) 
WRITE (*,*) 
WRITE(*,*> 
CALL  RETURN (CLR) 
WRITE(23,*) 

WRITE'23,*)'  YIELD  STRENGTH  OF  WINDOW  MATERIAL  ', SIGMAX 
WRITE(*,*)CLR,P0S2 

WRITEt*,*)'  FROM  FIGURE  6.   ENTER  THE  FOLLOWING:' 
WRITE(*,*) 

WRITE(*,*) 'ENTER  TL  FOLLOWED  BY  A  SPACE  AND  THEN  L  IN  INCHES' 
READ(*,*!TL,TLL 
WRITE(23,*)'  TL  AND  L  FROM  FIGURE  6.' 


s: 


WRITE<23,*)TL,TLL 

WRITEC23,*) 

WRITE(+,*)CLR,P0S2 

WRITE(*,*!'  ENTER  C,  THE  LOADING  LENGTH  IN  INCHES' 

READ(*,*)C 

WRITEC23,*)'  C,  THE  LOADING  LENGTH  IN  INCHES' 

WRITE(23,*)C 

WRITEC23,*) 

WRITE(*,*)CLR,P0S2 

WRITE(*,*)'  ENTER  THE  YIELD  STRENGTH  OF  THE  LIP  MATERIAL  IN' 

WRITE(*,*)'  KSL' 

REAB(*,*)  YS 

YS=Y5*1000. 

WRITE(23,*)'  YIELD  STRENGTH  OF  LIP  MATERIAL' 

WRITE(23,*)YS 

WRITE (23,*) 

WRITE(*,*)CLR,P0S2 

VA= (YS* (TL**2. ) ) / (4*TLL-2. *C) 

WRITE (23,*) 

write;*,*) 
write(*,*)'  va=  ',va 

WRITE(*,*) 

WRITE(*,*) 

CALL  RETURN (CLR) 

WRITE(23,*)'  CALCULATED  VALUE  OF  VA   ',VA 

WRITE(23,*)'  VHAX  FOR  THIS  WINDOW  IS  ',  WAX 

WRITE(23,») 

IF  (VMAX.LT.VA)  THEN 

WRITE(23,*)'  THE  ENCLOSURE  IS  ACCEPTABLE' 

ELSE 

WRITE(23,*)'  THE  ENCLOSURE  IS  UNACCEPTABLE' 

ENDIF 

WRITE(*,*)CLR,P0S2 

WRITE(*,*)'  DO  YOU  WANT  TO  SCREEN  ANOTHER  WINDOW?' 

WRITE(*,'(A\)')'  ENTER  YOUR  ANSW£R<Y>' 

READ(*,310)ST 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 

NN=NN+1 

GO  TO  2222 

ELSE 

ENDIF 

WRITE(*,*)CLR,POS2 

RETURN 

END 
c    ************************************************************ 
C    *  * 

C    *       SUBROUTINE  WELDCH    (WELD. FOR)  * 

C    *  * 

c    ************************************************************ 
SUBROUTINE  WELDCH ( ID, CLR, P0S2, AAA) 

DIMENSION  AAA (12, 4) 

CHARACTERS  ID 

CHARACTER*!  POS2U0)  ,CLR(4) 
393     WRITE(*,*)  CLR.P0S2 

i!  ITY' 
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WRITE(*,*> 
WRITEt*,*) 
CALL  RETURN  fCLR) 
160=2 
NUKB=2i 

CALL  DISPLAY (160. NUMB, CLR) 
160=3 
NUMB=2i 

CALL  DISPLAY (160, NUMB, CLR) 
CALL  WELDE(CLR) 
CALL  RETURN (CLR) 
WRITE(*,*>  CLR.P0S2 

WRITE**,*)'    ARE  YOU  READY  TO  EVALUATE  THE  ENCLOSURE  FOR 
I WELD' 
WRITE!*,*) 'QUALITY  AND  CLASSIFICATION?  IF  YOU  ANSWER  "N",  YOU' 
WRITE (*,'(A\)')'  WILL  BE  SHOWN  THE  PREVIOUS  SCREENS, (Y)' 
READ(*,310)ST 

IF  (ST.EQ.'Y'.OR.ST.EQ.'Y'.OR.ST.EQ.'  ')  THEN 
i30  TO  395 
ELSE 

60  TO  393 
ENDIF 
395   WRITE(*,*)  CLR.P0S2 

WR1TE<*-*)  'DOES  THE  ENCLOSURE  HAVE  ANY  CLASS  VI' 
WRITE (*,'(A\)')'  JOINTS?  (Y)' 
READ(*,310>  ST 
WRITE(*,*)CLR,P0S2 

IF  (ST.EQ.'Y'.OR.ST.EQ.'  '.OR.ST.EQ.'Y')  THEN 
WRITE(23,*5'    THE  ENCLOSURE  FAILS  THE  WELD  QUALITY  SCREEN' 
WRITE (23,*) 'BECAUSE  IT  HAS  A  CLASS  VI  WELD  WHICH  IS' 
WRITE(23,*) 'UNACCEPTABLE  FOR  DYNAMIC  LOADING. ' 
WRITE(23,*) 
ELSE 

WRITEf*,*)'    DO  THE  VISUAL  AND  DIMENSIONAL  INSPECTIONS' 
WRITE (*,*) 'OF  THE  WELDS  MEET  THE  REQUIREMENTS  IN' 
WRITE(*t ' (A\) *  j "  TABLE  4.2?<Y>' 
READ(*,3i0)  ST 

IF  (ST.EQ. 'Y!. OR.ST.EQ.'Y'. OR, ST, EQ.'  ')  THEN 
WRITE(23,*)'    THE  ENCLOSURE  FAILS  THE  WELD  QUALITY  SCREEN' 
WRITE (23,*) 'BECAUSE  OF  VISUAL  INSPECTION  OR  LACK  OF' 
WRITE (23,*)' INFORMATION  ABOUT  THE  QUALITY  OR  TYPE  OF  WELD.' 
WRITE(23,*) 
ELSE 

WRITE(23,*)'   THE  ENCLOSURE  PASSE5  THE  WELD  QUALITY  SCREEN.' 
ENDIF 
ENDIF 
310   FORMAT (A) 
END 

c   m**»*******************  *******♦**********#**  *************** 

C    *  * 

C    *        SUBROUTINE  WELDE   (WELDE.FOR)  * 

C    *  * 

SUBROUTINE  WELDE (CLR) 
CHARACTER*!  TEXT (80) , CLR (4) 


S4 


OPEN  (21, FILE=' WELDEF.DOC') 
WRITE(*,*)  CLR 
DO  3103  IX=1t21 

READ  (21.3101) (TEXT(JX),JX=1, 80) 
WRITE(*,3102HTEXT(JX),JX=1,79) 
3103  CONTINUE 
REWIND (21) 
CLOSE (21) 
RETURN 

3101  FORMAT (80A) 

3102  FORMAT (1X.79A1) 
END 
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